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Atmospheric  Transmittance/ Radiance: 
Computer  Code  LOWTRAN  4 


1.  INTRODUCTION 

This  report  describes  the  development  of  a Fortran  computer  code  LOWTRAN  4 
designed  to  calculate  atmospheric  radiance  and/or  transmittance  for  a given  at- 
mospheric path  at  moderate  spectral  resolution.  The  present  code  is  based  on  the 

* 1 

current  LOWTRAN  atmospheric  transmittance  code,  LOWTRAN  3B  (and  its  pre- 
decessors LOWTRAN  3,  2 LOWTRAN  2,  3 LOWTRAN  l4).  All  the  capabilities  of 
the  LOWTRAN  transmittance  codes  have  been  preserved  in  the  conversion  of  the 
computer  code  to  atmospheric  radiance  calculations.  The  LOWTRAN  4 program 
will  perform  either  atmospheric  transmittance  or  radiance  calculations.  The  mode 
of  the  calculation  is  determined  by  a single  input  control  parameter. 


(Received  for  publication  27  February  197  8) 
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PRECEDING  page  blank 


Continuity  with  the  LOWTRAN  3B  transmittance  code  has  been  retained. 

Changes  in  the  program  required  for  radiance  calculations  have  been  minimized. 

Only  two  new  subroutines  have  been  added  to  the  program.  Although  a new  card 

sequencing  has  been  established  in  the  LOWTRAN  4 program,  a comparison  with 

the  previous  LOWTRAN  3 B code  has  been  provided  (see  Appendix  A)  to  facilitate 

user  conversion  of  this  program  to  the  radiance  mode. 

In  this  report,  we  will  describe  the  changes  made  in  the  transmittance  code. 

These  include  a discussion  of  the  method  used  in  the  calculation  of  atmospheric 

radiance,  the  addition  of  UNO.,  as  an  atmospheric  absorber,  and  the  modification 

a 

of  the  empirical  transmittance  functions  for  small  absorber  amounts.  In  addition, 
instructions  for  using  LOWTRAN  4 are  given  in  Section  6.  A listing  of  the  com- 
puter code  and  data  is  given  in  Appendix  A (note  data  is  the  same  as  LOWTRAN  oB). 
A flow  chart  of  the  program  is  provided  in  Appendix  B,  and  a definition  of  symbols 
in  Appendix  C.  Sample  output  from  the  program  is  provided  in  Section  7.  Examples 
of  atmospheric  radiance  calculations  from  LOWTRAN  4 are  presented  in  Section  8. 
Comparison  with  measurements  are  given  in  Section  9.  Some  comments  on  the 
program  use  are  discussed  in  Section  10. 

The  latest  errata  sheet  (September  1 0 7 7 > to  the  LOWTRAN  3B  (1976)  code  has 
been  reprinted  in  Appendix  I).  I'sers  of  the  LOWTRAN  3B  program  should  verify 
that  the  corrections  listed  in  Appendix  D have  been  made  in  their  present  codes  and 
in  the  conversion  of  their  cedes  to  the  radiance  mode,  as  described  in  this  report. 

I'sers  of  the  current  LOWTRAN  313  transmittance  code  will  find  this  report 
sufficient  in  itself  for  performing  atmospheric  radiance/transmittance  calculations. 

For  those  unfamiliar  with  the  LOWTRAN  band  model  type  calculations,  more  com- 

1 -4 

plete  information  is  provided  in  the  previous  LOWTRAN  reports,  as  well  as  the 

5 

Optical  Properties  of  the  Atmosphere  Report. 

The  computer  code  has  been  tested  for  atmospheric  radiance  calculations  for 
various  atmospheric  paths.  However,  no  attempt  has  been  made  at  this  time  to 
optimize  either  the  computer  code  or  the  execution  time.  If  any  discrepancies  are 
encountered  or  problems  occur  in  the  use  of  the  code,  please  notify  F.X.  Kneizys, 
AFGL/OP1,  Hanscom  AFB,  MA  01731. 

The  LOWTRAN  4 card  deck  will  be  made  available  from  the  National  Climatic 
Center,  Federal  Building,  Asheville,  NC  28801  for  a charge  of  $20.  00.  (Please 
address  requests  to  Mr.  R.  Davis. ) 


5. 


M 


eClatchey,  R.  A. , l'onn,  R.W.,  Selby,  J.E.  A,  , Volz,  F.  E.  , and  Caring,  J.S. 
(197  21  Optical  Properties  of  the  Atmosphere  (Third  Edition)  AFCRL-72-04  97. 


8 


2.  ATMOSrHKRlf.  RADIAN  OK 


The  LOWTKAN  SB  transmittance  program  has  been  modified  to  calculate  at- 
mospheric and  earth  radiance.  A numerical  evaluation  of  the  integral  form  of  the 
equation  of  radiative  transfer  lias  been  added  to  the  program.  The  emission  from 
aerosols  and  the  treatment  of  aerosol  and  molecular  scattering  were  considered 
only  in  the  zeroth  order.  Additional  contributions  to  atmospheric  emission  from 
radiation  scattered  one  or  more  times  has  been  neglected.  Local  thermodynamic 
equilibrium  was  assumed  in  the  atmosphere. 

The  average  atmospheric  radiance  (over  a 20-em  ' interval)  at  the  frequency. 
F.  along  a given  line-of-sight  in  terms  of  the  LOWTKAN  transmittance  parameters 
is  given  by 

1 

hr)  = f dr  B(“.  T)t  + B(F.  T.  It*’  ill 

a . a s b t 


where  the  integral  represents  the  atmospheric  contribution,  and  the  second,  term  is 
the  contribution  of  the  boundary,  (for  example,  the  surface  of  the  earth  or  a cloud 
top)  and 


a 


T 


b 


average  transmittance  due  to  absorption, 

average  transmittance  due  to  scattering, 
average  total  transmittance, 

average  total  transmittanees  from  the  observer  to  boundary, 

average  Planck  (blackbody)  function  corresponding  to  the 
frequency  r and  the  temperature  T of  an  atmospheric  layer, 

temperature  of  the  boundary. 


The  emissivitv  of  the  boundary  is  assumed  to  be  unity. 

The  LOWTKAN  band  model  approach  used  here  assumes  that  since  the  black- 

body  function  is  a slowly  varying  function  of  frequency  we  can  represent  the  average 

value  of  the  radiance  in  terms  of  the  average  values  of  the  transmittance  and  the 

bluckbodv  function,  r . r . and  r.  varv  from  1 to  t*L  F*\  and  r along  the 

as  t as  t 

observers  line-of-sight.  For  lines  of  sight  which  do  not  intersect  the  earth  or  a 
cloud  layer,  the  second  term  in  Kq.  (l>  is  omitted. 

The  numerical  analogue  to  iCq.  (1)  has  been  incorporated  into  the  LOWTKAN  A 
computer  program.  The  numerical  integration  of  the  radiance  along  a line-of-sight 
for  a given  model  atmosphere  defined  at  N levels  is  given  bv 


0 


N-l 

1(7)  = Tj  (t  (i) 

i=l  ' 


/-  T(i)+T(i+ 1)  \ 

( Ts(i)+Ts(i+1)\ 

\1'.  2 / 

\ 2 / 

B(i/.Tb)Ttb 

(2) 


Thus,  the  spectral  radiance  from  a given  atmospheric  slant  path  (line-of-sight)  can 
be  calculated  by  dividing  the  atmosphere  into  a series  of  isothermal  layers  and 
summing  the  radiance  contributions  from  each  of  the  layers  along  the  line-of-sight, 
that  is,  numerically  evaluating  Eq.  (1).  This  can  be  clearly  seen  from  the  follow- 
ing simple  example. 

Neglecting  scattering,  consider  a three-layered  atmosphere  characterized  by 
temperatures  T,,  T,>,  and  T^  as  shown  in  Figure  1.  Let  Ty  T 2>  and  Tg  be  the 
transmittances  from  the  ground  to  the  boundaries  of  each  of  the  layers  respectively 
(see  Figure  la).  Figure  lb  shows  the  corresponding  case  for  an  observer  in  space 
(distinguished  by  primed  t values). 


OBSERVER  IN  SFftCE 


T T 1, 

* X 

1 T3 

* 1 

T 

l2 

, 1 

t2 

, ^ 

J .. 

7 ] 

;3  T, 

t 

GROUNO  OBSERVER  GROUND 


(a)  (b) 


Figure  1.  Upward  and  Downward  Atmospheric  Radiance 
Through  a Three-Layered  Atmosphere 


Then  from  Eq.  (2)  the  total  downward  spectral  radiance  for  an  observer  on  the 
ground  (looking  upwards)  is  given  by 

1(7)  i = (l-r1)B(7,T1)  + (Tj-fgiBfoTg)  + (r2-r3)B(7,  Tg)  . (3) 

Similarly  for  an  observer  looking  down  from  the  top  of  the  atmosphere  (see  Fig- 
ure lb),  the  total  upward  spectral  radiance  is  given  by 

1(7)  f = (1-T  J)B(7,T3)  + (r  J-f')B(7,T2)  + (F'-f')B(7,T1)  + r'  B (7,  Tb). 

(4) 

A comparison  of  Eqs.  (3)  and  (4)  shows  that  in  addition  to  the  boundary  con- 
tributions to  the  total  upward  spectral  radiance,  the  total  downward  and  the  total 
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upward  spectral  radiances  from  the  same  atmospheric  layers  are  not  the  same  but 
depend  on  the  position  of  the  observer  relative  to  a given  atmospheric  slant  path. 

In  the  LOWTRAN  4 radiance  program,  the  position  of  the  observer  is  always  defined 
by  the  input  parameter,  HI  (see  Section  6). 

3.  MODIFICATION  OF  TRANSMITTANCE  FUNCTION 

In  the  LOWTRAN  3/3B  transmittance  model,  * “ the  average  transmittance  r 
over  a 20-cm  * interval  (due  to  molecular  absorption)  is  represented  by  a single 
parameter  model  of  the  form 

F = (5) 

where  is  a wavelength  (or  wavenumber)  dependent  absorption  coefficient  and  u> 
is  an  equivalent  absorber  amount"  for  the  atmospheric  path,  which  is  defined  in 
terms  of  the  pressure  P(z),  temperature  T(z),  concentration  of  absorber  AL.  and 
an  empirical  constant  n as  follows 

n 

(6) 


If  Eq.  (6)  is  substituted  in  Eq.  (5)  and  n is  set  to  zero  and  unity,  respectively,  Eq. 

(5)  reverts  to  the  well  known  weak  line  and  strong  line  approximations  common  to 
most  band  models. 

The  form  of  the  function  f and  parameter  n was  determined  empirically  using 
both  laboratory  transmittance  data  and  available  molecular  line  constants.  In  both 
cases,  the  transmittance  was  degraded  in  resolution  to  20  cm'1  throughout  the 
entire  spectral  range  covered  here.  It  was  found  that  the  functions  f for  HgO  and 
the  combiner  contributions  of  the  uniformly  mixed  gases  were  essentially  identical, 
although  the  parameter  n differed  in  the  two  cases.  Mean  values  of  n were  deter- 
mined to  be  0.  9 for  HgO,  0.7  5 for  the  uniformly  mixed  gases,  and  0.4  for  ozone. 

For  sufficiently  small  values  of  the  argument  C^u/,  the  transmittance  in  LOWTRAN 
3B  and  earlier  models  was  set  to  unity. 

Since  the  LOWTRAN  4 program  now  calculates  radiance  as  well  as  transmittance, 
the  transmittance  functions,  f,  were  modified  for  radiance  calculations  from  at- 
mospheric layers  of  small  optical  thickness.  For  cases  where  (0.  999  * F * 1)  the 
transmittance  functions  now  have  the  analytic  form 
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w pli  a U,  OMM  and  l)  0.  Ml  fur  H.,1)  and  the  uniformly  mixed  puses  and  a 0.  Ohh  ami 
I'  I . 0!5  lor  it.'oiu',  Tins  pseudo -linear  approximation  lit  l\q.  (»)  IS  licit**  l ill  the  coin 
put  or  |hn»}'.imiu  lor  li'.iiisimll.'im't’::  Im'I  wrrn  0,  Odd  and  1. 

The  para  inclrrs  a a in  I b Wi-rr  dot  rnnuuMl  from  a least  .square:;  lit  of  tin*  elli- 
j-iriealls  derived  l runsiUllbi tire  funelion  ui  l\q.  (!>), 


1.  MTU  If. 

lUousu reiiu'Ul s tnaile  from  balloon  flijdils*',  have  shown  the  existence  ol  nitric 
acid  m the  earth's  atmosphere.  Althonp.h  tip  rie  arid  is  of  only  minor  importance 
m at  mo, spheric  transmittance  ca  leulat  ions,  it  has  been  shown  to  lie  a sip.ni  lira  lit 
source  of  si  r dosphorir  emission,  pa  rtienl  a r Iv  m the  at  luusphcric  window  rep,  ion 
from  10  to  P’.  Jim,  Therefore,  mlrie  neiil  has  been  added  to  the  l A >W ‘THAN  pro- 
pram  as  a.  separata'  at  mosphc  ri  e absorber. 

The  Iransmii.t  »nee  Pile  to  UNO.,  has  been  assumed  to  lie  in  tin*  weak  line  or 

' - 1 

linear  rep, ion.  Absorption  eoeltieienls  dipili/.ed  at  I>  em  intervals  lor  the  h.  5>/un, 
7,  hpni,  ami  ll.Mjim  bauds  of  I1NO,,  have  been  incorporated  iuU>  tlie  ]X)\VTI\AN 
prop, ram  as  a subrout  ine  (SU  MKOUT1 N K IINOM).  1'hese  eoi'l l ieienls  wore  obtained 
by  t lolilmnn,  Kyle,  and  Itononio'  by  fiUiup,  their  experimental  results  with  the 
statistical  haml  model  spproxi mat ioii. 

The  eoneent  rat  ion  of  atmospheric  nitrje  acid  varies  with  altitude  and  also  ap- 
p,-  ,-s  to  PepenP  on  latitude  ai\P  season.  I’ipnre  M shows  the  volume'  mixinp,  ratio 
profile  of  atmospheric  nitrie  aeiP  as  a funetion  of  aUituPe  from  the  luca yn I’ement s 
of  Kvans,  Kerr,  and  WarPle.'*  Pur  the  purpose  of  this  report,  we  have  ehosen  this 
profile  to  represent  a menu  nitrie  aeiP  profile  m the  LOWTHAN  prop, ram.  'This 
profile  appears  in  a data  statement.  in  the  main  propram.  It  a more  definitive  nitrie 
aeiP  profile  for  a piven  latitude  raid  season  is  available,  the  user  can  ehanpe  the 
nitrie  acid  concentration  by  : nply  replacuip  the  data  statement  piveu  in  the  main 
p v\  »p,  ra  m. 

The  inclusion  of  lulrie  acid  as  an  additional  absorber  will  modify  somewhat 
t ransmittaneefci  from  those  caleulateP  with  llu*  l.OYYTUAN  MM  (IdVlO  code  in  tile 
spectral  rep,  ions  desrribed  above.  Differences  in  transmit  I anee  values  will  only  be 
s ip,  ui  fie. nil  for  lonj»  atmospheric  limb  paths  passing  through  the  peak  of  the  nitric 
arid  profile. 

t».  Murerav,  h,  P.  , Kyle,  T.  lb  , Murrray,  l*  . II.  # ami  Williams,  W.IP.  (15M1M) 
Nitrie  aenl  and  nilrie  oxide  in  the  Iowan*  stratospliere.  Nature 

V.  liolilmiiu,  i\t  , Kyle,  T.P.,  and  Motiomo,  P.W.  ( I !’7  1 ) Statist  ical  band  model 
parameters  and  inirpruted  intensities  for  the  h.hji.  V.Ii-p,  and  11,-M/i 
hands  of  UNO  vapor,  Appl,  Opt. 

M.  Kv.uis,  \\  . P.  , Kerr,  ,!.M,  . and  Wardle,  P.  1.  tthVh)  The  APS  St  ratosphorie 

balloon  Measurements  Project:  Preliminary  ItesnltsT  Atmos|ilierK*  l'lnvTron- 
hiTmiI  "Service,  TJownsview,  Onbuuo,  I’nnaTla,  Tteport . TMo.  APliH  .10  \ 4. 
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5.  riuiutt  vm  modii  u aiion 

1 1 i\'h  i -Un  li>  (hr  I i'\\  1 U Ai\  S / .i  l ; t V ms in i| l ,imv  an. las  mvass.ii  \ iat  i i i 

aan\  arsiail  la  a ( anus  i tntt.uiaa  l .nhaitaa  aatln  Warn  m.ula  ill  snail  .i  w in  .it.  to  pan 
Sana  ilia  initial  l U>\\  in  tin-  aaijpunl  plan  i inis,  Appall. h\  A .vnt.uns  ,i  latina  al 
tin'  lA'Iil'UAN  aamput  aa  aa.la,  l OAA  I’ll  AN  -I  lapalhaa  with  tin-  lane  siihrant  mas, 
I't'IMT,  ANtll  It,  UMiVI,  ami  1 1 \ III.  A 11  lnnij.;h  ,i  nan  sat|Uallt'Uif;  has  haan  m.nla  in 
tin'  aa.la,  tin'  l.iyt  I'l'liiiuu  ui  Ilia  listtiij;  at  Appall. Il  s A lii.lmalas  Ilia  aainaspan.lain'a 
With  Ilia  1 1 '\\  rl(  \ N ;m  no.  la.  I 'halloas  m tha  an.  I a ,i  la  lli.lt  a.  it  at  I h\  Ilia  w a a,  I,  M )\\  , 
Hi  tha  hint  ts'limm. 

Twa  naw  slihi'i'iit  mas  h.iva  l>aan  a.hla.l  m tha  I I'WTII  AN  A aa.la.  Sithraul  i na 
i’.\  I'll  . lata  rm  un- n tha  aiunnnil.it  iv  a ..1  vs  a aha  a mnmnt  thian^h  a.iah  a I tha  la\  a is 
1 lit  a inai't  a>t  hv  tha  la'H'lia.t  atmasphai  a si. Hit  path.  Ilia  niiiainila  nra  stara.l  in  tha 
hint  ns,  \A  I' A I’ll,  tin-  a.iah  al  tha  nlsai'hlii|',  spaalas.  Siihi'.mt  tna  UNO.'  a.  a, ilia. I la 
lin.l  t ha  nit  rta  ha  t.l  nt  nw'i'jit  lan  'at  I laiaiit  s .is  ,i  tnnatn'uat  I i'ai|iiatia\  . 

In  tha  I'li'Sailt  pi'i'jt I'htu,  ,i  oiimpla  paalila  tal'  tin  UN''.,  talnina  linsilij;  i.itia  is 

J'l-avi.l.'.l  til  Ilia  .hit. I st  Itamalit  ut  tha  in.im  paat;rnm  U'A  I’ A IISMIX  ’ /!.  Tha 

a It  1 1 II. las  t'l  tha  IniMlip,  i it|."i  t'ai'i'.'spaii.l  ta  thasatarlha  ma.lal  .il  maspha  aa  alt  it  lutas , 

1 'ill'  twa  nan  a.mllal  p.i a.i mat aas  tll'AI!  Till  lUNIO  h.it  a haan  .i.Ulail  la  tha 
I'l't'Hl'.ilil.  Tlta-Si'  I'.ii'hinalai's,  main. la. I an  tha  Hast  ''antral  a, ml,  will  ha  lullt  as 
plalim.t  m Ilia  mst  fiiat  inns  p.ivan  in  'math'll  h. 


<i.  INSTRUCTIONS  FOR  liSINO  I.OWTKAN  4 


Tlu;  instructions  for  using  l.OWTRAN  -1,  with  the  exception  ol'  a change  in  a 
single  control  card,  are  essentially  the  same  as  those  lor  the  HOWTRAN  3,  3B 
transmittance  codes.  In  an  attempt  to  make  the  instructions  as  clear  as  possible 
and  to  provide  sufficient  information  for  utilization  of  the  computer  code  in  this 
report,  Section  5 of  the  LOWTRAN  d'  report  has  been  repeated  here  in  its  entirety. 
Changes  in  the  instructions  required  for  radiance  calculations  arc  indicated  by 
bars  in  the  margins. 

The  input  data  for  l.OWTKAN  4 arc  given  in  Appendix  A.  In  general,  it  is  only 
necessary  to  change  the  last  four  cards  (referred  to  here  as  No's.  1-4)  in  order  to 
run  the  program  for  a given  problem.  The  formats  for  tile  last  four  cards  and  their 
application  will  be  discussed. 


6.1  Input  Until  anil  Formats 

The  data  necessary  to  specify  a given  problem  are  given  on  the  hist  four  cards 
as  follows: 


CARD  1 

MODKl,, 

111AZK 

RO, 

T HOUND 

CARD  2 

1U, 

112, 

ANCil,K 

CARD  3 

VI, 

V 2, 

DV 

CARD  4 

1XY 

1TYPK,  LliN,  JP,  1M 

KAN Ci  1),  11KTA,  VIS 


Ml,  M2,  Ml,  Ml,,  IK  MISS, 
FORMAT  (1113,  2F10.  3) 

FORMAT  (OF  10.  3) 
FORMAT  (UK  10.  3) 
FORMAT  (13) 


Definitions  of  the  above  quantities  will  be  discussed  in  Section  l>,  2. 

If  the  quantity  MODKl,  given  in  CARD  1 is  set  equal  to  0 or  7 (which  is  the  ease 
if  meteorological  data  are  used  input  to  the  program),  then  the  above  card  sequence 
(and  format  for  CARD  2)  is  changed.  These  cases  will  be  described  in  Section  0.  3. 


6.2  Ikmii'  Instruction* 

The  various  quantities  to  be  specified  on  each  of  the  four  control  cards  (sum- 
marized in  Section  6.  1)  wild  be  discussed  in  this  section. 


i ' MiM  | 


MOM  1 , III  A . !• , I1A  I'l  . l.l'N,  IT,  I M.  Ml,  l\i:!.  Ma.  Ml,  Il'MISS. 

KO,  I IIOI'NI) 

1 hr  I > . 1 1 ■ ,1 1 ; n - 1 < ' I - MOM  I sele-ts  mu  i>l  llir  Sin  eco.;  raphienl  nuclei  ,il  inesplicees 
er  specifies  lii.il  mcteereleeical  .lata  arc  te  l'c  usc.l  in  place  el  tin'  slnn,lnr,l  iii.'.K  Is, 
1 1 1 A , I spec  ill  cm  u I ict  her  nc  reset  atteliil.it  len  is  te  he  iiielii.letl  m l lie  ea  I e ill  at  leu  er 
not.  I'.' r any  prehlem  the  at Iiiii,s|i|iei'ie  path  I'liist  he  specified  as  Cue  el  three  ivpcs 
aeeer. Ini;;  te  IlYl'l  an. I I IN.  II-'MISK  aiul  T HOI  :hi  I ' are  new  central  parameters 
ter  r.i.li  ai  ire  ea  I e nl.it  lens,  'I  lie  rest  n!‘  the  .plant  it  i es  riven  en  O Alt  1 1 I tiv  lit  ell  Oil 
he  letl  blank  u net  required)  previ.te  the  user  w it h apt  iens  te  suppi  ess  priutiue  v-M’). 
te  intermlA;  the  sis  stan.lar.l  me, let  atmespheres  (Ml,  M'.l,  Mtt>  an.l  te  input  a new 
me, Id  ati'iesphere  \1M,  Ml  ).  The  upturns  ler  the  aheve  parameters  all, I their  nr.es 
are  state, i an, I described  m detail  below  ; 

MOUI-'I  ~ 0 it'  meteetfeler  teal  .lata  are  specified  tier  heri/enl  al 
paths  nnU  ' ' . 

1 selects  TtiOl’KAl  MOMM.  A I MON  I 'll  tslt  h\ 

1’  selects  MIM  A I 1 I'MW  ShMMI'li. 

- It  selects  MIM  Allll.ll l‘  W1NTK15. 

* -I  selects  SlMtAllOTU'  NltMMKH. 

= a selects  SMtARi'Tk'  WIN  I'l'K  . 

- t>  selects  lt'tt'.l  i'.S.  ST  AN  DA  HI) 

• V it  a iu'u  me, lei  .tmesphere  (er  r.niiesetuie  data) 

is  te  he  inncrlcd, 

111  \ . . 1-'  - (i  means  tie  aereset  attenuation  inelmle.l  in  tin-  c.ileulat  iens. 

1 1 1 A / 1 1 ■'  I er  t’  it  arresel  attrnnat  ieu  i;t  required  (see  also  t'AUI'  -D, 

It  111  \ . ' I ' is  set  equal  te  1 er  l!  atul  visual  ramp'  (\  IS ) is  net  specified  an  t’AUI'  2, 
then  the  pre,;r.iiu  Mill  automat  icailv  Select  visual  ranees  el'  tilt  lun  er  h kin  respec- 
t ivel\  , 

ltlAAl1'  V Head  ethrr  arresel  nieih'l  into  the  propram, 

t 1' \ I’l-  I ter  a he ri ■■  enta I (eenstant  pressure)  path. 

2 ler  a vrrtiral  er  slant  path  between  tun  altitudes, 
a ter  a vertical  er  slant  path  te  spare. 

^Iii  ilicse  cartes  the  fermat  ler  lard  chances  isee 
Sect  ieu  tk  a, 


no  t is  t a in  ta  ril  I'em  111  iens) 


I 


Tin-  TYl'K  1 path  should  not  bo  confused  with  a lo up.  90°  path  whore  tin*  local  height 
of  the  end  of  the  I rujeetory  is  at  a significantly  different  height , In  snob  a on  so, 
specify  the  path  according  to  ITYl’K  = 2. 

LKN  = 0 for  normal  operation  of  program. 

LKN  - 1 selects  1 lie  downward  TYl’K  2 I, GNU  path. 

The  parameter  I ,KN  can  be  ignored  (tho.t  is,  left  blank)  for  the  majority  of  cases. 

It  need  only  be  used  for  a downward  looking  path  (H2<H1)  when  two  paths  are 
possible  for  the  same  input  parameters.  In  such  a ease,  a computer  printout  state- 
ment will  be  given  indicating  that  the  user  lias  two  choices  for  the  problem  and  that 
the  shorter  path  has  been  executed.  Set.  LKN  = 1 for  the  longer  case. 

,11’  = 0 for  normal  operation  of  program 

JP  - 1 to  suppress  printing  of  transmittance  table/or  radiance  table 

1M  - 1 when  radiosonde  data  are  to  be  read  in  initially 

IM  = 0 for  normal  operation  of  program  or  when  subsequent  calculations 

are  to  be  run  with  MGDKI,  - 7 

ML  - Number  of  levels  to  be  road  in  for  M011KL  - 7 

Note  that  IM  and  ML  are  only  used  when  MGDKI,  = 7 and  then 
only  on  the  first  calculations  when  the  data  are  road  in. 

Ml  - M2  = M3  = U for  normal  operation  of  program, 

The  parameters  IVll,  M2,  and  M3  can  each  take  integral  values  between  0 arid  6 and 
are  used  to  modify  or  supplement  the  altitude  profiles  of  temperature,  water  vapor, 
and  ozone  respectively,  for  any  given  atmospheric  model  specified  by  MODKL. 

For  example: 

Ml  = 1 selects  the  TROPICAL  temperature  altitude  profile 

Ml  - 2 selects  the  MlDLATITllDK  SUMMER  temperature  altitude  profile 

Ml  = 6 selects  the  1962  U.  S.  STANDARD  temperature  altitude  profile 

M3  = 1 selects  the  TROPICAL  water  vapor  altitude  profile 

M2  = 2 selects  the  MIDLATITUDK  SUMMER  water  vapor  altitude  profile 

M2  = 6 selects  the  1962  U,  S,  STANDARD  water  vapor  altitude  profile 

M3  = 1 selects  the  TROPICAL  ozone  altitude  profile 

M3  = 2 selects  the  MIDLATITUDK  SUM M HR  ozone  altitude  profile 

M3  - 6 selects  t lie  1962  U.  S,  STANDARD  ozone  altitude  profile. 

The  control  parameter,  IKMISS,  determines  the  mode  of  execution  of  the  program, 
IKMISS  = 0 for  program  execution  in  transmittance  mode 
IKMISS  = 1 for  program  execution  in  radiance  mode. 

A message  is  printed  to  the  user  on  the  output  file  indicating  the  mode  of  program 
| execution. 

HO  = radius  of  the  earth  (knO  at  the  particular  geographical  location 

at  which  the  calculation  is  to  be  performed. 
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1 

I 

1 

I 

I 

I 

1 


I 


1 1 I JO  is  left  lila  ilk,  (hr  pm  am  w ill  m-e  I In'  i.  it  I latitude  v.iluf  of  till I . I'D  kin  1 1 
MODl'I.  is  sot  equal  In  0 or  7,  Otherwise  tin'  earth  radius  far  the  appropriate 
standard  model  at  mospho  re  (sp.'eilied  hr  MOItl  | I will  he  u-  od. 

T1UH  !\l'  - teeipar  iMire  .if  the  earth  t 'h)  at  tilt'  toeation  at  ir  tilt  h tlra 
ealeulatiou  a-  In  he  pet  formed. 

I lUH  Trli  is  1 1 1 1 1 \ used  in  the  rath. Hire  mode  v > l the  proprnni  for  "hint  paths  whirl' 
iulrrsort  the  earth.  K TiU'l'NI'  is  left  blank,  tin'  program  will  us"  the  temperature 
of  the  first  at  mosphorie  (aver  as  t he  heuiitla  ry  temporal  lire. 

In  tin'  ease  where  MODl'I  * V,  the  new  at  niosphot'O  vmodol  t'r  rati  iosoudo  data' 
is  insert od  between  OAlil'K  1 ami  L’  tso"  Koetion  tk  li', 

h An  1 1 f in,  ir.  a n i ; i i ' , u a n i ;i  :,  n r.\,  \ is 


(A  AU  I > f!  is  tiseil  tit  define  the  seoluet  i'tt',:il  path  parameters  for  a siren  penhtem, 

111  initial  aliitiitle  tlou) 

II'.’  ! ina I alt  i I tale  tkin ) 

It  is  important  In  emphas i .’o'  here  tltat  in  tin'  r.nlianee  mode  of  prop,  rani  e\eenlit'n 
til 'MISS  n.  111,  the  initial  altitude,  always  defines  the  position  of  the  observer 
t . ' t'  seitst'i'f.  lit  ami  III’  eaiinot  be  used  i lit orohuneoabi y as  in  the  Iransmitt.ineo 
niotle. 


AN  i a r ~ Hitt  iai  /enith  auqle  ideon'es)  as  measured  from  111 
H AN  til'  ~ path  lenrt  It  tkin  1 

lil-'TA  earth  t-enter  anple  subtended  by  Ml  ami  III’  hlerroes' 

VIS  sea  level  visual  ranee  tkm). 

It  is  at'*  noeossary  to  speeifv  every  quantity  qiven  abor  t':  only  those  that  ;ub' 
quately  tleseribe  the  prohlriu  aeoordinp,  to  the  parameter  I l k Ml'  tas  desoi-jbed 
below); 


(11  Horizontal  Maths  tl  ft  I'll  -•  1 1 

hi)  spenty  111.  HAN'illainl  \ IS  only 

thl  If  nonstandard  inoteorolojpea!  data  are  to  ho  vised,  that  m,  if 
MvlDh'l  •••  0 on  VAlil'  I,  then  the  to  1 low iifr  parameters  must  be  spivifh'd  on 
t'AKh  111.  1’,  T.  I'M,  lx  11,  Mil.  Wf',  UK.  l;  A Nr  if  aooordin,:;  to  KOU  MAI' 

vM-'Hl.ii,  ’.J  !•■?*.  1.  iih'lO.H,  I'KItl.in.  Wlter.'  t I',  MM.  lilt.  Mil.  amt  UO  ai-o  file 

pressut'i'  (nib),  temperature  t"v').  view  point  I emperat  are  t'V\  relative  livisilidit y 
t "'■■).  II. ,1)  density  (sin  in  "')  an. I o'/mio  .lensitv  vqiu  m '')  rt'speet ively . 

Note  that  it  is  noeossary  tospeeify  all  of  the  quantities  uutlerl mod  rvitli  a full 
line  and  vine  of  the  quantities,  mnlet'liuod  with  a dashed  line.  If  the  ozone  density 
(WO)  is  not  known,  a value  oau  be  ehosen  from  one  of  the  standard  at mosplierie 
models  by  tisins  flu'  partuiirter  Ms  e i v'AKI'  1. 


! 


I 

I 


t 

i 


(2)  Slant  Paths  to  Spai  n (ITV  1 ' K 2) 

(a)  specify  111,  AN i U .K  tttiil  V IS 

(h)  specify  HI,  11MIN  and  VIS  (for  limb  viewing  problem  whet'*’  HM1N  is 
tlm  required  tangent  height  or  minimum  altitude  of  the  path  trajectory. 

(3)  Slant  Hatha  lU'twocn  Two  Altitudes  UTVHK  ill 
(a)  specify  111,  Ilf!,  ANtll.K  and  VIS 

lb)  specify  111,  ANHdd'!,  KANHK  and  VIS 
(i-l  specify  HI,  112.  HANOI'!  and  VIS 

Kor  rust's  (hi  and  (el,  tin*  program  will  ralrtilatr  112  and  ANiildi  respectively, 
assumniing  no  refract it'll:  thru  prorrrd  as  for  raso  (a).  This  niothoii  of  defining 
tin'  prohU’m  should  bo  list'd  wlii'n  rrfrai’t ion  olffots  art'  not  important;  for  example, 
for  ranges  of  a low  tons  of  km  at  7.0 nit li  angles  loss  than  80’'.  It  ran  also  bo  usod 
for  larger  angles  (including  00 J)  provided  that  tho  path  lit'S  within  ono  atmospheric 


lave  r. 

<t  1 > Spt'i’ify  HI,  H2,  HKTA  and  VIS.  l.oavt'  ANlU.K  and  KAN  til'' blank  in 
this  oast'.  This  method  i'an  bo  used  when  the  geometrical  configuration  of  tin* 
sourer  and  receiver  is  known  accurately.  but  the  initial  yrntth  angle  is  not  known 
precisely  due  to  atmosphorii*  refraction  efforts,  Kota  is  most  frequently  deter- 
mined bv  the  user  from  ground  range  information. 

In  the  oases  of  2(b)  and  2(d)  above,  the  subroutine  ANtildf  Is  cubed  in  the 
program  to  determine  the  appropriate  input  venith  angle  by  an  iterative  technique 
taking  into  aeeount  atmosphorir  refraetion. 

In  the  ease  where  IUOPKI.  - 7,  the  new  model  atmosphere  (or  radiosonde  data) 
is  inserted  between  HAKlfS  1 and  7, 


cAiinn  vi,  v 2 , nv 

'l'tu*  spectral  range  over  which  t eausmittauee  data  are  required  and  tho  spec- 
tral increments  at  which  the  data  are  to  be  printed  out  is  determined  h.v  OAU1)  3. 

V 1 = initial  frequency  in  wavenumbers  (ctu  ) 

V2  - final  frequency  in  wavenumbers  (cm  *)  where  V 2 ' VI 

1 IV  * frequency  increment  (or  step  si7.e)  (cm  ‘ ) 

(Note  that  i<  = 10  /X  where  r is  the  frequency  in  em  and  X is  the  wavelength  in 
microns,  and  that  1 IV  ran  only  take  values  which,  are  a multiple  of  (1.  ) 

HAH  1)4  IX  Y 

The  control  parameter  1XY  can  cause  the  program  to  recycle,  so  that  a series 
of  problems  rati  be  run  with  one  submission  of  l.OWTKAN.  Hive  values  of  1XY  01111 
be  used  to  provide  the  options  given  on  the  following  pages. 


Ill 
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IXY  = 0 or  blank  rani  to  niii  of  program 

= 1 to  srlort  a mow  YARD  3 an. I YARD  I only  (assuming  oilier 

para  I Motors  a or  unchanged) 

:!  to  soloot  a now  data  sequence  (YARDS  1,  Li,  a,  atiq  ■! ) 

3 to  soloi‘1  a now  DARI)  Li  an.)  t'AUI)  -I  only  (assuming  nt  1^-  i - 
par;  motors  arc  unchanged' 

= ■!  to  soloot  a now  l'AHl>  1 amt  d.\!!P  ■!  only  (assuming  otl;cr 
parameters  arc  unchanged' 

Thus,  if  lor  tlto  satpo  moito]  atmosphere  ant  typo  of  at mosphorR.  palp  ttto 
1 realtor  wishes  to  mako  further  1 ransmit  t anoo  caiculat  ions  in  diiTeroqt  spectra]  jn- 

f tervals  Y i 1 to  Yd1  etc.  and  lor  a different  stop  si'.’ o tip  ' etc.  1,  thoi\  1 XV  is  sot 

t 

1 equal  *a  1,  In  this  ease,  the  card  sequence  is  as  follows  and  can  lie  pi'peatcd  as 

| many  times  as  required, 

! 

1 YARD  -1  IX Y - 1 

t'AUI)  f>  Yl'  YL’1  DY 

CAH|)  (i  IX Y = 1 

YAH  |)  7 Yl"  v;>"  nv” 

YAH  DH  IX Y * 0 

i 

l Tito  1 1 rial  IXY  card  should  always  hi-  a blank  or  aav,  When  using  the  h\Y  ~ 1 option, 

ttie  wavelength  dopotulonee  of  the  refraotive  index  is  not  changed  (uSe  IXY  1!  option 
if  this  is  required!. 

To  make  successive  transmittance  computations  whore  just  the  geographical 
model  atmosphere  is  changed  and/or  with  or  without  aerosol  atl emtatirin,  set 
IXY  ” -1  and  eonstruot  a data  oard  sequonee  along  the  same  lines  as  giv-g  above. 

This  sequence  of  recycling  can  ho  rep.eateii  successively. 

When  ;i  series  of  problems  is  to  he  oxocntod  (with  one  submission  of  1 .OWTUAN) 
involving  the  standard  atmospherie  models  (MODKI.  - 1 to  t>)  as  well  as  rases  in- 
volving MODKI,  - 0 and  MODKI,  = ",  then  the  order  in  which  the  data  are  sol  up 
becomes  very  important.  Note  the  following  sequence. 

1.  linn  all  problems  using  MODKI,  1 through  !!  first. 

3.  Secondly,  run  all  problems  involving  the  use  of  MODKI.  = 0, 

3 . Run  all  problems  involving  the  use  of  MODKI,  ’ 7 last.  The  reason  for 
running  MODKI,  = 7 cases  last  is  that  when  a new  atmospherie  mode]  is  road  in, 
the  altitudes  may  not  correspond  with  those  given  in  the  standard  nujilola  and  the 
program  will  erase  thorn.  Similarly,  if  a MODKI,  “ 0 ease  is  run  following  a 
MODKI j = 7 ease,  the  first  level  of  MODKI.  7 is  erased, 

<i.;t  Mm-Shimlaul  ('.omlilionn 

Three  options  are  available  if  atmospherie  transmittance  calctil at ions  arc  re- 
quired for  non-standard  conditions.  Here  non-standard  refers  to  conditions  other 
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th.ui  those  li od  by  Hu'  six  model  atmospheres  provided  hv  I OWTHAAJ,  which 

art'  selected  bv  the  parameter  MOMKl.  on  v’AKI'  1.  The  llnvc  options  enable  the 
reader  U'  insert : 

it'  Hi  sown  model  utniosphoi'C|sl  ii\  place  of  any  for  all'  of  flu'  six  standard 
models,  provided  that  tho  data  at'o  in  exactly  tho  satin'  format  and  an,'  specified  at 
tin-  s tun'  altitudes  as  tho  latlor.  In  this  raso  tin-  appropriate  print  statements  in 
1 f'WTKAN  tthat  nlmtifv  tho  atmospheric  niodol  usedl  must  ho  changed  correspond- 
ingly , 

iLM  An  additional  at ,i nosphoi’ie  modt'l  (1UOHK1,  71,  whioh  ran  bo  in  tho  form  of 
vndiosopdo  data,  Tho  data  nood  not  bo  sporifitnl  at  flu'  satno  altitudes  as  in  tho 
standard  mo. tots, 

td'  Meteorological  conditions  for  a given  horizontal  path  calculation  (MOHl’h  ■ 

0 ease'. 

Tito  Host  of  those  out  tons  requires  tho.  most  effort  and  needs  no  further  dismis- 
sion hero,  other  than  a reference  to  Appendix  A for  a summary  of  the  standard 
model  atmosphere  pat'.tmotoos,  units,  and  formats. 

AliPlTlONAl.  ATMOSl'llh'UU'  MOlil'l,  tMOOKl.  - 71 

Now  niodol  atmospheres  rati  bo  inserted  between  o'. A HI  IS  1 and  LI  provided  the 
parameters  MOIH-1  and  IM  arc  set  equal  to  7 and  1 respectively  onv'AliP  1.  The 
number  of  atmospheric  levels  to  bo  inserted  tMl.l  must  also  be  specified  on  t’AHD  1. 
The  appropriate  met oorologioal  parameters  and  format  for  tho  atmospheric  data 
are  given  below. 

11’  !’•  !•  !>!*.  till.  W U,  WO,  AHA/h  I KOU  MAT  tdl'Ul.d,  LtKfi.  1,  L'blO.d, 

LL  t*'  10 . a)  1 \v here 

- nltit udo  t km ' 

1’  pressure  fmh' 

T ambient  temperature  1'f'l 

l'l’  dew  point  temperature  ('V> 

lilt  relative  humidity  r’  l 

Wll  water  vapor  density  lean  m ' 1 

WO  ozone  density  tgm  m ' I 

AM  A OK  aerosol  number  density  tom  ' ) 

Note  that  it  is  only  necessary  to  specify  those  quantities  underlined  with  a full  lino 
and  either  of  tho  quantities  underlined  with  the  dashed  lino. 

II  the  ozone  density  tWOl  is  not  known  then  a value  can  be  obtained  from  one  of 
the  standard  atmospheric  models  (for  the  appropriate  latitude  and  season!  by  using 
the  parameter  MA  on  I’AUl)  1, 

Lid 


It'  t iii'  ,u  iv  i'l  > mi ul v i'  .li’iis it y i\  as  not  moasnrod  as  a fun .'t  ion  of  a ll  it  ndo  an, I 
tin-  i,  a, ior  ’..ishos  to  mohi.lo  aoiosol  aUonuaUon  In  Ilia  ,':tla\it  ation,  sot  1IIA.  P 1 

on  t’AR  l>  l.  In  this  ,■  a s * ' tas  w it  It  tho  Ml,  \'.t.  and  Mil  I'ptioi'.s)  IOW  ] UAN  will  nso 

ttio  aoconol  mo, lots  alroady  oont.ii in-, I in  tho  proor.nii  ami  uitoopol at o to  pivo  aoiosol 
inimhor  donstlv  vnltios  at  tho  sat-.io  altitudos  as  iho  rydh'Snndo  too  now  mo, lot  at - 
mospliorot  itaia,  Tho  program  will  thon  loot.  for  a sou  lovrl  visual  rai1.^'  f\  tSI  to 

t'o  spooifiod  on  o'  A it  1 1 '.t.  If  Vis  is  not  spooi  find,  a.  I'l'  kin  son  tovol  visual  nauj’.o  w dl 

ho  assumo.l.  If  aoi'osol  attonuation  is  iivii  roipti ivd.  Sot  lllA.'P  0 on  P.\KI>  t a* 
l't’foi’o. 

HOUIPONTAL  PA  ft  IS  |MOl>KI,  - 01 

il  tni'l  oorolopionl  data  a tv  to  l'o  list'd  for  liori/ontat  path  atiuosphorio  trans- 
tnitl  .moo  onhuilat  ion*,  thon  sol  MOl'll  = o'  on  PAlll'  l,  'tho  loltowino,  pa  ranu'tors 
run  tln'n  ho  spooifiod  on  PAivU  J; 

< Ai; 1 1 m.  v,  v.  pp.  hjj,  \vu,  wo,  us.  itAMPf  irt'nm-vr  tit,  kip,  n. 

it  Pa.  1,  '.M'lO.o,  'Jl'lO.I'lj  win' ri’  Iho  nt'Ovo  paranmtors  rol'or  to  alutu.lo  ikm', 
prossuro  pnM  ainliii’ul  loinporatin  o r'p),  dow  point  lon'd'rraturr  ('pi,  rolalivo 
humidity  v''.  >.  w ator  vapor  donsity  (.;m  tn  w.-iu,,'  .lonsity  fom  m'0),  visual  rano.o 
iknt)  an.i  path  lonrth  (Uni)  rospootivoly . 

’I  Itr  format  for  iho  al-ovo  oard  is  similar  to  that  for  inputting  raiiiosomto  data 
iMPUtl'l  " it*  A pa  in,  it  is  only  iii'ocssiirv  t.'  spoo'fy  tlio  ipiantitios  umt'.'rl  inod  witli 
tho  solid  tiro  and  ono  ol  tho  ipiantitios  imdorlinod  w it  It  tho  dastiod  tino,  Tlio  o.'ono 
dons  it  \ WO  oan  t'o  spooifiod  using  tho  paramotor  j\Ul  on  v ’ At?  1 > 1 if  mrasuromoiits 
aro  not  availat'lo.  In  tlio  lattor  oaso,  a valno  will  t'o  ralonlatrd  at  ultitndo  lit  hast'd 
on  tho  appropriato  modol  atmosphoro  soh'otod  hv  Mil, 


7.  l\  \Ml*l  If  Of  nUHdl  \M  ISI. 


7.1  l*t  obli'iii 


Paioutato  tim  i raiis n i it i ain’t*  ft' own  .toT'O  to  'Jthp  om  ' in  stops  of  a oro,  ^ for  a 
slant  path  from  It,  S km  to  11.  S tun  at  a '/onitti  anglr  of  t>fo',  for  a suharrtir  vvintor 
modol  dinosphoro.  and  a 'Jit  Utn  visual  rangr.  Hr  pout  tlio  oaloutat  it'll  for  tlio  saint' 
oonditi.oiis  ovoout i tig  tho  program  in  Hit'  radianor  iiiotio* 
t'AtiP  1 r>  !+♦:! 


o' ap,  n a 


P \K  I ' It 
v’  Alt  1 1 -1 
o'  AH  I ' h 
PA  I!  I ' il 


+ » + . r,p,,  .tvp  p.a  ip  r, 0 V SO  V 0 V ti f.  a 

* ***7'A50,  ot'VM'2-ir'O.  O’*  + < t"t  'Sif'.  0 

v'  t 

v s to'  o t i 0 p-:,  + + v o r v o v p -v a -,\  \\  a a v s p t ' . , hi : n n I' o' 
Pd  A hits 

span*  on  tin-  oard' 


:*  t 


7 .2  Out  pul  from  I.OWTKAN  4 

Tin-  output  lor  this  problem  is  given  in  Table  1.  A message  indicating  the 
mode  of  execution  of  the  propram  is  printed  as  the  first  line  of  the  output,  h'or  this 
problem,  the  first  ease  will  be  executed  in  the  transmittance  inode. 

The  parameters  defining  the  atmospherie  path,  model  atmospheres  and  fre- 
quency range  are  next  printed  out,  Pol'inving  the  heading  HOltl/.ONTAl,  PHOKJl.KS 
there  are  IS  eoltnuns,  The  first  eolmnn  gives  a running  integer  ussoeiated  with 
eaeh  level  (level  indieator).  The  seeond  eolumn  gives  the  level  altitude  in  km.  The 
next  it  eolumns  give  the  equivalent  absorber  amounts  per  km  Tor  the  following  ab- 
sorbing species;  water  vapor,  uniformly  mixed  gases,  ozone,  nitropen,  eoutimuun, 
water  vapor  eoidinuum  (10  pm),  moleeular  scattering,  aerosol  extinction^  and  UY 
ozone,  respeetively . The  next  three  eolumns  pivo  th.e  mean  refraetive  index  modulus 
from  that  level  to  the  level  above,  the  equivalent  absorber  amounts  per  km  for  the 
water  vapor  eontinuum  (4  pm)  and  for  nitric  arid, 

A heading  V Ell TICAL  PltOKIEKS  is  then  printed  followed  by  If'  columns.  The 
first  and  seeond  columns  give  the  integer  associated  with  the  levels  traversed  by 
the  path  and  the  height  of  the  level.  Then  follow  ft  columns  which  give  the  integrated 
equivalent  absorber  amounts  from  the  initial  altitude  to  the  level  above  (in  the  same 
order  as  indicated  above).  The  next  4 columns  are  labelled  PS1,  Till,  BETA,  and 
TH11TA,  and  correspond  to  the  angles  similar  to  i p,  0,  jl>  and  d described  in 
l.OWTUAN  2,  “ Columns  PS1  and  BETA  give  the  accumulated  values  of  ip  and  ji  io 
the  level  above.  Columns  THETA  and  PHI  give  the  local  zenith  angle  dj  correspond- 
ing to  that  level  and  the  angle  of  arrival  at  the  level  above,  respectively.  The 
accumulated  slant  range  is  printed  out  in  the  last  eolumn  under  KANC1K. 

The  total  equivalent  absorber  amounts  for  eaeh  absorber  species  are  then  sum- 
marized below  in  their  appropriate  units. 

The  seeond  line  in  the  total  equivalent  absorber  amount  table  gives  the  water 
vapor  eontinuum  amount  (4  pm)  and  the  nitrie  acid  amount. 

A transmittance  table,  containing  12  columns,  now  follows.  The  first  2 eol- 
umns give  the  frequency  (cm  wavelength  (pm),  and  total  transmittance.  '1  he 
next  7 columns  show  the  individual  transmittance  due  to  water  vapor,  uniformly 
mixed  gases,  ozone,  nitrogen  (4  pm)  eontinuum,  total  water  vapor  eontinuum,  mole- 
cular scattering,  and  aerosol  extinction.  The  last  2 columns  give  absorption  due 
to  aerosols  and  the  cumulative  integrated  absorption.  The  latter  quantity  can  be  used 
to  determine  the  average  transmittance  over  any  given  speetral  interval  within  the 
spectral  range  covered  by  the  calculation.  Kitially,  the  total  integrated  absorption 
from  VI  to  V2  is  printed  out  (units  are  cm  ) together  with  the  average  transmittance 
over  the  band. 

^ f or  all  radiance  calculations  in  this  report,  the  average  continental  aerosol  model 
was  used. 


Table  1.  Typical  Output  of  I.OtvT RAN  4 
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The  second  case  is  now  executed  for  the  sumo  conditions  in  the  radiance  mode. 
Idle  output  of  the  program  is  identical  to  the  transmittance  mode  up  to  and  including 
the  printing  of  the  equivalent  sea  level  absorber  amounts. 

Two  parameters,  .11  and  .13,  are  then  printed  out.  These  parameters  control 
the  loading  of  the  cumulative  absorher  amounts  into  the  matrix,  WFAT11, 

A heading  OUMIH  .ATIVF  ABSORBKU  AMOUNTS  KOI?  THK  ATMOS  I’ll  KKIC 
PATH  is  then  printed  followed  by  13  columns.  The  first  column  gives  an  integer 
associated  with  the  layer  traversal  by  the  atmospheric  slant  path.  The  following 
10  columns  which  give  the  cumulative  absorber  amounts  for  the  following  species: 
water  vapor,  uniformly  mixed  gases,  ox  one,  nitrogen  continuum,  water  vapor 
continuum  (10  pm),  moleeplar  scattering,  aerosol  extinction,  UV  ozone,  water 
vapor  continuum  (4  /mi)  and  nitric  acid.  The  last  column  Is  the  average  tempera- 
ture of  the  layer. 

A radiance  table,  containing  six  columns,  now  follows.  The  first  two  columns 
give  the  frequency  (cm  *)  and  the  wavelength  (pm).  The  next  two  columns  give  the 

2 „ j 2 

radiance  in  units  of  W / cm" •■ster -cm  and  W/ cm  -ster-pm.  The  next  column  gives 

•> 

the  cumulative  integrated  radiance  (VV/oin“-ster).  The  last  column  is  the  total 
transmittance. 

Finally  the  maximum  and  minimum  radiances  and  their  frequencies,  the  inte- 
grated absorption,  the  average  transmittance,  and  the  total  integrated  radiance  are 
printed. 

8.  F.XAMI’LKS  Ob  RADIANC.F,  Sl’F.CTRA 

Some  examples  or  radiance  spectra  obtained  from  BOWTRAN  4 are  presented 
in  Figures  3 through  0,  Figures  3 and  4 show  the  atmospheric  radiance  as  seen 
by  an  observer  at  the  ground  looking  straight  up  to  space  (1 1 1-0  km,  zenith  ANtil  ,14-0'’ ) 
for  the  six  model  atmospheres  with  a 33  km  visual  range.  Figure  3 is  for  the 
spectral  region  from  400  to  3000  cm  * a Figure  4 for  the  spectral  region  from 
2000  to  3000  cm  *. 

Figures  0 and  6 show  the  atmospheric  radiance  as  seen  by  an  observer  in  space 
looking  straight  down  to  the  ground  (111s  100  km,  zenith  ANGI  = 1110'’)  for  the  six 

model  atmospheres  with  a 23  km  visual  range.  The  temperature  of  the  ground  for 
these  plots  is  the  appropriate  boundary  temperature  of  the  first  layer  in  die  model 
atmosphere.  Figure  f)  is  for  the  spectral  region  from  400  to  2000  um  * and  Fig- 
ure 8 for  the  spectral  region  from  2000  to  3000  cm’*. 

Figures  7 through  !)  show  atmospheric  radiance  spectra  from  400  to  4000  cm 
calculated  using  the  U.S,  Standard  Atmosphere,  1002,  with  a 23  km  visual  range 
for  three  different  types  of  atmospheric  paths.  Figure  7 shows  the  zenith  radiance 
for  an  observer  at  altitudes  of  0,  20,  and  40  km.  Figure  11  shows  the  atmospheric 


ra.liauoo  ai  a /onifh  ani;le  of  4!i'  ;is  soon  by  an  observer  at  altiUulos  of  0,  20,  ami 
■10  km,  1'ii.nn'o  0 shows  tin1  ooinparisons  of  flu-  limb  railianro  as  soon  from  spare 
lor  tangent  heights  of  0,  20,  ami  -10  km. 


HftVt  NUMBER  (CM  1)  lMO'1 


Kirure  2.  Atmosphorie  HaOianoo  for  a Yortioal  Path  to  Spare 
from  (.ironnil  l.ovol  for  Six  lYlotlel  Atmosphoros  (400  to  2000  om-J) 


• i ' 
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Figure  6.  Atmospheric  Radiance  for  a Vertical  Path  to  Ciround 
from  Space  for  Six  Model  Atmospheres  (2000  to  0600  cm'*) 
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Figure  7.  Variation  of  Atmospheric  Radiance  With  Altitude  for 
Vertical  Paths  to  Space  and  the  1062  H.S.  Standard  Atmosphere 
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Figure  8,  Variation  of  Atmospheric  Radiance  With  Altitude 
for  Slant  Paths  io  Space  (Zenith  Angle  = 4S°)  and  the  1862 
U,S.  Standard  Atmosphere 


Figure  0.  Variation  of  Atmospheric  Limb  Radiance  With  Tangent 
Height  for. the  if>62  U.S.  Standai  d Atmosphere 
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Kigures  It)  through  show  some  comparisons  of  I.OWTHAN  -I  calculations 
with  measurement s of  atmospheric  spectral  radiance  from  both  balloon  anil  Satellite 
plat  forms. 

h'igure  10  shows  a comparison  of  the  calculated  upward  atmospheric  radiance 

with  an  interferometer  measurement  from  a balloon  flight  over  northern  Nebraska 

o 

by  Chaney  at  the  University  of  Michigan, ' The  measurement  was  taken  at  a float 
altitmle  of  1 1 1.700  ft.  The  calculated  radiance  used  the  midlatitude  winter  model, 
with  a dd  km  visual  range,  and  a ground  temperature  of  dBO'K. 

Figure  11  shows  a comparison  of  an  interferometer  t.  asnrement  made  from 
the  Nimbus  3 satellite  looking  down  over  the  ilult'  of  Mexico  with  the  ealeulatod 
ntmospherii'  radiance.  The  resolution  of  the  interferometer  was  P cm  as  com- 
pared to  the  dd  cm  ' resolution  of  l.OvVTR  AN.  Two  theoretical  models,  the  tropi- 
cal and  midlatitude  summer,  wore  used  for  comparison,  as  shown  in  Figure  11, 
and  are  displaced  two  divisions  above  and  below  the  measured  radiance  for  clarity, 
both  models  assumed  a dd  km  visual  ratine  and  used  the  temperature  at  OKM  in  the 
model  atmosphere  as  the  boundary  temperature. 

Figure  I'd  shows  the  comparison  of  atmospheric  radiance  as  seen  from  space 
between  the  I.OYVTHAN  calculation  and  measurements  from  the  Nimbus  -1  satellite^ 
for  three  different  geographic  locations,  t he  spectra,  obtained  with  a Michelson 
interferometer  of  resolution  d.  tl  cm  *,  were  measured  over  the  Sahara  Desert,  the 
Mediterranean,  and  the  Antarctic,  The  calculated  l.OWTRAN  radiances  used  the 
midlatitude  winter  model  and  a ground  temperature  of  3 d 0 '' K for  tlie  Sahara;  the 
midlatitude  winter  model  and  a ground  temperature  of  dBT'K  for  the  Mediterranean; 

and  an  arctic  winter  cold  model  taken  from  the  AFCRL  Handbook  of  Geophysics  and 
1 11 

Space  1'uvironments  and  a around  temperature  of  100'K  for  the  Antarctic  com- 
parison, All  three  calculations  assumed  a d3  km  visual  ranee  for  aerosols, 

Figures  id  through  dd  show  a detailed  comparison  of  calculated  and  observed 
atmospheric  spectral  radiance  vs  wavelength  in  both  the  ii  to  H pm  and  the  IB  to 
'.!7  pm  spectral  regions.  The  measurements  were  made  on  a balloon  flight  launched 

0.  Chaney,  1..W.  (1000)  An  experimental  Fourier  Transform  Asymmetrical  Inter- 
ferometer  for  Atmospheric  RadiafToTTriTensuromoids,  University  of  Michigan 
Technical  Heport  htditid - llt-'r, 

10.  ,’onrath.  11..1.,  llanel,  H.A.,  Kundc,  Y.C.,  and  I’rabhakara,  C.  (10701  The 

Infrared  Interferometer  experiment  on  Nimbus  3.  Goddard  Space  blight  Center 
T r e cube  It,  Maryland,  Report  N -(120  -70 -li  13. 

11.  Hanoi.  H.  A.  , and  Gonrath,  11. .1.  ( 107  01  Thermal  Emission  Spectra  of  the  Forth 

and  Atmosphere  Ohtaino , 1 from  the  Nimbus  4 kTichclson  Interferometer 
experiment,  Goddard  S | >are  Flight  Center,  (“ireenbolt,'  Maryfand,  lfepo r t 
TT-R20 -70-244. 

Id,  Valiev,  S,  D,  , ltd.  (1000)  Handbook  of  Geophysics  and  Spare  Knviromnents, 
Ah’CUl ,. 


from  Holloman  A KB,  Now  Mexico  by  Mureray  et  al,  University  of  Denver,  The 

instrument  used  for  these  observations  was  a Idle  grating  spectrometer,  operated 

in  tlie  first  and  se<  end  order  of  the  grating.  The  resolution  was  0.03  /mi  in  the 

(t  to  1-1  i im  region,  and  0.  0G  /tin  m the  ill  to  27  /un  region.  Tlie  ilata  in  these  figures 

are  pi'eseuted  as  a function  of  altitude  and  as  a function  ot  zenith  angle.  Figures  13 

through  ltl  rover  the  11  to  14  /Jin  region  and  figures  ltl  through  22,  tlie  Hi  to  27  /mi 

reeion  The  DOWTUAN  radianee  calculation  used  the  pressure,  temperature, 

*'  i;t 

ozone,  and  nitric  arid  profiles  from  tlie  Mureray  report,  and  the  niidlatitude 

winter  water  vapor  profile  eontained  in  l,tl\VTHAN, 

Tlie  dominant  spectral  features  in  Figures  13  to  18  are  the  !>.  0 and  8 /tin 

o/.one  hands,  the  11.3  pm  hand  of  nitric  arid,  and  the  earhon  dioxide  hands  between 

12  and  14  p. m.  hi  the  long  wavelength  region,  the  spectral  radiance  is  due  primarily 

to  water  vapor  rotational  transitions.  Deeause  the  resolution  of  the  DOW  TUAN 

radianee  ealeulation  is  inferior  to  the  measurements  in  this  region,  the  comparisons 

here  can  only  serve  to  verify  the  level  ot  tlie  calculated  radianee. 


' ; f '-'v 

\\  t h'V, 

V \!i !: 

i?  % 


t i-iiW, 


I . I i 1 I . I 1 I ; I 1 1 I 1 1 1 , 

500  600  700  800  900  1000  1100  1200  1 500  1400  1500  1600  1700 

WAVENUMBER -cm-' 

Figure  10.  Comparison  Hetween  DOWTUAN  Prediction  ami  University  of  Michigan 
Balloon  Measurement  of  Atmospheric  Radiance  Over  Northern  Nebraska 


13.  Mureray.  D,  G,  , Brooks.  N.  , Goldman,  A,,  Fosters,  .D,D,  and 

Williams,  W.,1.  <10771  Water  Vapor  Nitric  Arid  niiii  O?.oilo  Mixing  Ratio 
Height  Profiles  Derived  from-  Spectral  Radiometric  Measurements. 
Tlnlvor'sity  of  'Heaver,  Denver,  Colorado  itoiioil.  Contract  Report  No.  332. 
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Kitfuro  12.  Comparison  Between  l.OWTHAN  Predict  ion  ami 
NIMBUS  4 Satellite  Measurements  of  Atmospheric  Hiulianee 
over  the  Sahara  Desert,  the  Mediterranean,  and  the 
Antarctic 
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Figure  1H.  Sample  Spectrum  of  Short 
Wavelength  Region  Observed  at  an 
Altitude  of  0,  5 kin  and  a Zenith  Angle 
of  6 3 ■'  on  19  February  197  5,  and 
LOWTRAN  Comparison 
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Figure  14.  Sample  Speetruin  of  Short 
Wavelength  Region  Observed  at  ail 
Altitude  of  13.5  km  and  a Zenith  Angle 
of  610  on  19  February  197  5,  and 
LOWTRAN  Comparison 
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Figure  15.  Sample  Spectrum  of  Short 
Wavelength  Region  Observed  at  an 
Altitude  of  19.  0 km  and  a Zenith  Angle 
of  63"  on  19  February  1975,  and 
1 iOWTRAN  Comparison 
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Figure  16.  Sample  Spectrum  of  Short 
Wavelength  Region  Observed  at  an 
Altitude  of  24,0  km  and  a Zenith  Angle 
of  63"  on  19  February  197  3,  and 
LOWTRAN  Comparison 
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Figure  17.  Sample  Spectrum  of  Short 
Wavelength  Region  Observed  at  an 
Altitude  of  251,  i kni  and  a Zenith  Angle 
of  93°  on  19  February  1975,  and 
LOWTRAN  Comparison 
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Figure  19.  Sample  Spectrum  of  l.ong 
Wavelength  Region  Observed  at  an 
Altitude  of  9,  5 km  and  a Zenith  Angle 
of  63°  on  19  February  1975,  and 
l.OWl'RAN  Comparison 
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Figure  111,  Sample  Spectrum  of  Short 
Wavelength  Region  Observed  at  an 
Altitude  of  29.  1 lun  ami  a Zenith  Angle 
of  111"’  on  li>  February  1975,  and 
LOWTRAN  Comparison 


— — MIIRCKAT  II  .il.  HO  II 0 MAN  Alp.  N!IV  MfXICO 
|i  I I RROARN  I'm 
1.0  ft  1 R A N 


W.m  LirNC.T  H l,,ml 

h’iguro  20.  Sample  Speetruni  of  Long 
Wavelength  Region  Observed  at  an 
Altitude  of  20.0  km  and  a Zenith  Angle 
of  63"  on  19  February  197  5,  and 
LOW  111  A N Comparison 


— M;iKi*K&N  \ 1 Al.  MO l i.0 MAN  A<  B.  NFW  MEXICO, 
IQ  M BRl'AKX  I *>/*> 

- 10\\  I KAN 


Figure  22.  Sample  Spectrum  of  Long 
Wavelength  Region  Observed  at  an 
Altitude  of  US',  1 km  and  a Zenith  Angle 
of  ill'"  on  10  February  107  5,  and 
1 ,OWT  R A N C omparis  on 
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Figure  21.  Sample  Spectrum  of  Long 
Wavelength  Region  Observed  at  an 
Altitude  of  29.  1 km  and  a Zenith  Angle 
of  93°  on  19  February  1975,  and 
LOWTRAN  Comparison 


III.  HHIMIM'S 


TIh‘  I OWTIi  AX  model  \\ originally  I’anrciv cd  as  an  atmospheric  transmittance 
mo. Id.  In  tlio  conversion  to  a radiance  model,  additional  constraints  on  both  the 
\ all. lit  v of  tho  model  as  well  as  the  ranpc  ol'  applieabi  I il.v  arc  i tit  reduced  lor  alim>s- 
pherte  radiance  oaleulat  ions.  It  shotthi  hr  remembered  that  I hr  dipiti/ed  molecular 
absorption  coefficients  in  l.OWTH.AN  were  obtained  for  conditions  represent  stive  of 
meliorate  atmospheric  paths  anil  wall  leiiil  to  mule  rest  i mat  e the  transm  ittanee  for 
vers  Iona  paths  and  overestimate  the  1 ransmittanee  for  verv  short  paths.  The 
lnoilifiration  of  the  (ransmittanee  function  ileseribosl  m Section  H was  inaile  to'pivo 
some  improvements  to  the  raitianee  calculations  for  the  short  path  eases.  In  addi- 
tion, the  radiance  calculations  assume  local  thermodynamic  e.jUiHbriuni  exists  in 
each  layer  of  the  model  atmospheres.  This  assumpt  ion  will  break  down  for  tadianee 
calculations  in  the  upper  atmosphere , Therefore,  beeause  of  the  limitations  in  the 
l.OWl'U'W  model  for  short  paths  (or  small  absorber  amounts!  and  deviations  from 
thermal  equilibrium,  belli  renditions  whieh  oeeur  in  the  upper  atmosphere  it  is 
recommended  'hat  the  lOWTHAN  radiance  calculations  he  restricted  to  altitudes 
below  -10  km. 

b'or  the  shorter  w av oie n:;t hs  (<  5 pm),  scattered  solar  radiation  becomes  an 
important  soinvt  of  barkprouud  radiation.  Since  this  lias  not  been  included  in  the 
l.OWTRAN  model,  radiance  calculations  at  the  shorter  wavelengths  with  ;i  sunlit 
atmosphere  should  be  made  with  caul  ion. 

With  the  obvious  limitations  in  the  1 .OWTR A N radiance  code  described  above, 
t'ne  agreement  between  measurements  and  radiance  calculations  shown  in  l-’it'tires 
10  throuph  'J2  is  remarkably  pood  for  the  eases  considered,  l-'urther  comparisons 
with  measurements  are  planned  for  other  spectral  regions  and  other  peometrios  to 
verify  the  radiance  model, 

An  additional  note  should  he  made  here  on  the  calculation  of  transmittance. 
Although  the  code  will  calculate  total  1 ransmittanee  for  a riven  atmospheric  path 
in  either  mode  of  prop. ram  execution,  the  time  is  increased  by  a factor  ol  N in  the 
radiance  mode,  where  N is  the  number  of  atmospheric  layers  alotu;  a piven  path. 
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\ppendix  A 

Listinu  ot  Piograni  and  Data 


A ! i .s 1 1 1 g of  the  Po  rt  ran  prog  re  m 1 .OW  T K AN  -1  1 1 ' l{<  hi  [f  A i\l  1 .OW  K M I is  pp  veil 
m Table  A 1 together  with  the  four  .subroutines  PA  I'll,  ItNOd,  POINT,  and  ANiid!.. 

The  input  ihitii  lor  the  |uw,r.im  is  pnen  it;  Table  All.  A general  flow  eliurt  for  the 
main  program  is  presented  in  Appendix  P,  atnl  definitions  ol  the  symbols  used  in 
the  computer  codes  are  summarized  in  Appendix  C. 

The  subroutine  1‘OIN  T has  ,t  twofold  purpose.  When  the  subroutine  is  called 
for  a given  altitude  N,  it  is  used  to  determine  the  mean  refractive  index  ( 1 ) and  in 
the  layer  between  \ and  the  level  above.  T\(b),  and  i:‘l  in  the  lay<  r between  \ and 
the  level  below,  AN.  In  addition,  an  interpolation  scheme  is  used  to  determine  the 
effective  absorber  amounts  per  km  at  altitude  \ for  each  absorber.  When  the 
parameter  If  is  set  equal  to  zero,  only  the  mean  refractive  index  above  and  below 
altitude  \ is  detorniiuod  from  ft  11  NT. 

The  subroutine  ANtll  is  used  solely  for  the  purpose  of  calculating,  the  initial 

zenith  angle  ((I  or  ANtll. ) by  an  iterative  scheme  taking  into  account  refraction, 

given  (1)  the  initial  and  final  altitudes  ol  the  path  (111  and  lid  respectively)  and  the 

angle  subtended  at  the  earth’s  center  tp  or  fl\  1' \1  1".  the  trajectory,  or  idl  the  initial 

altitude  and  tangent  height  till  and  lliillN  respect  i\ ely ),  A more  detailed  explanation 

-> 

of  subroutine  ANtll . is  given  in  Appendix  t'  of  the  l.OWTUAN  a report.  “ 

PA  I'll  and  lINOu  are  new  subroutines  which  have  been  added  lo  the  program. 
.Subroutine  PATH  is  used  to  determine  mid  store  the  cumulative  absorber  amounts 
through  the  layers  of  the  atmospheric  slant  path.  Subroutine  ilNO'.i  is  called  to  find 
the  nitric  acid  absorption  coefficients  as  a function  of  frequency. 


■ i i ■ 


vaararac  -i 


TIh>  last  column  in  tin-  listing  of  tin*  program  iti  Table  A1  is  tin-  corresponding 
listing  of  tin-  hOWTUAN  dH  computer  code.  Changes  in  the  previous  code  arc 
indicated  by  the  word,  NKW. 

It  should  he  noted  that  in  the  main  program,  card  l.OWlC  73110,  the  temporary 
fog  correction  has  been  commented.  The  fog  model  should  only  be  used  for  trans- 
mittance calculations  as  discussed  in  Appendix  D.  It  should  not  be  used  for  radiance 
calculations. 
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I’K'f'I.PftM  l PHI  4(  InPUT  » OUT  PU ! »T  APF  S » TAPTh-OUTPUT  ,T  API'  7 I LOWt 

common  M 3*»)  , »•  < 7 , 34 ) , 1 < 7 , *4)  « CH  ( 1 1 1 14)  ,WH(  7,  3 4 1 , M , NL  , Rl  , 0 W , 0 0,  P I L OMl 

0 I NSI  ON  WO  < 7,  A 4 I , M?  i < 3 4 1 , MTC<  h ) , AH  A3T  U4>  , AH7C  < ^ 0 ) LOHL 

1U  It  N SI  ON  'I’M  *.n  , » W <0  7>  » H / ( ?)  ,TX(tl)*VH(ll)  LOW* 

O I NSI  ON  HST  0^(14)  LOW 

common  tk  .■**•!) i;> , c.m  ,cx  < o i ,c.4  < 1 s,\)  ,cr(  ir>)  ,ca  ( x u / ) lpm 

COMMON  V*  (1.4)  ,f'7  (4rD  ,CT»  (44)  1.0 HI 

OlMINSION  PMIX(  <4)  , Hi  ft  Y ( U,  m I.OHI 

IMHI'NMO*  NPATh  <Sf.  ,1  n (Ml  LOW! 

< on  MON  Ml/  IfM  I SS  , kMA  X,  A NO  L f*  » U N , H M 1 M , T J , J 1 , J?  , JM  1 N , Jf  \ 1 WA  , IT  Y p»  LOW! 

0 0MMnN  /{**?/  T i , 1 kM  A\  , 10  NS  T ilP  , NIL  t W ( 1 l ) ♦ M 1 1 ) 10  Hi 

~OMMON  /r M.</  H 1 , H NP  1 , M00{  L LC’Ht 

HMiXI  IIOIMM  VMUI**I  Ml\]NI.  KA1I0*  T 1 m r q £*qu  FROM  LVANS  CHlHLl.  l OWl 

PAT  A HP  1y/Q*i).  *0  , t ,T  , A3  i 0.  ft  , t 1 . 4»  1 .f»,  l . H,  1 .0,  ?,  0 , P.  1,  ?.  J,  3.  I) » 3, LOWt 
J 7,  A,  t,  •0.0/  LO HI 

n (1  ,v>  --1 . (v**ro  /Ctxrd,  4 A87«i*v/i  )-j  . ) vow 

WATTS,  ('M-0  SF-t  MICRON- 1 LOHI 

r ( n ) * > l»  ( l /i , 1 7 »i  0 - 1 •#  . s*T ».  * A - r . 4 T ' * A * A > * A l 0 Wl 

XATA  Hf’|U/r>M.'j  K M/  , H/  ( 0 I / r>  H «»  KM/  LOHI 

L OH  T PAP  1 n»  f 77  LOW! 

I OWi 

AUT  HOI  S t OWl 

J.  1 . A-  S(  l HY,  l OH I 

f , Y .AN  r I 7YS  , lo  Hi 

J,  M.  r HI  thyuh  .W  • LOWt 

P.  A*  MCCIAlCHtY  LOHI 

LOHI 

PKp.-.PAM  LOW!  ^AN  CAlCUlATLS  IMfc  I'i  AN  SM  J T TANC  { ANP/OK  K'AOIANCl  LOW 
0 1 T Hr  ITMOSpm  pf  LOWt 

r p pm  rn-i  To  woonn  rn-i  To  m t f vons  i at  .‘j  om-i  iom 

SP'.LTPAl.  f?l  SOI  HI  T ON  ON  A L TNf  AM  NAVrNUM^FR  SCALt.  lOWt 

^IhHACTlON  AND  I A P T H roRVATO««  |MTDTS  API  INCLUOIO*  A T MOSPMI  Pt  LOH I 
IS  i AYf  t?rp  IN  OMt  KM.  INT'-^vM.S  OlTHpIN  0 ANO  PS  KM,,  «3  KM,  1NTik»-lOH| 
VAIS  Tp  »‘,r  <M , , A TWCNTY  KM , TNTf‘>VAt  TO  7 0 KM,,  ANO  A THIRTY  KM,  LOHl 
IMIi'VU  TO  150  KM,  LOHI 


P ‘JOUR AH  AOTTVATPP  HY  SpOMTSSlON  OF  I- OUR  CARO  SFOUENUL  AS  FOLLOWS  LOHI 

LOHI 

C ARII  1 MufU  L » I HA  Z I , ITYFl,LFN*  JP,  I M , M 1 , M , M 3 , ML  , HMISS,P0  ,T  HOUNO  LOHI 

I ORMAT  < 1 IT  3 ,.D*  1U,  3)  LOHI 

CAR1'  .*  Ml  »M.'f  ANI.L*  » WANC.I  ,01.1  A,  VIS  lOH^HJI  1U..O  LOWt 

('ART  -X  VI,  V.'*  HV  I ORMAT  C 71  1 n , ,U  t OWl 

CARO  4 IX  Y FOKMATU.xi  LOWt 

l OW 

MilOM  “ 1 , J,  4,  s OR  h SILt-OlS  ONI  nT  T Ml  I OLL  OH  T NO  MODI  L A 1 MOSPHt  Ri  l 0 HI 

TPOPIOAL.MIDL  All  HUH  S0M’1rR,  HIIU  AT  I 1 UOI  HI  NU  P.  SUH  AMD  1C  SUMhU  H»  LOWt 
SOIUPCT  IC  W1NUP,0R  T Ml  V\h?  II. S-  STANOARH  R FSP  FOT  1 VI  l v LOHl 

» pt?  H n p I 7 . path  wh>  n *«.' IF  oW(U  « CA  T&  VSfc  F * I N V 1 1 An  p»  CAk'P  :‘LOHI 
Rl  All  HI  ,1’C  MO)  , T ( nto  Oflirw  I*  T , T F HP  < 0 f 0 C ) , v.  R F L HUM  IP  I T Y*  H CO  It  NSITLOHt 


i n 
CO 
^ 0 
ml 
*,,1 
M) 
7 0 
rtf- 
"R  u 
10  0 
HU 
1 CO 
l AO 
JAO 
lr>0 

i U V 
1 .*li 
\ ,Hl, 

ID  C 
COO 
CIO 
C CO 
C.10 
Cm1 
C s c 

CuO 
C 7 0 
CaO 
C‘lu 
TOO 
J 10 

1 *>  J 

3 J II 
3 4 0 
.5  SO 
Tin' 

3 7i 
3U0 
391' 
400 

4 10 
4CV 

4 m 

4 4y 
4 ‘,0 
4 0 0 
4 70 

4 fit) 
4 l!  0 
‘•0  0 


Nl  W 
Nl  H 
A J* 
N t W 
Nl  W 
Nl  W 
Nl  W 
Nl  W 
N‘  H 
N:  W 
Nl  W 
N t W 
Nl  w 
Nl  W 
Nf  W 
Nl  H 
Nl  W 
A SC. 
A P 
A 7 
Nt  H 
N*  W 
Nf  K 
Nl  W 
Nl  W 
Nl  4 
Nl  W 
Nl  H 
Nl  W 
Nl  W 
A O 
A 1 C 

A 1 1 

A J C 

A A J 

A t 4 

A 1 ‘> 

A \ i. 

A 1 l 

Ni  H 
Nt  W 
A 1 M 

A CO 

A C 1 

A CC 

A ’ 3 

A C 4 

A C‘. 

a ,'♦•* 

A C7* 


•*)  n 


c 

c 

c 

r 

c 

r 

C 

c 

c 

c 

c 

r 

c 

c 


c 

c 

c 

c 

r* 

c 

c 

r 

c 

r 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 

c 


PHIS  PAGE  IS  BEST  QUAY, TTY  FRACTtCA&LS 
OlXPX  E.OKNlblliUJ  XV  W>Q  


'i 'all U1  /\  1.  I •istint;  nl  Kortcun  Coilv  1 ,OWT U A N <1 

<r-N  .M-J)  , 0 J riFNS  JT  Y(GH.  N-7)  , V IS < KM)  , RA KGE < KH1  WITH  FORMAT  4?S. 

LO  Wf 

MO 

A 

ra* 

NO  Of  L =7  HHFN  N t'H  MODEL  4 1 NOS P HERE  (E. . G.  PADlOSONOL  DATA)  USED. 

LOME 

‘i  30 

A 

?9A 

DATA  CAROS  ARE  RE  4 0 IN  IPTHFEN  CARDS  1 ANT  ?,  AND  SHOULD  CONlAINI 

LOWt 

JO 

A 

3 9n 

ALT  TTIIOF  (KM.  ) .PRFSSURf  . ICMP.CFH  PI .TFhP.REL . t-UM  III  I 1 Y , HDO  OENSITYjLOMF 

St  0 

A 

rsc 

0?  OF  NS  IT  Y , A FPOSOL  NO.  OEnS  1 1 Y ( C H-3 ) ACCORDING  TO  FORMAT  i*  ?9  • 

LOWE 

5SQ 

A 

39U 

NOT  f THAT  FITHER  Ot'W  pT,  TEMP. OR  RrL . HUMIUHY  CAN  BE  USED. 

LOWE 

‘ibO 

A 

3 St 

LDWt 

STO 

A 

?9F 

ARF  USED  To  CHANGE  TFMP,H20,  AND  0?  ALTITUDE  PROFILES. 

LOME 

SM0 

A 

3 96 

1FNIS$=0  = TPAN«»mjSS  JON  MODE  ( IE  MI  S S*  1 = L M I S S I ON  HOOF 

LO  KF 

S 90 

Nt  w 

TOOUNOsl FMPFRATORF  OF  FAKTH  IN  DEGREES  KELVIN 

LOWt 

MOO 

NEW 

IF  T BOOK'D  = 7ER0 1 ASSUMES  AIK  TFMPEKATURf  OF  M(JD  f L ATMOS. 

LOWf 

U13 

Nt  H 

LO  HF 

6 30 

A 

30 

IF  THA?r=o  NO  AEROSOL  SCATTERING  IS  COHFUTFP 

LOW? 

6 30 

A 

31 

1HA/F  =1  Ic  AEROSOL  AfTcNMATION  RtOUIRFO  (THIS  IS  USfcO  IN 

Lour 

640 

A 

33 

CONJUNCTION  XT  IH  VISUAL  RANC.E(SFE  card  ,’)) 

LOWt 

hSO 

A 

3 3 

iha7f  = i or  ? also  gtvt  atposol  attenuation  for  ?3km  anu  skm  vis. 

L (Y  WF 

6 60 

A 

34 

HA7F  Mpprts  Wc  SFFCTIVlLY  X F VIS  =0  ON  CARO  ? 

LOVJF 

6 70 

A 

3 G 

IHa;t  ~ 7 FOR  OTHER  AfrtflSOL  MODELS  U.C.,  MftRITlHF  LCT>  WHICH 

LOwF 

f»an 

A 

39  ♦ 

Akf  Rf A 0 INTO  PROGRAM 

lore 

690 

A 

3r*t 

L 0 HE 

70^ 

A 

3 6 

ITVPF.  = 1«?  op  j INDICATES  thc  type  OF  AIMOSPHTRIC  path 

LOWE 

710 

A 

37 

I TVPFe?  tVFRT  KAL  OP  StftNT  PATH  TO  SPftCF 

lo  wr 

A 

3 8 

iTVPF-r  »Vf  RT  IPAL  OR  SLANT  PATH  nriWtFN  TWO  ALTITUDES 

LOHf. 

A 

39 

I T Y FF  =1  * CORRESPONDS  TU  4 HORIZONTAL  (CONSTANT  PPlSSURD  PATH 

L OWL 

7 4i. 

A 

1 

LO  Ml 

7 00 

A 

41 

HlrOOSfTVFR  ALTITUnt  (KM) 

LOHF 

7 bO 

A 

43 

H.-'sjnUPr.F  ATTITUDE  ( <M) 

LOht 

/7U 

A 

4 3 

ANGLFs  7 F M I T H ANGl  F AT  HI  < D t F, R E h S ) 

LOWt 

7 ft  U 

A 

44 

RANGE- T ATM  LENGTH  (Km) 

l OWf 

?9C 

A 

4 6 

ftrr  A = F APT  H CENTRE  A Hr.  lF 

LOWT 

5 0 J 

A 

4t> 

VI'  s VISUAL  RANGE  Af  SfA  LFVEL  (KM ) 

LOWt 

bin 

A 

4 7 

( I R I TY  Pr  s % READ  HI  A NC)  K4  NGET  I F I T Y P fc>  ? RF  AD  HI  AND  ANGLE. 

LOWE 

A 

4 d 

IF  lTTri=C  READ  Ml  AND  THO  (HHtR  “AVAKLDKS  e ,b  . HD  ANT  ANGLE) 

LOHF 

o 3 0 

A 

49 

LOHF 

640 

A 

GO 

V 1 s INI  Y TAL  FVOUFNFY  < HA Vf NHMHt F CM-1  ) 1NTFQFR  V A LIU 

LOWf 

vH‘30 

A 

•>1 

V?=FINAL  FPEOUFMCY  (WAVTNHMREH  CH-1  ) I NT  f Of  P VAUJ‘ 

LOHL 

A 

r.  7 

pv=  frequency  intervals  at  which  transmittance  is  printfo 

LCWt 

*■70 

A 

b 3 

N DT  t t DV  MUST  Of  A NULTI»LF  OF  «i  CM-1 

LOWF 

6,-)0 

A 

64 

LOWF 

690 

A 

4 

I V Y = 0 TO  GNU  DATA  ,=  1 FOR  NEW  V1,V?,0V  ONLY  , a.-*  TO  CONTINUE  DATALOWf 

900 

A 

9 6 

Iky  = 3 ro^  new  card  p only,  =4  for  nfw  CARn  i only. 

LOWF 

910 

A 

6 7 A 

LOHF 

S3I3 

A 

*T/tt 

HHIX(9) =HMIK ( ?P) = 1 . 0 f.KO 

LOWF 

930 

NE  W 

IKY:  0 

LOHf 

94  0 

A 

6 7C 

*NAX  = ii 

LllWt 

9r»u 

Nf  W 

READ  (0,4  0 0)  I ft  TM  , NL 

LOWF 

Sb(l 

A 

9 8 

KT  All  (>■  ,1,01)  ( H 7 1 (X  t , I *1  .ML  » 

LOWF 

970 

A 

G9 

RTAIT  ( S i 4D  1)  (M?2  ( I)  , Tsl  ,r> ) 

LOWf 

9 d 0 

A 

60 

h;2<m  = hzi  (n) 

LOWF 

990 

A 

6 0 t 

DO  1 J«I,3 

LOWFlftUO 

A 

M 

>*»wh^|.i  :$+-.  -n;  ■*•• 


THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
FROM  CO!  Y FURR  I SUED  TO  DDC  


TabW1  A 1 . l.intinir  ot  Fortran  1'oilr  I.OWTHAN  4 (Cont) 


LOHf  J 0 10 

A 

b? 

K1=K?-1 

LO*F 10?0 

A 

Gv3 

no  i i*i,nl 

l OHF  1 G 30 

A 

b4 

i 

RT  AO  <r-,402)  ?U>  , (PCK,  I)  ,T  IK,  1)  ,HMIK*  II  ,W0(K,  I)  , K=K1,  K?) 

LOWL 1040 

A 

REAP  (*,43i>  (V*(I  > ,C7<I>  ,C7MT)  ,1  = 1,44) 

LOHF  10  SO 

A 

b G* 

REAP  1^,4031  (IPU)  ,FH(I)  ,roai  ,1  = 1, 67) 

LOHLlOGQ 

A 

6/ 

RT  AP  (6,404)  iCl  ( I)  ,1=1  , ) 

LOWE  1070 

A 

G* 

REAP  16,404)  ( C ? 1 1)  *1  = 1,  ) f>  7 r« ) 

L OHF 1 0A 0 

A 

PS 

REAL  (5,404)  ( C 3 ( I ) *1=1,640) 

LOHF  Hi  HQ 

A 

7 H 

REAP  (6,405)  (C4<I)  ,Ixi,  i33) 

LOWE  110  0 

A 

t\ 

REAP  (6,404)  (Hri  < I)  ,!  *1  * 1 *3  > 

L O wt  1 1 i o 

A 

7 c 

RE AO  («-,  40r>)  <GM  !>  , 1=1 1 10^*  1 

LOHt.U?0 

A 

?3 

JST  t‘R=P 

LOHF  1 14  Ij 

NF  W 

PI*?,0*ASIM(  1,0) 

LOME  1 1 4 r 

A 

7 4» 

CAsPI/lflOt 

LOW!  IlhO 

A 

7 4 

TP=0 

LQWF  u?o 

A 

7 b 

2 

nONlTMtir 

L OWL  11 ti  0 

A 

7 7 

Vt*b37l  ■ 5 3 

LOHt  1 19  J 

A 

1 F lK-n  = p 

L0M  1700 

A 

79 

C 

Ju  NF  C SUPRfSS  PRINT 

Lourin'* 

A 

REAP  4 0 0,‘1Or>n.,lHft7f  , lTY*>cf  if.n,  JP,  IM,  1EM1SS,  RO,TBCUNO 

L0  4f  1 JL-C 

M:'  w 

0 

iemiss-=o=t RANsrt i ssi on  Pont  7 imi'ism-f mission  mope 

LOHl  1 ? 3 0 

Ni  W 

I F ( I E P 1 SS , F O • 1 ) PRINT  11  M 

U>Hi  1P40 

Nr  U 

IFUrniSS.  rQ»  0)  PRINT  HM 

LOHf  lPA*’ 

Ml  H 

lEN^TOE=LFN 

l OHF  1 r b G 

Ni  W 

PRTNT40&, MOOlI ,1 HA  7f  , I TYPE ,LfcN, JP, 1M,H1 , NT, M3, ML , If MISS, RO,T  UUUNO 

L 0 KHU 

N(  W 

:oo 

H sHOn^l 

LOHF  1 ?o  U 

A 

8,’ 

IF  (M,t  0.  1 ) Pr  = f>37A  , IQ 

L O H(  1?  )0 

' A 

8 3 

IF  (M.f0.4)  QF=P35fc.«ll 

LOHF  1^00 

A 

v84 

IF  ( M.f  O.  r>>  RF=b36b. 91 

LOWf  1310 

A 

8‘-A 

IF  (IMA  7F  ,t|f  .7)  G'J  TO  2r-0 

LOStA3pu 

A 

6b-J 

REAP  4?1,  (VX(1),0  7(1)  ,r7A(H  ,1=1,4  4) 

LOW t UJ0 

A 

»‘>C 

PRINT  4.U,  CVX(T)  ,07  (I)  ,C7A(I)  ,1=1,44) 

L UHL  1 34  0 

A 

br-D 

T H A T F = 1 

LOHf  U5i* 

A 

I F < PO  « f.  T • 0 ) R E = R 0 

LOHf  1 3b  0 

A 

\*,r 

I F ( M, FP . 7 . AN O.tM *Nt  .0)  GO  TO  4 

LOHF 1 370 

A 

'« i,G 

IF<  TXY,GT,3)  GO  10  ? 

LOHfc 1380 

A 

• »;>< 

IF  CHOP EL*  El*  0)  GO  To  4 

LOHF  1390 

A 

tt- 

300 

RF  A 0 406*  H1,H7,  ANGir, RANGF, HKT A ,VI^ 

L OWE  1 4 o 0 

A 

b7* 

PRINT  4?5,  *1, HP, ANGLE, °/\NGr,HETA  ,VIS 

LOWl 1410 

A 

X1  = PF  FHl 

LOHf  1 4 ? 0 

A 

89 

IF  UTYPF.FQ.3)  GO  To  560 

LOHF  14J0 

A 

90* 

IF  CITVPF.EO.U  GO  TO  * 

LOHF 1440 

A 

41 

XP=PE+H? 

LOHf  145U 

A 

h: 

IF  (PANff.n  .0  .)  GO  TO  6 

LOW)  1 4 b J 

A 

M3 

PRINT  4?fl,  Hi  ,MP  , ANGLf  » RANGE  , PFT  A ,V  IS 

LOHF  14  7 0 

A 

M4 

IF  (HO. FO.O* AND. ANGLE. NE, 0)  GO  IP  3 

LOHL  1 4 d C 

A 

1J5 

ANGLE  =A COS  l 0 *5M  . ♦XT/ XI ) / RANGE -KANG E /T  1)  ) / CA 

LONt  1490 

A 

G C TO  7 

LOHf 15u0 

A 

9 7 

IwjUaiiiLja***,* utl»^-ii  iif'liilirii-'artMiYiil  Iilif  'ii  >fliril^h-flk vj^is 


son 


r>i*o 


fRl S PACE  IS  BJ.ST  QUALITY  PBACTICJlBM 

ftW#  WTX  fTJl'vNliiu^u  10  DUG 


Tabic  Al.  Listing  of  Fortran  I'oiic  LOWTUAN  ■>  (t  out) 


X ?=5QRT  < < X 1 / P ANGF  ♦ R ANGr  /XI  *?.  ,Q*C0S  IANGLE*GA)  )*X1*KANGE  1 

10  41  IMP 

A 

MM 

H 2 = X 2 -F*  E 

tout  lr.  JO 

A 

HM 

CO  TO  7 

LOUt  1 1 * S 4 

A 

i v : 

CONTINUE 

LOW*  1^40 

A 

10  1^ 

IF(KL.Lf  . (1)  NL*1 

LOHf 1550 

A 

lCv* 

no  540  K = 1,ML 

LOHE  lr»bb 

A 

1 0 3A 

AH47E  (K)  =0,(1 

LOW!  1570 

A 

n-3H 

IF(M,EO,Q  > RE  AO  4?<)>H1  ,P<7,1>  .THP.HP.RH,  WH!?,  *>,N0<  7.K)  »VIS#rANGf; 

LOME  l‘id0 

A 

n jc 

IF <H,EP.0> PRINT  4JD9H1,P(7»1)  vTHPfOP«RHfHH<7.  "I  ftfn(7fK)  .VI$vRAngE 

lone  ir,sn 

A 

1 P JO 

IF  <MtGT , C) RE  AO  4 7R, 7 C K) f°  C ?,K) ,T  HP  ,0  P , RN  , WM  ( 7 ,K*  f NO <7, K)  , AHA  7 f U) 

LOHf  ll. OP 

A 

1 13f 

IF <N,EO*n>7<K>»Hl 

LOME 1510 

a* io jr 

J= IFI>  < 7< K> *1 ,0£-M  *1  . 

LOUtthLu 

A* 

1 4 Jl» 

IF  (7  IK)  .GF.2  5.  t)>  J«  <Z  <K) -?5, 4> /*«  0+?6* 

LOUt  lfvTp 

A 

1 I-  JH 

ip(74k».gf.5‘).q)  jMzm-so.n/zo.j  + si. 

LOWE  li.4 D 

A 

11  J1 

TF<7<K|  »CF  .TQ«0)  J*  (Z  <K  1-70.0)/  3 0.1  + 3?, 

LOME lb50 

A 

1 0 JJ 

IF< J.OT.33J J»53 

LOWfc  Uujw 

A 

lv  TK 

FAC*Z (K)- FLOAT (J«l) 

E OWE  lb To 

A 

1 l»  3L 

IF(JtlT.?FI  GO  TO  500 

LOUf iOdl 

A 

i r jrt 

TAC* (7(K) -G.Q^FLOAT  < J-2fc> -?5.>/*5, 

LOW! 

A 

U?N 

IFIJ.Gf.M)  FAC= <7<K> -50.01/20. 

LOHE  1700 

A 

1 G JO 

I F < J , GE .3  2)  F«C*  <ZtM -70.0) /J9« 

LONE  1/10 

A 

1 0 JP 

IF<FAC.GT , l, 0)  FACsi.O 

LOHt  1 7 J li 

A 

i ja 

l-J*l 

LOW? 1 7 JO 

A 

1 C 3*t 

T ( ?,KI»TMPf?7 *. 

LOW  1 74G 

A 

It  JS 

IF  (Nl.CT.Q)T  ( 7t  K)  *T  <M1,  J)  ^(T(Nl  VU  / T (Plt  JM  v»FAC 

LC1WE  175  0 

A 

ll  J1 

T Ts?73. lS/Tt 7,K) 

LOME  I7b^ 

A 

1 - JU 

If  (»'H.U  .0  .0  ) TT07  3 . 15/  < ?7  j.  imqp) 

LOU!  17  70 

A 

It  'V 

IF<ViH(7fK),LF.0.'i)  VH  (7  , O = r(TT) 

10  Wf  1 7.,  j 

A 

1 ju 

TM*2.CT,0)wn<7,k>  - PH  (*Vt  .))  * ) / uu(  m(i(  j>  j **rAr 

L t • W , 1 /m 

A 

ir('H,L,T,  0,1)  H»>(  7t  It)  =c.  1 (7,*() 

L'1  Hi  ).  j i 

A 

liU 

IF(M3,r>T.riWO<7,*C>*KMM7,  J)»(WrtCHf  LI/WCCMJ,  J>>  **FAC 

LOUt  lc  1 

\ 

1 V 3 7 

HFT  OF ( K)  = ?, 

1.0  Ht  18.\ 

N W 

IF(»*HIK  (Jl-LF.  D.l  CO  TO  5 V 

1041  1 F 30 

Nt  w 

MFTOf*  J)  MHNIX  (L)  /HM  I X ( J } } **FAC 

LOWE lt4C 

N 1 H 

CONTINUE 

Ei' 4i  ior»^ 

Nf  W 

I F ( 2 (<) .CE .5 . 0) GO  TO  5?0 

LOHt  1 PbO 

A 

1 if  4 A 

lF(AHA?t(K)'t'0«0.0)  *N2?(K>sH72(J)MHZ?(l)/H7;MJ)  )**FAC 

LOHf  lb  70 

A 

1C  4 1 

IF(AH47c(K).FT,C,0>AHA7E(»O=HZl  (J)  MHZl(L)/H7t  (j)  ) * *F AC 
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NEW 

GO  TO  1024 

LO WE3940 

ME  M 

1 0 22 

N < K ) * f V 

LOhf  39  50 

urw 

iFtji.Nr.jr*  r.o  to  to?*. 

LONE  3 9b0 

Nt  H 

,K>  »W(K! 

LOME  3970 

NEW 

W(K)C||. 

LOWE  3980 

NEW 

Jf  X TR  fl*  l 

LOME  3990 

NE  W 

10  3*, 

CONTINUE 

LOME  40  00 

NEW 

THIS  PAGE  IS  BEST  QUALITY  PIUCIICABUI 
J^OttGOO 


Table  Al.  Listing  ol  I'ortran  Lode  l.OWTHAN  4 (Cont) 


If  (JP.n.Q)  PRINT  435,  I,X1,  <VH(L)f  L=if6>  ,PSl*PHT,af1At  THETA.SR 

LOWE  40  10 

A 

2 b 2* 

IP  (I  . C- E *Nl  ► GO  TO  rr> 

L C H f 4 0 2 0 

A 

2b  3 

If  <I*1.FQ,J2)  FH(q  »I+i(  = YN 

LOhl-UG  3 0 

A 

2 04 

n u.fq.jii  eH<q,i)*Txj 

LOHfc  40  4 0 

A 

?b5 

RN  = fH(S.lM)/FH<9,j) 

LGWF4050 

A 

?b<' 

SPHI  = SPh7  *px/rn 

LCWf  4v.G0 

A 

?17 

IF  <SALF.Gf.RN)  5PHI=?ALP 

'.OHf -.070 

A 

2F6 

25 

CONTI  Niff 

LQWt  4 0 6u 

A 

2t?  9 

GO  TO  47 

LOHL4090 

A 

2 70 

c **  ** 

horizontal  path 

LOWf  4 LOO 

A 

27  1 

-?6 

OO  ??  Ksl.KHAX 

LOHE4H  0 

NE  W 

H<K)sPAN&F*f H<K, 1) 

LOWE  41 20 

A 

2 73* 

IF  (MODEL. CT.P)  W(k) aPANC»f*TX<K) 

LOME  4130 

A 

2 7 4* 

V H < K>  = W (K) 

LORE  4 140 

NEW 

77 

CONTINUE 

LOWF 4190 

A 

2 75 

GC  TO  4 9 

LOWF  41  GO 

A 

:7b 

:*8 

CONTINUE 

LOU6  41 70 

A 

27  7 

r * * *» 

DOWNWARD  TRAJLCTnPY 

L OHt 4160 

A 

2 7 6 

r 

LO  RL  4l  S», 

A 

279 

IF  (NP1  .FQ.1  > J 1 =Jl " i 

L0HF42U v 

A 

230 

J?  = JIM 

LOME  4210 

k 

>M 

J'JIM 

LOWF422 o 

A 

763 

YN 1 = YN 

LOWL4230 

A 

76? 

IF  ( M?.  GT*7(J1*1)  * pff*  HI . CO. H?)  GO  TO  30 

LOWf  4240 

A 

2 8 4 

IF  <NF1  .EO.l  .ANO.N?  ,GE.7  tJlMM  GO  TO  30 

LOWE  42l>  0 

A 

76  5 

CALL  POINT  <Hr,YN,N*NF?f TX,Tp) 

LOME 4j b0 

A 

26b 

no  ?9  K=1,KMAX 

LOWE  42  7 0 

NE  H 

29 

N ( K ) - T X ( K ) 

LOHfc  42  6 0 

7 66 

TX2tTx<  9) 

L OH f 4290 

4 

26  9 

YK?tYN 

iowe:43oo 

A 

793 

IF  <M?,  lT.HI)  H*H? 

LOWf 4310 

A 

7 9 1 

J?-N 

LOWf  4?  7C 

A 

292 

IF  iJl-ifG.J?)  TY?=7X1FYN?-E|MT,N> 

LOnE433C 

A 

79  3 

IF  {H3.GT.HD  TXL=r>? 

LOWfc.4340 

A 

294 

IF  ( Jl» fQ,j3 , ANO.M? .LT.MI)  YN1=TX2 

LOVE  43*50 

A 

2 95 

30 

A 0=  (RF+ HD  •SPHDYNi 

LOWf  4361 

A 

2 9b 

IF  lH2.GE.Hl)  XN2-YN1 

LOWE  43  7 0 

A 

74  7 

no  n T = 1 ,J1 

LOWE  43  60 

A 

2 98 

HMI  N« AG/EH<9 , J) -»t 

LOWF43N0 

A 

299 

IF  (I.EO.JJ)  HMIN^flO/YNl-RF 

L O WE  44  0 0 

A 

30  0 

JMIN*I 

LOWf  44 1 Q 

A 

301 

IF  (HMlN.lEi  Z<IU>>  l»0  TO  32 

LOWf  4420 

A 

307 

31 

CONTINUF 

LOWE  4430 

A 

3 03 

3? 

X=HMIN 

LO  Wf 4440 

A 

30  4 

IF  (HhlN.LF.O)  GO  TO  3U 

LOWE  445  0 

A 

305 

CALL  PO IN  v < X,YN,N,NPfT> , IP) 

LOWf  4460 

A 

30b 

J M I K “ N 

L OWL  44  7 0 

A 

3 07 

T XJ  = T X(  9) 

LOHfc  4460 

A 

3 U 8 

IF  (J2.  EO.N.  OR.  Jl.ETJ.  N)  T YN2  + TX  <9  ) -EH  ( 9 , N) 

LOWf 4490 

A 

30  9 

IF<  TX3.LT.0.  0)  TX3-T  *<  9) 

LOWfc  4500 

A 

3 0 9» 

affiiSPAa£is  bestiality  practicable; 

pggx'flutolb'iicui  10  W),c  — ’ 


Table  Al.  Listing  of  Fortran  Coilr  FOWTHAN  4 (Coni) 


IF  <Jl.tQ.N*ANn»H2.GF  .HI)  00  TO  31 
VHI  N=  AO  /T  X3- PE 

If  .GT.Q  .PH  0 1»  00  TO  3? 

33  'll’  <J1.EQ.N,  AN0.H2.  Cf  ,HD  YN1=TX* 

IK  < J2.EG.N.  ANtT-Jl.NE.tJ2>-  YN2=TX$ 

IK  TX?p;(X3 

IF  .( H2  • OF  * HI  ) J?*N‘  V 

If  fM?»GE .Hl.GR.H7. LT.HHIN)  H=HMIN 

PRIM'  436,  HMIN 

If  <H?«LT.HHIN)  J?=N 

If  (H?h.lT  ,HMIN)  M«;TNT  440,  HMIN 

GO  TO  !f 

34  PRINT  VJ6X,  HMIN 

IK  <h?.lt;hd  go  to  is 

IF  <ITVPF,tQ.3.0R.H2,GE.Hi>  PRINT  437 
ITYPE *2 
TX?=EH<9, 1) 

J M T M=  0 
02  = 1 
N2=  0.  0 
H = 0 .0 

C****  NOW  DFFTNF  V FRT iC At  PATH  QUANTITIES  VH(l-8) 

35  IF  C JF.  EQ.  C ) PRINT  4?Q 

JST0R=0-1 

00  40  1=1, NL 
J = J-1 

PFF=FH(C?J) 

if  u.FQ.n  ref=yni 

IF  <1  .FQ. 1. AN0.K2.E 0. 1)  Rff=YN2 

IF  < J.FQ.  J2.  AN0.K2.  EQ.O  REF  =T  X ? 

IF  < T »NE  * 1 ) Xl  = ?Un> 

X 2=  7 ( 0) 

IF  (J.EQ.J?.  ANO.K?.  FQ.O)  T2-H 
IF  (J.FO.JMIN. AN0.K2.E0. 1J  X2=HHIN 
HN=(RE*X1) ^SPHl-RF 

IF  <HM,GT.7(J>  .ANO.HM.GT.X?)  X?=HH 
RX  - ( RE*  XI)  7<RF*X2> 

OS=  X1-X2 
ALP=90. 0 

TKFT=ASIN (SpHI>/CA 
S ALP*  RX*$PH  I 

IF  (AES  (X2-HH) .GT-1 .OE-5)  ALP= ASIN<3ALP)/C A 
PFT^ALF-THET 

IF  (SPHI.GT.i.OE-10)  OS= (R£*X2)*SIN1BET*CA)/SPHI 

ThFTA  = 1 80.0-THF.T 

BET  A=8ETA*RET 

PSI=BET A- ALP* ANGLE ♦180#  0 

SR=  SR  +PS 

OO  1039  K = 1 , KMA  X 
AJafH(K,J) 


L0W*  4 510 
LOWE  4520 
LOWE  45  30 
LOW*  4540 
LOW*  4554 
LOWL4500 
LOWf  4570 
L 0 WE  45  6 0 
LOWf  4590 
LOWE  46  0 C 
LOWF  4b 1 0 
L 0 WL  462  0 
1 0 WE  4h  3Q 
LOW*  4b 4u 
LO  HE  46^0 
LOWf  46b0 
LO wf 4670 
LOWF  466  0 
LOWf  ho  9$ 
LOWf  47  00 
LO Wf 4710 
L 0 WE  4 7 2 0 
LOWf  4730 
LOWf  4740 
LOWf 4750 
LOWt  4760 
LOWE  4770 
LO WE4780 
LOWF  4 790 
LOWF  48  00 
LOWF4610 
LOWF  4820 
LOWf  4 8 30 
LO  WF  4o  4 o 
LOWF 48  5L 
LOWF 4 860 
LOWF4870 
LOWF  48tj0 
LO WF4O90 
LOWF 49  0U 
LOWE  4910 
L 0 WF  49  2 0 
LOHE4930 
LOWF  494  j 
LOWF  49  5 0 
LOWE4960 
L0WE497D 
L OWE  49  8 C 
LOWE  4990 
LOWE5000 


A 310 
A 3ll 
A 312 
A 3 13 
A 314 
A 315 
A 31b 
A 317 
A 318 
A 318* 
A 319 
A 3 2 0 
A 321 
A 322 
A 32  3 
A 324 
A 325 
A 32  6 
A J?/ 

A 328 
A 329 
A 330 
A 3 31* 
Nt  w 
A 33? 
A 3 33 
A 3 34 
A 335 
A 3 36 
A 337 
A 338 
A 339 
6 340 
A 341 
A 342 
A 343 
A 3 4 4 
A 34  5 
A 346 
A 347 
A 34b 
A 349 
A 350 
A 351 
A 352 
A 353 
A 3 54 
A 355 
NEW 
A 357 


54 


**>  *kv-sVj*^  •'  ■*  ' . « .!>  ■r*:A"b$ljtfl?1~-  f 


THIS  PACE  IS.PE5T  QUALITY  PRACTICABLE 
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Tahir  Al.  Listing  ot  Fortran  (’otio 

l ,OWT KAN  <1 

(Oont) 

nj*fH<K* j*d 

LOHE501U 

A 

350 

IF  (J.FO.Jl)  9J*F1K) 

LOWE502G 

A 

369 

IF  (J.F  Q,  J?.ANn«H?fc  IT  •H1.AND»H?«GT.0*Q)  AJ=jW<K) 

LOHt  5 0 30 

A 

360 

IF  (J.fcC. JHlN.ANn.HP.Cf «Hi)  4 J*TX (K ) 

LOWf  5U40 

A 

3 61 

IF  U,FQ.  JMIN.ANO.AnS<H?-HM)  .LT.U0E-5)  AJ*T><k> 

LOHF60S0 

A 

362 

IF  <K?.EQ.C>  GO  TO  36 

L 0 WF  50  60 

A 

363 

IF  9J=M(KJ 

LOWIGtf 70 

A 

364 

IF  tJ.rC.JMlNI  AJ*TX(K> 

LOHt  50  60 

A 

3 65 

36 

IF  1 AJ,  FQ.  0«  0*  OR«flJ  .f  0,0  , 0)  GO  TO  33 

LOWf 5090 

A 

3 66 

IF  ( AJ.  FQ.*U>  GO  TO  37 

LONF51QQ 

A 

367 

f V=  OS*  5 Aj-BJ)  /AlOG  ( AJ/nj) 

LOWE  5110 

A 

360 

G C T O 3 9 

LOHE  5120 

A 

3 69 

37 

LV=OS*AJ 

LOME  5130 

A 

370 

GO  TO  39 

L0MF514G 

A 

371 

37 

tV*CiQ 

LOHF5150 

A 

37  2 

39 

VH<K)^VH<K»  *-F  V 

LOWE  5160 

A 

373 

10.59 

HI  A Y<  Jf  Kl  * F V 

LO  ME  6 170 

NF  W 

IF  (JP.EU.0)  PRINT  4?Sf  J,X1,  L=l»8)  < PSI,  ALP, 

etT  A,  T HfcT  A > SR 

LOWF518Q 

A 

3 7 4* 

IF  (J*f Q*JP. AND«H2,GF.Hlt  GO  TO  45 

LOME  5190 

A 

375 

IF  <J  .t' Q.JMI  N»  AND  .K?.t  P.  1 ) GO  TO  43 

LOME  52  00 

A 

376 

IF  TJ.NE.ll  RN=RFF/tH(9, J-l> 

LOHC5210 

A 

377 

IF  (J.fQ.JrH)  RN  *RF  F /T  X? 

LOWE  5220 

A 

3 78 

IF  IJ.f Q.J2. AND .K2,fO#0>  RN=RLF/YN? 

LOWE  5230 

A 

3 79 

IF  < J.FC.  ( JUNfl)  *A60.K2  ,E0.1)  RN=RE*VTX3 

LOWf  6240 

A 

36  0 

IF  UALP.GF.RN)  9N=1.0 

LOW! 5250 

A 

301 

SRH  I=SAIP*PN 

LOWE  5260 

A 

362 

IF  <J  ,f  Q.J2,  ANfl.K?.€O,0>  GO  TO  41 

LOWE  52/0 

A 

383 

40 

COMIMl'F 

LOWf  5?fl 0 

A 

304 

41 

IF  <HMlN. LF . 0 ) GO  TO  47 

LOHt  52  90 

A 

3 84 

IF  (IFN*FO,0>  PRINT  438 

LOWE  5300 

A 

3 06 

IF  tLEN,EQ.0»  GO  TO  47 

LOWF  5310 

A 

387 

IF  (LFN,FQ. 1 ) PRINT  439 

LOHE  5320 

A 

386 

ii 

LOHt  53  30 

A 

3 86 

XI*  X? 

LOWE  53-40 

A 

3 90 

IF  UPS  5X1-HNIN)  .LF%0.not»  GO  TO  47 

LOWf  5350 

A 

391 

H*  HMI N 

LOWF  63  bi) 

A 

392 

J= J2+1 

l OWE  53  70 

A 

393 

IF  lNP2.FQ*n  J'J-1 

LOWE5360 

A 

394 

H = 9 F T A 

LONF5J9U 

A 

395 

PH=  l«0. O-ASIN(SPHt) /GA 

1 OHE5400 

A 

396 

TS= 

LOWF5410 

A 

34  7 

PS* PS  I 

LO WF542D 

A 

390 

OO  4?  K=1,KNAX 

LOWE  5430 

NEW 

42 

FT  K) (<» 

LOHt 5440 

A 

400 

GO  TO  39 

LOWE  5450 

A 

401 

43 

BET  A=2.  *BETA-n 

L OWE  5 4b 0 

A 

402 

PS  W»*-  FST-PS 

LOHt  54  70 

A 

403 

SR=?.*SR~TS 

LO  RE  54  00 

A 

404 

r 

LONG  'ATH  TAKEN 

L 0 Wfc  649  0 

A 

4C5 

PH  I=PN 

LOHt  5500 

A 

4 06 

JMS  PAQE  Xi  BEi.IT  QUALITY  FRACTlCABSiS 


Table  A l,  lasting  of  Koetr 


no  <.4  K*l, kh  a x 

44  VMM  W*.*  VH(K)-t  (K> 

CO  10  4 7 

4**  CIO  4h  KM.KNAX 

4b  VH<  K»*?.0*VH(K) 

WE  I A*  ?•  C * OF  Y A 

S R-  ?.  0*SR 

IF  CHr.fQ.Hl)  GO  10  4/ 

RNM  X1/VN1 
SPH  I * S I N f ftN0tt*CA) 

I f CSPHT.LT.RN)  SPH  l*  SPH I/RN 
GO  TO  11 
47  C ON  1 1 NOE 

IF CANGIE*GT.90. 0>  PRINT  40b. MM 

00  4B  * * 1 f KH  ftx 
WCK)*VH(K) 

40  CONTINUE 

49  WR1TF  <b,419) 

WRITE  fF.il  *5) 

WRIT?  <b,4m  CHU),I=i»a)tW(H>(WUil 
WH I TE (7 | X 0 0 0 I N , IHA7t ,1T  VPf ,H1, ANGLE » HM I N » 
.1000  FORMAT  C,U3.6Ftl  .4) 

1 3 1 
L-  1 

IVlsVl/5.0 
ivr=v?/9,+.9s 
lvi*5MVl 
IV 2***1  V? 

IF  1 1 V 1 »L  T , 34Q  > I VI -5  40 

if  uv?.gi ,'icjooo)  ivr*4oo)o 

IF  COV.LT.*.)  OV*r*. 

irv^ov 

iv«tvi-iov 

icoont*o 

I F < IF  HISS.  FT.  01  f.O  10  *0 
RA0SUN* 0. 0 
FACTOR* C.  S 

CAU  PATH  ( ML  AY  »W  PATH,  TOOT) 

PRINT  115b 
PRINT  1157 

if cuwiss.ro. A)  ikmax*iklo 

n**»*  DEGIHING  OF  1 RANS^  IT  T ANCr  CALCULATIONS 

go  iv=iv*mv 

SUMV=0. 

TLOLOM,  4 TSOLO*1. 

I KL  0= 1 

IF(  IFHISS.FQ.O)  I KL  Os  1 K M A X 
00  1050  IK*IKU).  IKMAX 
IFUFH1SS.FQ.0I  CO  TO  1055 
no  1052  K*1,KNAX 


Ciuie  l.OWTHAN  4 (i\ml) 


1.0  HI  ' 

•510 

NL  W 

LOW  » 

->  *»  ? 0 

A 

40  A 

l OWE 55  JO 

A 

404 

LOHt.  ‘>5  40 

NL  W 

LOW  554U 

A 

411 

l OKI  45  bU 

A 

4 1 2 

L 0 H l 4 5 7 0 

A 

413 

LO Hi  4 5 Hi) 

ft 

4 1 4 

LOW*  * 

»59ti 

A 

4 15 

Ll'Hl  M’OU 

A 

4 1 l« 

LOW  : 

>b  1 1) 

A 

41  7 

L UWt  5bJ0 

A 

41  tS 

LOW!  r 

•b  JO 

A 

4 19 

LOW)  ‘ 

>b4C 

A 

4 l *H 

LOW  '>br'0 

Nl  H 

l OHt  ‘*bhO 

A 

4 1 

LOW  u b 7 0 

A 

•* r 

L OHt  5b  60 

A 

4 ?s 

LOW  5090 

Ni  W 

LOHt  ‘ 

‘700 

Nt  H 

LOW  5?  10 

Nl  5 

L 0 Wf  1 

‘7E0 

Nl  H 

i OH)  ' 

‘>730 

A 

4?  4 

low  ' 

>740 

A 

4 2 b 

LOHt  ' 

5 7 415 

A 

•*ct 

LOHt  * 

A 

4 

LUH  t ' 

l>7/0 

A 

4 2 9 

LOW  1 

‘■740 

A 

‘•3  0 

LOHt  5 7 90 

A 

431 

LOHt  ' 

:>  an  o 

A 

4 32 

LOW  50  10 

A 

433 

LOHt  1 

48i*0 

A 

4 H 

LOW  « 

4 O 3 0 

A 

4 3 4 

l OHt  ' 

5^40 

A 

H 3 h 

LOHt  ‘ 

Nl  H 

LOHt  5*1*0 

Nt  W 

LOHj  5H  7 1) 

N t H 

l OHt  1 

Nt  H 

LOHt  5A90 

Nt  W 

LOHt  4 90 0 

Nt  H 

LOHt 

*•910 

Nt  W 

l ONf 59  0 0 

A 

437 

l OHt 

553  0 

A 

4 3 oA 

LOHt  5940 

Nl  H 

LOHt 

49  5 0 

Nt  H 

LOW 

5950 

Nt  H 

LOHt 

49  7 0 

Nt  N 

LOHt  4980 

Nt  W 

LOHt  59% 

Nt  H 

LOHt  bflOO 

Nt  H 

THIS  PARE  BEST  QUALITY  nSACTlCAl^' 
EEiAW  (iyi'X  EUiilUSlUiD  WUUC 

x • i-  ■ 


Jabli1  A',  l.istinj;  ot  Fortran  Codr  l.OWTHAN  1 b'ont) 


H ( K ] *WPATH( IK, K) 

LOWtbOlO 

Nf.  W 

1 0 5? 

CONT  I Nil F 

L O 60  ? 0 

NL  H 

iogg 

i j=  n: 

U)WFb030 

NF  w 

IF  < JP  «N f . U ) GO  TO  5? 

L O Ht  b 0 4 0 

A 

4 3 00 

if  (Icount ,ro.o>  go  to  si 

LOWFbl)  50 

A 

43  9 

ir  ( icmiNt .eo.^o)  oo  to  <h 

LONLbJbO 

A 

4 A 0 

GO  TO  T-? 

LOHF60  70 

A 

44  ] 

5i 

ICOUNTi 0 

LO  HF  60  6 0 

A 

44? 

1FUCMI  SS.EQ.O)  PRINT  4?? 

LOHfbOSO 

Nf  W 

5? 

no  53  K=1|KHAX 

LOH[ 6100 

Nt  H 

T X<  K) *0  * 0 

LOWt  GUO 

A 

445 

IF  (K  >l  Ta  4)  T X ( K ) = 1 . 0 

LOHEM20 

A 

44  b 

53 

rCNlINUE 

LOWFbl 30 

A 

44  7 

icounTs irouNTf i 

LOWf bl40 

A 

44  0 

5 U N s 0 * 0 

LOHF  bl 50 

A 

44  9 

V*  IV 

LOWf  M60 

A 

460 

1=  UV-36Q)/5M 

LOHf  6170 

A 

46  1 

C ** 

LOHL6160 

NEW 

C 

HN03  AOPfRPTION  CALCULATION 

LOWF  619  0 

NF  H 

CALI  HN03  <V,HAnS) 

L OWE  6?  0 0 

Nf  W 

TX< 11)=HAPS*H< 11 ) 

LOWf 6? 1 0 

Nf  H 

slm=suh*tx<  t i) 

LOWrp.120 

Nf  H 

IF  < IV  .IT,  670  > GO  3 0 7? 

LOWf b? |0 

* A 

4 60* 

IF(  I V.L  F .3000)  GO  TO  01 

L O HI  6?  4 U 

* A 

4 3* 

c * * • * 

»*»*  MOLECULAR  5 C A T T FRING 

L 0Htb£5u 

A 

454 

LOWI  6?  6 0 

A 

4r<6 

T X ( b)  (5) 

LOW  ft?  TO 

A 

4 5 r 

SLM  = SIJM*TX  (b) 

LOWf  bL'iJU 

A 

45/ 

IF  <IV.LT.«?0Ci)  GO  TO  7? 

LOHf  6?  9 0 

A 

4 60 

IF  (IV.  LT.no  00)  CO  TO  69 

LOW! 6300 

A 

4 59 

c**  ** 

***♦  UV  07 ONE 

LOHf  6310 

A 

4 b 0 

IF  <IV,Lr«  ?34DQ>  GO  TO  *54 

LOWf  63 i?u 

A 

461 

ir  <rV.GF.??500)  r>3  TO  r»5 

LUNt  bH  5 0 

A 

4 6 ? 

go  to  er 

L0WL634J 

A 

4 6 3 

54 

XK  * 30 0.  0 

LOWE  63 5C 

A 

464 

XI =<  V-l 3000 . 0)  /XX  ♦ 1 . 0 

LOWfbJbP 

A 

4 6 5 

L 1 M 

LOWf  63  7vJ 

A 

4 6 6 

Ll1*  53 

LOWF  6360 

A 

46  7 

GO  TO  T*b 

LOHF6390 

A 

4 b ft 

55 

X X * 5 0 Q . Q 

LOWr6400 

A 

469 

XI= (V-t 7500.0) /XXf57.0 

LOWf  6410 

A 

4 70 

11=57 

LOWt  64?Q 

A 

471 

L 7*  10? 

LO  WF 64  3 0 

A 

4 7? 

5 b 

no  57  MSL1.L? 

LOWfoHWU 

A 

4 73 

X0*XI-FLOAT<  N> 

LOWt  t>4'U 

A 

4 7 4 

IE  ixn)  59,55,57 

LO  WF  0460 

A 

4 7 4 

*>r 

CONTINUE 

LOWL  64  7 0 

A 

4 7b 

56 

T X ( 6)  =wm  *C6<  N> 

L O Ht  04  30 

A 

47  7 

GO  TO  b0 

LOWf  64  9 0 

A 

4 7 0 

59 

Tx((i»>r«<n»  *xa*«C8(N)  -ca(N-in 

LOWf  65  0 0 

A 

479 

57 


THIS  PACE  IS  BEST  QUALITY  FRACXlCABfJS 
COll  KUiuti  XU  UUC  — 


TaLlo  Al,  lasting  of  l-'urt i-.tii  loih'  l.OWTIiAM  -1  u'onl) 


T XI 8) = W (6) *T  X!8) 

L O Hi  65  1U 

A 

4b  0 

60 

SUM^SUMfT*  <8  ) 

LOHt  652  0 

A 

4 61 

IFUV.GT.  14500)  GO  10  H7 

LOWE6530 

A 

48  2* 

c; o to  f-q 

L O HE  654  0 

A 

4h  3 

WATER  VAPOR  COM  f NlIlM  10  MICRON  REGION 

10 HF  6550 

* A 

484 

61 

IF(  I V.C T,  1350)  GO  TO  62 

L0WE656G 

* A 

405 

TX<5)  - < i«.  1 6*557  8.0*  EXP  (-7.6  7E-3*  V)  )*H  <S) 

LOWF65/0 

•A 

4 fl  6 

GO  TO  6b 

LOHFf.ltlO 

•A 

487 

62 

I M I V «l  T*  ',*3S  (| ) GO  TO  bP 

LOWE  655  0 

•A 

488 

HATER  VAPOR. CONTINUUH  4 MICRON  REGION 

LOWE  6600 

* A 

4b  S 

X I = (V*?  350 • 0 ) /50 • fl*  1. 0 

LO  WE  66  1 0 

* A 

4 SO 

N H a XT  *1*001 

LOWE 66 20 

* A 

4 SI 

XM=  Xl-f  LGA  T ( NH) 

10  hf  tibJO 

•A 

4S2 

T X 1 5)  =T  5<  NHI 

LOHF  6640 

* A 

4 S3 

64 

TX(5)®TX(5|tXH* ( C 5 l NH )-C5fNH-l)  ) 

L O HE  6650 

* a 

496 

66 

T X<  5) =T  X<5»*  H(10  ) 

L OHtbUbt 

♦ A 

VJ7 

66 

SIN  = ELM  ♦ T X < *5 1 

LOW  i t'7U 

•A 

4SO 

IF  ( IV.  IF.  l.Vi  0 ,(1K  .TV  .C.T  . ?r  <>0)  Oil  10  7? 

to  HE  66  80 

* A 

4SS 

c • * * * 

****  NITROGEN  CONTINUUM 

LO  Hi  66  SO 

A 

5 6 0 

66 

IF  CIV* lT„ 2080)  GO  TO  72 

L O HI  67  00 

A 

5 0 1 

K4= 1-34 b 

LOHt  6 710 

A 

6 02 

TX  < 4) =C  4<  K4)  * W ( 4 ) 

LOWE  6 7 2U 

A 

5 e s 

SUH*  SIJP  + T X<4  ) 

LOME  6 7 30 

A 

‘>0  4 

GO  TO  72 

L O HI  6 74  0 

A 

*>C  4 

f »*  ** 

HATER  VAPOUR 

IOWF6760 

A 

r-06 

69 

IF  (IVtLT.l,’«00,4NI),IV.r.C,qiU>-)  Oil  Til  70 

LO  Hi  6 7 60 

A 

5 07 

if  <iv,u  .iu5?n.ANn.iv.Gr,i.Ti.noi  r.n  to  tj 

LOW  ('770 

A 

r>o  a 

GO  TO  7n 

LOHT  6 7 00 

A 

4 os 

70 

l-  I - 1 3‘. 

l.(iHFh7S0 

A 

‘>10 

GO  TO  7 2 

LOHt  68  OU 

A 

5 1 1 

7i 

I = I “ 2 59 

LO  HE  6 8 10 

A 

■>  1 r 

72 

K 1=  1 

LOHI  68  20 

A 

M 3 

IR  (HU),IT.  1,01  -20)  GO  TO  7b 

LO  Hf  ht>30 

A 

4 1 4 

H51  =ALOlU0  <W<  llMCl(I) 

LOHt  68  MU 

A 

4 1 6 

IF  (WSI  ,LT  J4MTI  TXH1  sl.-.0«77tl/*EXP(V.a‘jF,,j'T‘,.*HSl) 

L0HL68‘*l 

NE  H 

IF  iWSt .1 T*-?.346P>  GO  TO  76 

LOHI  6 860 

A 

‘•16 

TF  ( WSl  .GT  , 3.  468?)  GO  TO  75 

LOHFt-e70 

A 

M 7 

IE  < H S 1 » G T . 2 . 0 ) Kir.  4 0 

LOKF  bHH  0 

A 

'.'18 

00  7 3 K = K1 »6  7 

LOW( 68 Su 

A 

5 1 S 

IF  (WS1 ,UE ,EH (K) ) GOTO  74 

LOWE  b^00 

A 

52  0 

73 

CONTINUE 

LOWE  6 S 1 0 

A 

5 2 1 

7 4 

TXU)=TR<K)^nR(K-l)-TFMK)  )*(EH(KI-HSl)/|EH(K)-EWlK-in 

L ONI  bV20 

A 

5 2 2 

GO  TO  7 G 

LOHt  6S  3 v 

A 

5 2 3 

75 

T X< 1>  sG  ,0 

LO  WE  6 v»40 

A 

5 2 4 

76 

CONI  1NUF 

LOHI  ti<3r>0 

A 

5 "■*  >i 

Q • ••* 

••**  UNIFORMLY  MlXfO  GASES 

i C Hf  69  6 0 

A 

426 

IF  UV,LT.flOt»n.ANO,IV,Gr.SPO)  GO  TO  77 

L 0 Wfc  6 9 7 u 

A 

4 2 7 

I f ( IV.IT.IUS0.  ANP.  IV,  C.T,  12570)  GO  TO  7a 

LOWC6S81; 

A 

5 2 8 

GO  TO  03 

1 O H( 69  SO 

A 

5 2 S 

7? 

J = T-3  0 

LOHI  7 u 0 U 

A 

6 JO 

fi  it 


i 

I 

i 


THIS  FA&£  I S BEST  QUALITY  rpACTlCJUJItf1' 
Frtuil  liul  i 1’JluN  1S1L&D  TO  LlDC  


Tubto  Al,  I,lBtii»n  or  Fortran  L'oile  1 .OWTHAN  4 (Font) 


go  TO  7i 

LOwl  7010 

A 

4 .51 

7ft 

J=(I  V-i  ?950>  /5*l{ilb 

10  ME  70 20 

A 

63  2 

79 

IF<H(2I  .LT.1.0f-20)  GO  TO  aS 

LOME  7 0 3 0 

A 

r>  3 3* 

K 1-  1 

L0HF  70  4U 

A 

63  4 

HS?=ALOGlO  CM  (?>>  *C?  (J) 

LOHt 70  5Q 

A 

536 

TF  THS2.1T  .-2. J9b«)  T *1  2>  = 1 . - . 0 S7  M 7 *f  XP  ( 1 . g 5);99  n<H  s 2) 

UIHF7D60 

NC  N 

IT  THS2.LT. -2.  J9b#>  CO  TO  8? 

LOHL  7070 

A 

5 3b 

IF  <W5?.GT.J .58821  GO  TO  8? 

LOHl 7080 

A 

63  7 

IF  <HS?.GT,2.0)  Kl=  40 

LOHt  70  90 

A 

G 3 8 

00  80  K«<1,G7 

LOME  7100 

A 

539 

IF  (HS2.LE.F  W<  K)  > GO  TO 

LOHf 7110 

A 

6 4 0 

AO 

continue 

LOHt  7120 

A 

54  1 

81 

TX(?)=TR(K)t  (1  R<K-11-TR(  K)  ) *(  FH(K>  -HS2)  / ( F H(  K > -F  HI  X-  1)  ) 

LOHt 7130 

A 

64? 

GO  10  83 

LOHi 7140 

A 

543 

8? 

Ta(?» *0.0 

LOHF7) GO 

A 

64  4 

83 

CONTINUE 

LOW*  71  o 0 

A 

54  5 

C** 

******  O70NF 

towr  7170 

A 

54  6 

IF  <IV.LT,576,0R,  IV*GT. 327(1)  GO  To  ft* 

Lonr  7 i a o 

A 

54  bf 

L=T-46 

LONE  7190 

A 

64  7 

<1*1 

LOHf 7200 

A 

54  f 

IF  < W ( 31 ,LT. t . 0 E- 20 ) GO  TO  87 

LOWE  7? ! 0 

A 

9 99 

HS3**LO(;iO  < WC3)  > *C3  CL  ) 

LflHE  7220 

A 

55  C 

TF  <WS3. IT. -1.6778)  T X ( 3 > = 1 . - . 0 59 19  4 ( ? * 3b  ?86  3*  WS  3) 

LOHT  7230 

NL  M 

IF  TWSJ.LT.-  1.67781  GO  TO  87 

LOHt  7240 

A 

561 

IF  (HS7.GT  .3  • 9 3 4 *»  1 GO  TO  Rfa 

LOMF7260 

A 

66? 

IF  (HSJ.F.T  . 1 .5)  *1=36 

LOHI 72c>0 

A 

56  3 

nn  (14  K=Kl,b7 

LOHt  7270 

A 

56  4 

IF  (HS3.IF.F01K))  GO  TO  8S 

L O HE  7280 

A 

65  5 

84 

CONTINUE 

LOWf  7290 

A 

66  b- 

85 

TX( 3)  =TR(K) -1 TR(K) -TK (K«l) ) *(KO IKI -WS3I  /(F0(K 1-F0TK-1I 1 

lOHF 73 00 

A 

56  8 

GO  TO  87 

LOHt  7310 

A 

56  9 

86 

TX( 3) =0.0 

LOHE7320 

A 

GhO 

87 

CUN  UNUt 

LOHF  73  31' 

A 

6o  1 

C** 

******  AEROSOL  extinction 

LOHt  rj*.u 

A 

5 b ? 

ALAH*1.0E*4/V 

LOHE73G0 

A 

66  J A 

X x-.  0 . 0 

LOWF73G0 

A 

bbJll 

yy=o. n 

LOHL  7370 

A 

66  30 

c** 

******  Tt’MPOPKRY  FOG  COPPFCTION  FOR  VIS  OF  L OH  2 Km. 

LOHE7JOO 

A 

6 b 3 U 

c 

IFIVIS.GT. 0.0. KNO.VIS.LT. 3.0)  XX  = 0.1911 

LOHt  7390 

Nc  H 

c 

TEMPORARY  FOG  SUPPRfSSFfl 

LOHT  74  0 3 

Nt  H 

IF  (I  HA?E  . FO  . D.OP  .XX.GT  . 0 . 0)  GO  TO  90 

L OH  £7  9 1 0 

A 

4 64* 

CO  88  M=l,44 

LOW  792u 

A 

6 6 5* 

XO=  ALAH-VX  CN 1 

LOHT  79  3 0 

A 

5 0 b* 

IFCXD)«9, BB, 88 

I.OHf  71.90 

A 

5 67* 

88 

CON  T T NUF 

L OHf  7950 

A 

4 t>  8 A 

fl9 

XX  = (C7TM -C7  tN-1) ) *XD/< VXTN) -VX(N-l) ) «C7TN> 

LOHt  79b0 

A 

4 688 

YY*(C7MHt-C7ATN-l  ) )*  XCVIVX  TNI  -VXTN-1)>»C74  Ht) 

LOHT  79  70 

A 

5btJC 

90 

TX(10)=YY*W<  7) 

LOHE7480 

A 

6bOO 

T XT  7>  *X  X*W(Y) 

LOHt  7990 

A 

6G9* 

SU^=SIINTTX  (7  ) 

LOHI  7r>  0 0 

A 

5 70 

i 

J 


jtfML  OOi'*  KU**^1**0  1 


Table  Al,  l.isttiij;  of  I'Virtran  C'odo  1 ,0\VTI{  \N  -1  U'onl) 


T XfO)«FUM 

LOW*:  7510 

A 

4 7 1 

DO  94  K*4,KMAX 

LOHt  75?tJ 

Nt  W 

1 F ( T X ( K i .n  .0.0)  GO  To  9? 

LOWF  79 JU 

A 

67  3 

IF  •LE.0.1)  GO  TO  *U 

LOWF  7940 

A 

5 7 4 

IF  <T  X (K)  , GT  . ?0.  ) GO  TO  91 

L 0 W(  7950 

A 

6 7r' 

T X I K 5 =T  XP  C-TXtK) ) 

LOHE  7560 

A 

6 76 

GC  TO  94 

10U757U 

A 

677 

91 

T X CK) =1  »0-TX  (KM  0.5*TX(K)  MX(K> 

LOWt  7 5 80 

A 

578 

60  TO  04 

LOWf  7590 

A 

5 79 

9r 

T X <K)  =1  .0 

lOWfc  7b 0u 

A 

68  0 

r.  c To  Q<* 

LOHT  7b  10 

A 

58  1 

9.1 

tx  m =o . 

LOWl  7b_»0 

A 

58  C 

94 

CONTI NUF 

LOWl  7630 

A 

4 8 3 

T X f 9)  *TX< 1)*TX < ?) *1  X ( 3)  * TX  ( 9 > 

LOSt  7640 

A 

564 

if  tiV'Gt  .ucoai  TXU)*TX(M 

LOS{  765  Q 

A 

68  5 

IFUCilSS.fO.O)  GO  To  1210 

LOSr  7660 

N t W 

AlflMs  1.  0E«-04/ V 

LOWf  76  70 

Nl  W 

MfllK  =FF  ( T Q HX  ( IK)  ,V) 

IOM  7b9  0 

Nf  W 

1LNFN  »CTX(9)*TX(10H/(TX(7MTX<6)| 

t out  7690 

N ( N 

TSNTW  s<TT<7)»rx(b>l/rxilO> 

LOWt  77 00 

NEW 

OTAUsTLOLO-  TI.NF  W 

{ OWt 7710 

Nt  W 

IF  < OTAU.  LT.  1 ,CIE -9  .ANCl.TLNtW  .1 1.1.  OF -5)  GO  TO  1104 

LllWt  77C0 

Nl  W 

SUMV=SUMV»  0.  6*  BPIK*  pT  AU*  (TSPLP+TSNE  W> 

L O H f 7 7 3 0 

Nf  W 

TL  01  0*  Tl  Nf  H J TS0L0-T5NFW 

LOW! 7740 

N F w 

10r’0 

CONT iNUe 

LOWt  7 76  0 

Nl  H 

nn<. 

CONTINUF 

LOHt  7760 

NEW 

TAUC-0 

LOW*  77  7 0 

NF  W 

IF (NMlN.LE. O.O.ANO.Il.EO.l)  TAUG=TX (9) 

L OWL  7 7 90 

N(  H 

T 1*T (M*1 ) 

LOWf  7 7 90 

Nt  W 

I F CTROONn.r.T  .0.0)  Tl  = TROUNP 

L OWC  78  0 0 

NEW 

PBGeFF 1 Tl , V) * T AUG 

LOSE  7b  1 0 

NFW 

IF  (hhin.lf.  0>  suhv  *sumv  ainr. 

l.OWf  7820 

Nt  H 

SUMVV=SUMV 

LOWF  78  30 

Nf  H 

n UV.51.1V1)  HL.10K  = 1.1) 

LuWt.  7 84g 

Nl  N 

IF  av.CF.lV2)  FACTORS. 9 

LOWt  7850 

Nt  W 

SCMV*  <1 .0E*0  4/V**?> *SUM  V 

LOWf  7860 

NFW 

R Ap SUM* RAnSUN»nV* FACTOR* SUHV 

LOWL  7070 

Nf  H 

IF(JP.FC.O)  PRINT  llbO,  V, ALAMjSOMV > SUNVV » R APSUM, T < (9) 

LOHt  70  8 0 

NL  N 

IFISUHV .GC .RAOHAX)  VRMAX*V 

LOWt  7 6 % 

Nf.  W 

IF  ISUMV.GE.HADHAX)  RAriNAXsSlINV 

LOHL790C 

Nf  H 

IF! SUMV.IF.RAOMIN)  VRM1H=V 

LCWE79 1 0 

Nf  W 

IFISUNV.LE  .R  AOH  IN)  RA  DM  [ N|-  SUHV 

LOWl  79 CO 

Nt  H 

WRITE! 7,3010)  V.SUHV , SUWV V ,RA0SUN , TX ( 9) ,TX ( 1 ) 

LOWL  793  U 

Nt  W 

^ n l n 

FORMAT  < FlO.t , 1 IF  10. ?> 

LOHt  7940 

Nl  W 

icip 

TX<  10M1.-TX<10> 

LO  W f 79  5 0 

Nl  W 

AO* 1 , -TX (9) 

LOWE  796  0 

A 

6b*'U 

IF  < tV.F  Q. 1V1  .OR. IV.FO.I VC)  AB=0.4*AR 

towr  7970 

A 

6 86C 

SUMA-SUHAtA9*nv 

LOWF  7980 

A 

56  61) 

IFUEWISS.fQ.l)  GO  TO  1C?0 

LOWf  79  90 

NF  H 

IF(JP.tQ.O)  HRI  TF <b  *423)  IV, ALAM,T X ( 9) , tTX ( K) ,K= 1 , 7 ) , TX(  1 0) , 

, SUMA  LOHt  8000 

A 

5 8 7* 

tiO 


THIS  PAftEaS.BSST  QUALITY  FRA0X1CABU 
TOkJM  tt>F  i‘ JUStU,Sjl  IED  TO  DD  C 


Tallin  A l,  I,i^tin(T  of  Fortran  Code  l.OWTHAN  4 (Ormt) 


T CONTINUE 

1*  <1 V.  GL.  IV  2)  GO  TO  95 

r.o  To  r-o 

«EftO  «*0  0 t T<V 

JF  ( IE  HISS,  £0.1  > PRINT  1175,  VRMIN,R40MTN,\|R4«X,  BA  DM  AX 
6 FORM  AT < * RAOMIN  * ,F 1 ? , 3 , E i 2 . 5 * / » ♦ RAOMA*  *»  F 12  * 3, FI? *5) 

ENOFUf  7 

JSTOR-0  1 IFIN0*Q 
ARsl. 0-5UMA7 FLOAT < I V2-IV1I 
PRINT  424  , T VI  <IV?vSUHA,AB 
TP (IF RISS.FG.1)  P^INT  443.RADSUM 

F OrhaT  (*  INT  r GRATE  0 RAOIANCF  *•  , F 1?  .*3  ,*  HAT  T CH  -2  SR*) 

PRINT  400,1** 

iniXY'CQ,0)  GO  TO  100 

GOTO  *%,?»S7,98t 100 1, IKY 

RE  AO  406,  VI  ,v?,nv 

AVWrJ  00  00.  /VI 

At  AMs  10  000  ,/ VJ 

PRINT  418,  V X ,V2  *OV  ,IUAH  ,AVH 

SLM  A*  Q . 0 

GO  TO  «*T 

if (Hoaru, eo. 0)  go  to  200 
0 c TO  300 

RFAO  4C0,4OOFL,IHA7f>  ITYPl',LfN,  JP,  IH,H1,M2,H3,HL,  I EHISS,  RO  ,T  HOUND 
IPUEHISS.FO.II  PRINT  1170 
lF<IF*i«;s,Fil.OI  PRINT  1171 

lEnstop=l*n 

PRINT  40n,HOlFe.,  IHA7F,  ITT  PC  , LF  N , JP  , IM , Ml , MI*  ,M  , HL  , I EM  I SS,  RO  , TBOUN  P 
GO  TO  200 
^ TC F 

FORMAT (111 3, ?F10. 3) 

D fORMATlM  PROGRAM  WILL  H?  EXECUTED  IN  TMF  EMISSION  MODE*) 
l FORMAT (*1  PROGRAM  HILL  3d  F XtCUTEO  IN  1 HF  TRANSMISSION  MOOE*) 

FORMAT  (0619.31 

FORMAT  (F6.1  iIMFIC,  3,Ff»,  t .7F10.3H 
FORMAT  <4<F6.  3,  2F7.  4)  ) 

FORMAT  <1SF5.?> 

FORMAT  (6FP.?> 

FORMAT  ( 7 F 1 0 * 3) 

FORMAT  MORIi'ONTAL  PATH,  ALTITUDE  *,F7,M1H  KM, RANGE  = 

IF7,  J,  3H  KH) 

FORMAT  <//l3X,5C*H  SLANT  PATH  BET  HF  EN  ALTITUORS  HI  AND  H2  WHERE  Hi 
lr,F7,3,AH  K*1  H?  *,F7«3,16H  KM,?ENITH  ANGLE  =,F7.3,RH  OEGRFES) 
FORMAT  <//lOX,39H  S l ANT  PATH  TO  SPiCE  FROM  At 1 1 TUDE  HI  =,F7.3,19H 
1 < H , 7 CN I T H ANGLE  *,F7.3,8H  DEGREES) 

FORMAT  (/20X,ltJH  MODEL  ATMOSPHERE  ,T1,11H  » TROPICAL) 

FORMAT  (/?0XV1«H  MODEL  ATMOSPHERE  = MIDLATIYUOf  SlJHMFR) 


FORMAT  </20*,19H  MODEL  AYMOSPHFRf  ,I1??1H 
FORMAT  T/20*,18H  model  ATMOSPHERE  ,I1,?1H 
FORMAT  (/  2 OX , IBM  MOBFL  ATMQSPHF  RE  ,U,21H 


MTOt AT IYUOE  WINTER) 
SUB-ARCTIC  SUMMER  ) 
1 S62  US  STANDARD  ) 


L O HE.  «0  1 0 
LOWE8Q  20 
LOHE80  30 
LOME  8040 
LO WC  80  50 
LOHt  ^0  60 
LOWF90  70 
L OWt  00  8 0 
LOWE8090 
L OWL  8 1 0 0 
t OWL  8110 
LOHt  0120 
LOHt  0 1.30 
LOWE  8140 
LOWF  8150 
LOWE  8160 
LOWf 0170 
LOWf  0 1 80 
t OWE  81  90 
LOWf  fi?00 
L 0 WF  8 ? 1 0 
LDWF6220 
LOWE  0230 
LOWE  82  4 0 
LOWE  0260 
LOMF8260 
LQWL  8270 
LOH820Q 
LOHF029O 
LOWF  8300 
LOWF 83 10 
LOWF8320 

LO WC033D 
LOHt 8 34  0 
LOME  63  50 
LOWE036u 
LOWE6370 
LOME  83  80 
LO  HF  8390 
, LOWf ft^OO 
L O HE  fi  4 1 0 
LOME  0 4 20 
LOME  843G 
LOMf  0440 
LOH  E84‘i0 

LOWf  84 bu 

LOWE  84  70 

LOHt  04  8u 
LO  H f 3 49  0 
LOHF05  00 


NF  W 
A 588 

A 589 
A 690 
NEW 
NEW 
NEW 
NF  M 
NEW 
A 6 913 
NEW 

NE  W 
A 59  1C 
A 5910 
A 5 9 1 1 
A 592 
A 593 
A 594 
A 595 
A 69  fe* 
NEW 
A 098A 
A 5983 
NEW 
NEW 
NEW 
NE  W 
Nt  H 
A *9AF 
A 69  9* 
NF  H 
NEW 
NEW 
A no  j 
A 602 
A 60  3 
A 604 
A 606 
A 60  6 
A 607 
A 60  8 
A 6 U 9 
A blU 
A 611 
A 612 
A 613 
A 614 
A 616 
A 616 
A 617 
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Table  A1.  lasting  of  I’nrtruu  C'oiU'  UHVT1\AN  4 (rontl 


hI  5 

FORMAT 

416 

FORMAT 

417 

FORMAT 

416 

FORMA  T 

1 4 H CH- 

419 

F OR  MAT 

LOWt.  BboL 
' * L C HT  66  70 
LOWt  6b  Bo 
LOHf  6 b9Q 
LOWt  67  00 
Lnw{  67  10 
RANG  l OWE  6 7 20 
LOWf  ft  7 3U 
LO  Wf  6 7^4  C 


0X,lftH  HO  CL  L OTHOSuHfRt  fIl,,?lH  = SU0-A*CTIC  WlNTF  R ) LOHUIblQ 
OX  1 1 6 H HA?(  NOOfL  , 1 1 , 3 H = #A5,13H  VISUAL  RANGE)  LO  HL  6b2t 

5X*HA7F  MOOR  KH  VISUAL  RANGE  AT  SEA  LFVFL*)  IOWL6S30 

0X,?1H  FRFQllEtsi  f-Y  RANGE  Vl  = ,F7.1I13H  CH-i  TO  V.?=  i F7  ■ 1 , 1 LOHf 6b 40 
R OV  =,r6.l,°M  CM- 1 (,FG.?,*  - MICRONS  »*>  LOWlrt‘«bO 

0X,3flK  FOUIVAI.CNT  SFA  LEVEL  ARSORP6R  AMUUNTS//:i  XI  lOHWATLOWt  6560 
JEW  VAPOUR  CO?  t TC  # O70NE  N I T P OG  E h (CONT)  H7P  (CONT)  LOH  1&T* 

? HOI  SCAT  AEROSOL  070NL  (U- V ) / ?4  X , 7 MGM  CM-?t  1 IX,  2HKM,  u OH|  6‘joO 

30X,6HATM  CM  , LOX  , ?HKM,9X  , 7HGM  C M*  ?,  1 0 X , ? HKM  , 1 J X,  ,?  H K M , 1 0 X,  • .H  A 1 M CH|  LOWfBr)90 
420  FORMAT  (LHlt ///10X,  » VERTICAL  PROFILES  * , 64X  ,*PS  I*  ,6X,*PH1*  ,6X  ,*  LOWt  6b0(] 
1 RFT  a*  ,4X  ,*THfT  A RAKC>K*I  LOWtOblU 

4 ? 1 FORMAT  ( 7 1 0 X , 6H  U(1-B)*»(E14. 3) t 74X ,t 1 4* 3 , ? ft* ,L 1 4 . J / ) LOW  F 6 bJQ 

42?  FOPHAT  (IH1,/10X  f 3?H  havFLENGTH  TOTAL  H?0  , 5X4MC02  * , *i  X , bU'HF  Bb  3 j 

14H070NE  N?  CONT  H ?0  CoNT  HOL  SCAT  AFRO^OL  AEROSOL  I NTt  GR AT  f ULOwt  fib*.  0 
? / 1 1 X , 1 4H  CN-1  MICRONS,  T (4X^HTRAMS)  t4X  ,?0H  APS  AfiSORPTlON  J L O ML  665J 

U2\  FORMAT  ( 13X*IO»lLF9* 4 ) 

4?4  FORMAT  l*  ISTFGKATrn  ABSORPTION  FR0M*;l5,*  TO*, IS.*  CM-i  = * * F 1 0 
1*,  AVERAGE  TRANSMI  IT  AnCI  **,F6.4) 

4?5  FORMAT  (lux,  7F 1 0 * 3 ) 

426  FORMAT  </?PX , *AtftO>OL  SCATTERING  NOT  COMPUTE O , IH A2 F ~ Q* ) 

427  rpRMAT  «lHl,//710Xt TOM  HORIZONTAL  PROFILES/) 

438  FORMAT  (10X,*  H 1 = » , F 7 *3  t * K*  , M?=*  , F 7 , 3 , * KM . ANGLE8  ¥ • r t>  * 4 , * G F OH 
IE  = *,  F7,?,*KH, BETAs*f FA , 9,*,  vl S=*,Fb* 1) 

4?9  FORMATT  3F10.  3»FFS,1  ,?F10*  T,?F10,3) 

430  FORMA  T{  10Xf*  INPUT  M FT  TO  P OLOC I C Al  0 AT  A I • / i 0 X , * ?r.*  , F 7 , ? , » KM,  P=  ♦ , F7  L OH  F 8 7 l">  U 

I*?,*  Mp,T=**F%l  ,*  C,  De  w PT  . TF  HP*  , Fb.l  ,*  C,  RfL  HUM  I OI  T V=  * , F*j  , i , LOV.F6760 
7*  H3D  PFNSITYs^,  1PF9. 3,  • CM  H-T^/IQV,*  07OKt  ill-  NS  IT  Y=  * , M .?  , * ulCHE«77Q 
3 M - 3 , VISUAL  RANGF= * .OPFfr * lf*  KH, R ANC F= * , F 1 0 . J, ♦ KM  * ) lOXEoZBC 

431  FORMAUMFb.  ?,3F7.^))  LeHM?V0 

43?  FORMAT  <*  STARTING  FARArtF-TfRc  HI  ANU  AN  CL  F HAVE  Ucf  N R FJ  C F INF  0 1 H 1 -LU  Vi  880  i, 

1 * » T 1 0 • 1 , ^ANGL  c «%n(I.E)  LOWE  8610 

433  FORMAT  (*  TRAJECTORY  Ml^SFS  EARTHS  ATMOSPhFRF.  ClOSFST  TISTANCt.  OFLflWI  6630 
1 APTWOACH  IS*,F10  .3  il X, / , IT, »FNO  OF  CALCULATION*) 

47 4 FOPHAT  (lRXv  I4»F6«1 10.3)  ) 

435  FORMAT  M5,F7.1t6ri0.3,4Fg,<,,F7.l) 

4?G  FORMAT  (*  HMIN  = *iFin*3> 

43  7 FCRMA1  (*  PATH  InTERSFCTS  fcA^TH  . path  CMANGEO  TO  TYPE 
1 O.H  KM*} 

43ft  FORMAT  <*  CMOICr  OF  T HO  ° H HS  ^0  R THIS  CASF  -SHORT  F ST  PATH  T A KF  K . LCHLoo'lO 

1 FOR  LONG  F R PATH  SET  LtN=l.*>  L O HE  6 H D 0 

439  FORMAT  <»  CHOICF  OF  T Hpi  PATHS  FOR  THIS  CASE  -LONGFST  PAlM  TAktN 

l FOR  SHORT  °ATH  SET  L?N  - C **) 

440  FORMAT  (*  H2  HAS  SFT  LT  SS  THAN  HMIN  AND  HAS  BEEN  KFSf.T  f DUAL  To 
1 H HI N 1«  E • H?  5 *,F10,3> 

441  FORMATS*  HOOfl  ATMOSPhERF  MO.  7*,/  4X,*7  <KN)*,3X,*P  (MJ)*,4X, 

1 * T |C)  new  FT  XRH  H?0<OH,M»3)  G3(GM.M-3)  NO*  OCN.*) 

44?  FORMAT!*  FOG  CONDITION?  HAY  EXIST  AT  SF-A  L T V£L  F UK  THIS  VISUAL  RALOHt  6li70 
INGE**/,*  IF  SO  THFN  ASSUME  THE  TRANSMITTANCE  HUE  TO  FOG  IS  GIVFN  lOWF6k>&U 
26Y  THF  TRANSMITTANCE  A|  0,55  MICRONS*)  LClHEd990 

1109  FORMAT  (013)  LOWt  9u00 


LORE 863U 

L n uf pH40 

LOWt  rt.‘)5u 
LOWE  i>  0 
WITH  H ? =L O Wt  6676 
LO  Wt  60  6 C 


LO  ML  841  0 
LOWF6'l?o 
LOWt 69 3: 
LOWt 6940 
LOWf  6 Obi, 
l OWE  69 bU 


A b 1 6 
A b i 9 
4 h?0 
A 1%  ? 1 

a t;: 

A bJJ 
a b r 4 
A b *. 

A GJO 
A b?7 

m b ? 6 
Nt  W 
A b 3 0 * 
A bJl* 
A GJJ* 
Nt  W 
A U3  4A 
A b 3 4 1 
A 

A t-  3 
A b3  7 
A b,*fl 
a h.tq 
A b4  C * 
A b 4 1 • 
A hUO* 
A b4  3* 
A b**4* 
A b4f«* 
A t»4li 
A b4? 

A b4fl 
A b49 

A bb  n 
A b‘  1 
A b‘>? 

A bb.3 
A br4 
A b'5 
A Obb 
A b r>7 

A br>  6 
A o‘*M 
A lib  0 
A bbl* 
A bb?* 
A b h X * 
A t,  b 4* 
A bbb* 
N‘  W 
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LOHl 90 10  Nt  W 

IQRF902Q  Nf  W 

PER  MCRON*,*  L OHCSu  30  NtM 

LOHt  90  40  NEW 

10  Wf  90 9 0 Nt  W 

LOHE9060  A 6t>6* 


(5  3 


Ur>F  FORMAT  (11RX,  •NITRIC  AGIO*) 

11^6  FORMAT  <1H1, 30X , “RAni ANCE f H AT T S/HM? - ST E R-XXX ) * I 
119  7 FORMAT  l.TOX  ,*TR  (CM- 1)  WYt.<MICRON)  PF&  QM-1 

1 INTEGRAL  TRANS*) 

1160  FORMAT  < 30X,F8.  1 1 P‘13 .6  »?F  1 T.  5 ,Fl3.o  ) 
f N 0 
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SUBROUTINE  P»TM<RUY,Uf'ATH1TIIlW>  Bt  _ 

C01M0N  7 < ^ ,P(f  ,34)  ,1  (7,  31.)  ,EH  t 11,  JR)  . NH<  7, 3m  , M>NL  ,RE , c 
nmrNMuN  kiArusin  t'BBYim 

COMMON  /F  Ml/  KtiISS,KH»X,*N6tF«lEN,HHIN,lJ*J»»J7t  JKlN»Jr« 
COMMON  /FHT/  IL.IKHSX  .LENSTOR  , NIL ) X (1 t > » E < U > 

COMMON  /E  M3/  H1,M?,NP1  .MOREL  . 

IMITYPF  .EU.lt  r,0  TO  7000 
IF  ( I T Y Pi. . FQ  ,2. UNO.  HI  .f  Q «H?)  JS*J1 

IF  (H7.C-T.  HI.  »ND.  ANC.LF.C.T  . 9*.  . ANP.SP1  . EO.  t)  J1  = J1  -1 

IF  ( (ITT  FF!n  1 f)  .s  NO.  (HI  ,GT  . HJ  > . »NO.  (L  ENSFClR.  FO.  1)  ) JT-d!-  1 

IF(ITYFF.EO.Y)  J?  = Nl 

PRINT  HOT,  J1.J2 
1 1 OP  FORMEK  PI  I) 

pp  i mT  ^ 1 0 

9!0  r0RNHT  <//,2)X,*  CUHOUTI  it  AtlSORRFR  AMOUNT  S EON  THE 
1IC  PATH**// 

?tox,*Hco*  ,sx,*co?** » ♦<>*•. *x,  *n?*,  ax»  *h.,o  r*ff>y,*poi  s 

^ * Af  PC*  »S  X»  fO’'  UV*,7>t*H20  C*  t 7X  , *HNO  3*  >S*  * f T AVt4  ) 

no  iosr  TK=tffcfl 

TBPY«1<1  =0* 

(»0  IDS?  K=  1 ftTC 
HPATHdK,  K)  = 0. 
in  *;*»  CPnTINUF 
tt-  N=n 

NIL  =NL- 1 

IL=Jl*l 

IJ-II  *WLL 

!U'  lObfl  K = 1 i K M A X 

MKl  = 0, 

10H[t  CONTINUE 

IF  (ANGLE  *GT  *0  0,0)  GO  TO  iQbl 
U N - l • 

IL= Jl-1 
Hh I M-  i«  0E-6 
TJsNll 

10  61  CONTINUE 

no  lbSO  iKai.hd 

IF  <lEN.t0.IH  IL-IL-I 
TF  (LFfJ.EP.t)  I L = TL  +1 
U=  IJ-I 

IFUl.FO.C)  OO  TO  1050 
On  1055  K- 1. KMAX 
NIKlsEI*)  FML  AYUL.K) 

WPATH(IK.K)fW(K) 

IDEA  CONTINUE  „ ,„e, 

IF<IL. IE. 0. OR. IL.OE.NL)  GO  10  105J 

Tf)AR=<  T <M,  IL>  *T  ( M , U.  ♦ 1)  ) *0 .0  . . _ 

C JFXTRA  IS  1 ONLY  WHEN  PROGRAM  NEVER  LFAVES  ONE  EAY'R 

I F ( JE  XT  RA « EQ  • 1 ) t RA  R«  < T < M , Jll  *T  < M , J 1 *1 1)  »0 . 5 


PATH 

Id 

NL  H 

«A  TH 

:o 

NEW 

PATH 

3U 

Nf  N 

PMH 

40 

Nt  W 

path 

SO 

NEW 

PA  TM 

bO 

Nt'H 

PA  T H 

70 

Nf  H 

PATH 

00 

Nl  W 

PQ  T H 

9 0 

Nt  W 

PA  \ H 

10Q 

NfW 

PATH 

110 

NEW 

PA  I H 

iro 

Nt  H 

PATH 

130 

Nt  H 

PAT  H 

140 

Nt  H 

PAtw 

1‘iO 

Nf  H 

RPATh 

U)0 

Nt  W 

path 

1 7 C 

NE  W 

PATH 

160 

Nt.W 

PATH 

mo 

Nl  W 

P A In 

r o ti 

Ni  W 

° AT  H 

rto 

Nt  W 

PATH 

Nt  W 

PATH 

?J0 

N t W 

°aTh 

;'40 

NL  H 

P AT  H 

rr>o 

Nl  H 

PATH 

:so 

Nt  H 

path 

?7  G 

NEW 

PATH 

rae 

Nt  W 

PAT  M 

zno 

Ne  w 

r>  A T M 

3 00 

Nt  N 

PA  TH 

310 

NFH 

PATH 

3 00 

Nf  W 

PATH 

3 30 

Nt  N 

PA  1 H 

340 

N t N 

PATH 

J*>0 

Nt  w 

path 

3b  v 

Nt  N 

PATH 

3 70 

Nt  W 

path 

36  0 

N|  W 

pA  TH 

300 

Nt  U 

PA  TH 

430 

Nt  N 

PATH 

4 10 

Nl  W 

pA  TH 

4 < 0 

Nt  W 

pa  TH 

430 

Nt  W 

pa  ih 

440 

Nt  H 

PATH 

4 SO 

Nt  W 

pat  h 

460 

Nt  W 

pA  TH 

4 70 

NTH 

PATH 

460 

Nt  W 

PATH 

4 90 

Nf  W 

PATH 

i soo 

Nt  W 

64 


uyULkEifcii, 


1 D S3 

1103 


10  6;? 
10*0 

110* 

2 0 00 
20 


900 
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CONT INUF 

PATH 

510 

Nt  W 

T BBY  <IK>*  T6A* 

PA  TH 

$20 

Nf  H 

no  1103  K*  1 » KMA X 

PATH 

$30 

KFH 

F ( \ ) - M ( K) 

PATH 

$40 

Nf  H 

CONTINUE 

PATH 

$50 

NL  H 

IF  (ANGLE  *LF  ,90 . 0 . AND «IL  ,GQ*NLl)  00 

TO  1 10  >» 

PA  TH 

$b0 

Nt‘  W 

IF  (ITY  FE.Fl^.ANO,  ANGlF  .LE.90,0>  GO  TO  10b? 

PATH 

570 

NlW 

IF  (ITYPE.n.’.ANO.  Lf  K.Ea.l.ANO.R. 

EH.J?)  GO  TO 

1104  PATH 

$80 

NC  W 

IF  (ITYPF.EQ. J. A NO. LENSTOP. Ta. 0, AMD 

• U,EQ,J?)  GO 

TO  ltQ4  PATH 

590 

NF  H 

IF  < IL  . EC.  JMIN.  4NO.  HMIN.  CT  . 01  LFN^l 

P4  TH 

600 

Nf  W 

IF  ( II , EQ, 1. ANn.HMIM.EE. 0 .fl  GO  TO 

1104 

PATH 

610 

NEW 

IF  CLEN.EO.3C  GO  TO  1063 

PATH 

630 

Nf  W 

If  (IL.FQ.JMIN.  ANn.IJ.EO.  IL  + MLL)  IL 

= IL-  1 

PATH 

b 3 0 

NL  W 

IF  (T  TYPE . EQ.3.AN0 • IL .rQ. J2)  GOTO 

1 104 

PA  TH 

6 40 

N F W 

ECN  I IM'F 

PAIN 

650 

Nr  w 

PRlNTSnO, IK,  (WPATH(IK,K)  ,K=l»8t  , WPAT H ( I K, l 0 ) , WP A T H < I K , 1 1 ) , TORY ( I K> PATH 

660 

Nt  W 

continue 

PATH 

6 70 

Nt  W 

I KM  AX  = f 8 

PATH 

680 

Nt  U 

L IN-LFN^TO^ 

PA  TH 

690 

NEW 

Kf TURN 

PATH 

7 j Q 

NEW 

CCNT  INl’r 

PATH 

710 

NT  N 

I PIN  190  0,1*,  (WPATHIIK,K)  , K= 1 , 8 ), W P AT H ( I K , 1 0 ) , HP  A TH ( I K , 1 11 » T R R Y ( I K ) ° A T H 

720 

NEW 

IKf*AX  = I K 

PA  TH 

730 

NfW 

lcn=len*tor 

PATH 

74  U 

NL  W 

return 

PATH 

75  0 

Nt  K 

no  uOr»r  k-i,k-»ax 

PATH 

7b0 

N t W 

wpathii,k)-w<k> 

PA  T •« 

7 7 C 

N t W 

C ONT  I Nl*c 

PATH 

7 80 

Ni  K 

I F (MOr^L . F 9 , 0 ) JIM 

PATH 

790 

Nf  W 

J 2 - J 1 

PATH 

600 

NFW 

T H R Y ( 1) =1  ( M, ji> 

PATH 

810 

NF  W 

TK  HAX  = t 

PATH 

82  0 

Nt  W 

PRINT  11D9,  jl.Jc 

PATH 

0 3 C 

Nt  W 

PRINT  qic 

PA 

NFW 

IK-  1 

PA  I H 

0\ul 

Nt  W 

PKINTRO0, IK,  l WPA  TN ( IK ,K)  tK=l  ,0)  ,WP4TH(IK,10).HPATH(IK>  11) • T 9B  Y ( I K)  P A T M 

06C 

Nt  W 

HMJNr i ,QE-h 

PATH 

870 

Nt  W 

RETURN 

PATH 

0 BP 

Nt  W 

FORMAT (I*, 10  F 11 . 3, F 1 0 * 3) 

PATH 

090 

Nt  U 

End 

PA  TH 

S00 

Nf  W 

THIS  PAGE  IS  BEST  Q.UALITY  FRACTICABJUS 
15804*  COt  Y J-  LUuNti>tUuU  TO  DD,C  


Table  Al,  Listing  of  Fortran  Code  LOWTHAN  4 
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SUBROUTINE  HN03  <Y,HAHS> 

HN  03 

10 

NEW 

01  HEN  SI  ON  Hi  ( 151  ? H2I16I,  H3(1T) 

HNC3 

20 

NEW 

C ARRAY  Ml  CONTAINS  HN03  A9S,  CO EF  <CH  - 1A TM -1>  FROH  ftSQ  TO  920  CH-1 

HN03 

30 

NEW 

DATA  HI  72,19 7, 3. 911, 6. IS  4, 8. 150, 9. 217, 9. 46 1,11. 56, 11. 10,  11.17,12 

•4HN03 

40 

NEW 

10, 10. 69, 7. 50 9, 6. 136, 6. 899, 2.0667 

HN03 

50 

NEW 

C ARPAY  HZ  CONTAINS  HN03  ABS,  CoEF CCN- 1 AT M- 1 ) FROM  1275  TO1350  CM“1 

HN  03 

60 

NEW 

0 AT  A H2 72, 828, 4. 611, 6, 755, 0. 759, 10. 51, 13. 74, iB.Oq, 21. 51,  23. 09, 21 

.6HN03 

70 

NEW 

18, 21. 32, 16. *2,16.42, 17. 07, 14. 86, 8. 716/ 

HN03 

80 

NEW 

C ARRAY  H 3 CONTAINS  HN03  JUS,  COE  f ( CM-1 A T M-l  I FROM  1675  T01735  CN-1 

UNO  3 

90 

NEW 

OAT  A H3/5, 003, 8.803, 14.  12, 19.  83, 2 3.  31, 23*  58, 2 3,  22, 21. 0 9,  26. 99, 25 

.8HNC3 

100 

NEW 

14 ,24.79,17.68,9 .420/ 

HNC3 

110 

NEW 

HABS=C. 

Hn03 

120 

NEW 

IF  IV.GE.05O. Q. ANO. V. IE.920. Q>  GO  TO  1300 

HNOJ 

130 

NEW 

If  tV,GE. 1275.0. AND ,V.LE. 1350-0)  GO  TO  1001 

HN03 

140 

NEW 

IF  IV. GE. 1675, 0. ANO. V .LF,  1735,0)  GO  To  1002 

HN03 

150 

NEW 

GO  TO  1003 

HN03 

lbO 

NEW 

1000 

IslV-045.  ) /5» 

HNO  3 

170 

NEW 

H AB  S=H1  (I  ) 

HNC? 

180 

NEW 

GO  TO  1303 

HN  03 

190 

NE  M 

1001 

i«nr-i2?o.)/5. 

HN03 

200 

NEW 

HABS*H?(I) 

HN03 

210 

NF  W 

Go  TO  1003 

HNO  3 

220 

NEW 

100Z 

I*( V-1670*) /5* 

HN03 

2 30 

NEW 

HflftS-rH?  (I  J 

H N 03 

?4y 

NEW 

1003 

RETURN 

HNC? 

250 

NE  W 

eno 

HNO  3 

260 

NE  w 

66 


ooDoon^noonoo 


_,,.x  TTv  fi'»AnTlCAM^' 

SSSSSS*  - 
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1 


3 


4 


5 

6 


SUBROUT INC  POINT  (X , Y N , N , N P , T X , I P) 

POIN 

10 

0 

1 

COMMON  7<  34)  ,P|7,  34)  ,T<7,34>  ,£H(U  ,34)  , WH(7,34)  ,N,NL  ,RE,CW,CO,PI 

PO  IN 

20 

NEW 

COMMON  /tn\/  IENTSStKHAX 

POIN 

30 

NEW 

DIMENSION  TX  Ml) 

pOIN 

4 0 

NtW 

*PCIN 

50 

B 

4 

subroutine  point  computes  the  mean  rierative  index  above 

AND  BELOHPCIN 

60 

a 

5 

A GIVEN  ALTIUIOF  AnO  INTERPOLATES  EXPONENTIALLY  TO  DETERMINE  TNt 

POIN 

70 

n 

b 

EQUIVALENT  ABSORBER  AMOUNTS  AT  THAT  ALTITUDE. 

POIN 

ao 

9 

7 

POIN 

90 

p 

6 

*PClN 

100 

0 

9 

PylN 

110 

e 

1 0 

X I S THE  HEIGHT  IN  QUESTION 

POIN 

120 

a 

11 

T X < <3 ) ANO  YN  ARE  THE  MEAN  REFRACTIVE  INDICES  ABOVE  ANC 

0FLOH  X 

POIN 

130 

0 

1 ? 

N IS  THE  LEVEL  INTEGER  CORRESPONDING  TO  V OR  THE  LFVEL 

BELOW  X 

PCIN 

t 40 

P 

13 

NP  = 1 IF  X COINCIDES  NlTH  MODEL  ATMOSPHERE  LEVEL  ,IF  NOT  NP  = 0 

POIN 

150 

B 

14 

TXCJ-6)  ARE  ABSORB  FR  AMOUNTS  PER  KM  AT  HEIGHT  X 

POIN 

if.0 

G 

15 

***>«***« 

•POIN 

170 

a 

16* 

N = NL 

POIN 

iau 

p 

17* 

NP=  0 

POIN 

iso 

a 

1 0 

IMX„LT.O#0>  X = Z(1) 

POIN 

20  0 

0 

1 9A 

IF  (X.GT.7(NL))  GO  TO  4 

PQ  JN 

210 

B 

1 90 

DO  1 1=1, NL 

POIN 

c?o 

B 

70 

N=I 

POIN 

230 

e 

21 

IF  <X-?<m  2,4,1 

PC  IN 

240 

8 

■3  O 
<-  l- 

CONTINUE 

F>0IN 

250 

B 

23- 

J?=N 

POIN 

260 

a 

:r> 

N=N-1 

POIN 

c 7 0 

13 

?f> 

FAC  = C X-  U N)  ) /<Z<  J2)  - Z<N)  ) 

POIN 

2 6 0 

R 

2 ? 

PX1=P (M ,N)*(P(M,J?) /P (M,N))*»EAC 

POIN 

29  0 

B 

?e 

TX1  = T(M,N)*(T(M,J?)7T(M,NI)*»EAC 

PCIN 

300 

B 

29 

HX1  = MH(  H,N>  MWHf H, J2  ) /HH(M  ,N) )**FAC 

PCIN 

3 10 

fl 

? 0 

T X ( 3) =CQ»PXl/TXl-4.56E-6*WXi*TXl*CH 

POIN 

320 

9 

31 

TXTt)-rO^PIMfjr)/T%NfJc)-4#Sf»E-o“Wr{(i‘tJc)rTiH,JC>pQW 

POIN 

330 

o 

3 2 

TX  <l)=CO*P  (H«N)/T  <M,N>-4.  5 6t>6*WH(M,  N>  *TIM,N)  *CH 

POIN 

340 

8 

33 

7T<9>  s0.5E-6*ITXf  2)  ♦ TX(7|) 

PCIN 

350 

B 

34 

YN=n.5C-6»  (TX  (1)  *TX  (3) ) 

POIN 

3 60 

B 

35 

IF  (IF.EQ.O)  GO  TO  9 

POIN 

3 70 

a 

3b 

00  3 Kc i , KM AX 

POl  N 

380 

NEW 

IF<K.FO.<U  GO  TO  3 

P 0 1 N 

390 

* 0 

3 78 

TX(K>=0.0 

PQl  N 

400 

3 7C 

IE  <EH(K,NI.EO.O.C)  GO  To  3 

P0lN 

4 1 C 

e 

3 8 

IF  (EHf KtNl.GT. 1U00 .0>  GO  TO  3 

POIN 

420 

8 

3^ 

TX(KI=cH(K,N)*(EH(K,J2>/EH(K,N))**FAC 

PCIN 

430 

B 

4 0 

CONTINUE 

POIN 

440 

e 

4 1 

GC  TO  9 

PO  IN 

450 

B 

47 

NP’i 

POIN 

460 

B 

43 

IF  (IP. EO.O)  GO  TO  6 

POIN 

4 70 

e 

44 

00  5 K= 1 » KM A X 

POIN 

480 

N£  H 

TX(K)=EH(K,N> 

POIN 

490 

9 

4 6 

TX(0)=EH(9,N) “1# 

PQIN 

500 

3 

47 

67 
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TN*0.0 

PlllN 

510 

p 

4 e 

CAROS  B 24  AND  50  T HPOUGH  59  ARE  NO  LONGER  REQUIRED 

POtN 

520 

n 

4 

IF  (N.GT.l)  yH'l'H(q,N-J>-1.0 

POIN 

530 

n 

49 

9 

continue 

PCI  N 

540 

5 

60 

IT  (IP.tO.n  PRINT  400,  X,N,NP,TX  (41  ,tN,IPf  ITX(K)  ,K=1,M 

potn 

550 

B 

6 1 

TX  <S)sTX<0)U# 

PQ  IN 

560 

R 

62 

YN*  YN*’  1 » 

P PIN 

570 

a 

63 

return 

POIN 

560 

B 

64 

c 

PCIN 

590 

a 

65 

400 

FORMAT  </,*  FKOH  POINTI  HEI GHT *• ♦ F 10 . 4, * K* , N' *, I 3 , * , NF* * 

, 12 ,* , REFPOIN 

COO 

R 

65 

1.  I NOE  * AROVc  ^ B^EPH  X = *,?f 11,4 ,* ,IP-* ,I3,/,12X ,*EQUIV. 

ABSORBER  POIN 

610 

B 

6 7 

?flHOUNTS  PER  K*  AT  X-*,fF10.3> 

POIN 

R 

66 

ENP 

POIN 

b3U 

R 

69 
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SURPOUT INF  ANGL  <H1  ,HZ  , A NGLt  , 01 ,LE N, Ht > 

ANGl  10 

C 

1* 

COMMON  7 < 3A1  ,P  t 7,  3*. ) , T ( 7,  m>  , EH<  11 ,3  41  , MH  < 7,3*1  , H,  NL  ,RE,CW  ,CO,Pl 

ANGl  20 

NFH 

COMMON  /F*U/  TFMTSS,KMAX 

ANGL  30 

NF  H 

DIMFNSION  TXIU) 

ANGL  40 

NL  H 

* AN  GL  VO 

0 

4 

ANGL  oO 

c 

V 

THIS  SUflPOUTUIE  CALCULATES  THE  I N I T I » L ZENITH  ANGLE  (ONE.LE) 

ANGl.  7 0 

c 

b 

7AK1NE.  It.Tn  ACCOUNT  REFRACTION  EFFECTS  GIVFN  HI, HZ,  ANO  BETA 

ANGL  80 

c 

7 

( WHrRF  RFI  A IS  THE  EARTH  CFNT  DE  ANCLE  SURTENOrO  «Y  HI  ANO  HZ  ), 

ANGL  90 

c 

b 

ASSUMT'T,  THF  EFFRACTIVF  INtlFX  TO  BE  CONSTANT  IN  A G I Vt  N LAYER. 

ANGl  100 

L 

9 

FOR  GREATER  ACCURACY  INCREASE  THE  NUMBER  OF  LEVELS  IN  THE  HOOEL 

ANGL  110 

c 

1 0 

ATMOSPHFPE • 

ANGt  120 

c 

1 1 

ANGl  130 

c 

12 

THIS  SUBROUTING  CAN  B F RflOVEI  FROM  T HF  PROGRAM  IF  NOT  PCQUIRtO* 

ANGt.  140 

c 

1 3 

‘•ANGL  100 

c 

1 4 

TPrBB 

ANGL  ItiO 

c 

IV 

Cft  = r-i/i  *o. 

ANGt  170 

c 

1 f' 

X 1=  PC  *H  1 

ANGL  180 

c 

1 7 

x 2 = r F ♦ h 2 

ANGL  ISO 

c 

1 8 

L FN=  0 • 

ANGL  200 

c 

1 9 

TT  - C 

ANGl  ?ld 

c 

2 C 

oi-ni*CA 

ANGL  220 

c 

2 l 

T ANG  = X2*SlN(0i>  / ( X2*C0S<  Bl)  -XI) 

ANGL  24U 

c 

2 2 

THFTrAT AN( TAnG) 

ANGl.  2 V0 

c 

2 3 

IF  (THE  T.LT.0.0  l THFI=THET*PI 

ANGL  PM 

c 

2 4 

SFM  I = sim  THF  T) 

ANGl  270 

c 

2 V 

A NC  = THF  T/Cfi 

ANGL  ?80 

c 

2 b 

PRINT  404  * 0 1 » ANG  , T ANO 

ANGl  ?99 

c 

2 7 

TN= THFT 

ANGl  300 

c 

2 8 

TMsTN-r  * r'  • C A 

ANGt  310 

c 

2 9 

A N G IF  - T HF  T 

ANGL  3?0 

c 

3 0 

f 0 T = 0 . 

AN  Gl  3 Tg 

c 

31 

GET  A- 0. 

ANGl  34U 

r 

2 

(TFT  1 = 0 

ANGl  3 VO 

c 

33 

n£T?*0 

ANGl  3 bv 

c 

34 

F3T  1=0 

ANGL  37U 

c 

7 V 

f n t ?*  o 

ANGL  380 

c 

3 b 

FBI  3=0.  0 

ANGL  390 

c 

37 

!F( B1 ,L f .0.01  GO  TO  ? 

ANGL  400 

n 

37* 

PRINT  vOO,  IT 

ANGL  410 

c 

38 

Vs?  ,*THET 

ANGl  4 ? 0 

C 

39 

IF  <Y  -FI  ,GT  . 1 .0  E-e  ) GO  TO  A 

ANGL  A31) 

c 

4 C 

IF  <tP#  FQ.l  00)  CO  TO  b 

ANGL  4 40 

c 

4 1 

XMlN=X?*COSlRII-&f 

ANGL  4 r> 0 

c 

42 

IF  (XHJN-H1J  8,  A, A 

ANGt  A bO 

c 

4 3 

1 HMlN=H? 

ANGL  470 

c 

4 4A 

H2  =H1 

ANGL  480 

c 

440 

HI =HMir 

ANGL  4% 

c 

4 4C 

I ANQLF  = C*  *>*PI 

ANGL  V00 

c 

4 UU 
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T Mf  T-ANGLt 
s r h i = i , o 
A NO* A NO U /C* 

PRINT  U Q ti  f \\  1 , ANG  , STM  I 

I r * loo 

CALI  POINT  (HI  j VN,  N,MP,  T X,  IPI 
J 1*N 

T X l- T X < 4) 

CALI  POINT  <4r,YH,NiNPfTXf IP» 

IF  wr.IO.l)  NrN  “ 1 

jr  - n 

If  (Jl.lO.J^I  TXl  sTtl  ♦YN-rH((l*  Jl) 

00  7 J- J1 , JO 
X 1 = 4T  *7  <J> 
x;»*RF  * ? ,JM* 

ir  (j.FG.ji)  xi  = f?r.*Hi 

U (J.FO.JO)  Xf  * PL  ♦ Hu 

!5ALP=Xl*Sf'HI  /xr 
AIP-ASI N(SAL  P) 

RN*IN(4, J*l)  / IH 14  * J) 

IF  ( ( J*  1)  .f  1.  JD  RN'YN/f  H(4,  J) 
n (J.IG.Ji)  WN*I  M(4, 

IT  ( ( J+l)  .FT. JC.ANO  ,J  ,FQ.  M > PN  = tN/fXl 
n n - t wr  t - a i ° 
f H»-T  AN  Ul»'> 

n tJ  .Nf.Jl)  r HsFlHT  AN  o HFT  ) 

f 01 Tf  0T  tFO 

ttt  T A » PF  T A t ML  T 

T H 1 r.  T HI  1 f C A 

RtaRtT/CA 

C - A l P»*C  A 

POINT  4 1)',  J,/(J>  ,THf  T,ALP#HrT,PriA.rOT  »r,'»THt»'U  tl. 

jr  { V."1,  f Oa  R>:  ♦ H ‘ > Or  FT- AL*1 
IF  (SAl 0,01  • ONI  KN-  1 • 

^PH1=SALP/KN 
THF  l«  A$IN<  !TPHI) 

•:0N1  \ Nt'F 

in  hi  .i  i , o. o*  r.o  to  om 
GC  TO  Of, 
fllNf  INUF 
r ANr.«-T  ANT. 

A NOLI  =Pl“ANr»L( 

1 N = ANGl " 

ANf.sAW  u /oa 

r^INT  4 04  » <M  * ANf.  ,TANi. 

IT  (HI, It  * 0*0)  GO  T 0 J 
CONTINWI 
I p- 10  1 

l ; 4 LI  POINT  f HI, YN,N,NP1 ,TX, IP) 

T X 1 = T X(  01 


! 
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YN1  = YN 

IV  <NP1.E0.1>  N=N**1 
J?*NL 

if  <n.eQ.r>  jr-^i 

J 1=  N 
J-Jl*  1 

IF  <H2.GF.H1I  GO  TO  13 
ULL  POINT  (M3.  YN.N.NF.  T X,  IP) 
tx;*y  x (9) 

YN2=YN 

J3«N 

IF  <JUEQ.  J.’l  TX?#YM*TXHI-fHt‘i,  Ji) 

J'  J~\ 

x i-p  t ♦:  ij*ii 

X?=»E*/<v1) 

IF  (J.l  0. J! > Xi  ♦ Hi 
IF  < J , F G«  J2)  X2  = RF>H2 
SAt  P*X1*SPHI /X? 

H»-IN= 

PRINT  40?.  J . XI , ? < JJ , 5PMI ,SALP» HHlN.Rr 
IF  (CALP.Lt* 1. 0)  GO  U)  11 
SAl P*SPHI 

IF  (HHIN.GT.H?)  GO  TO  18 
ALP  s AS  I N<  SAl  P) 

THF  Tj AS JH (SPHI) 

»M  T *ALP  «T  HF 1 
nr  Tl  * OF  Tl  ♦ Ol  T 
f p-  T A N ( Al  »*> 

IF  (J.Nf.JU  Fn=rp-1  AN<THfT  > 

FttT  i*FBTt  +F9 
IH1=1HFT/CA 
itt  =0F  T/CA 
AL  = AL  p/C  ft 

PRINT  4 0?  , J.X3.THM  . Al  P »‘U  11  • ‘TFT  , RH  IN  .HNIN,  MU  1 , t HI  tlU  , AL 
iF  <V  • FQ  • Rr  ♦ H2  i 0=  FT  -At  ° 

KM  F H<  9,  J) 
ir  (J  « I Q » J I ) RTF-YN1 
U U.fC.J?)  »FF-TXr 
Tc  < J.F  Q» 1 ) GO  TO  1 r 
R N * F H < SI  . JI/CH  <9,  J-l  ) 

IF  U.K.Jtl  HNaYNl  /F  H ( 9 . J-l) 

I F < J , L'O  . J ? ♦ 1 ) R N = W 1 f / T X T 
IFUJ.FO.J2i  RN=MfF/YN2 
IF  <9AlF\Gt.RN)  RN=1, 

SPHIa  SALP*»'N 

IF  <?<J) GO  TO  12 
(-0  TO  10 
> 1 a X ? 

u < APS  <7  u > -HD  a.  T .1  .OF  - 10  . AND  . J.Nt  « 1 ) tiU  TO  13 
GO  TO  14 


AnGL  10  10 

C 

94 

angl 10  20 

c 

9*> 

ANGl 1030 

c 

9f>A 

ANGl 1040 

c 

9b‘l 

ANGl  lU'tO 

c 

9 7 

ANGl 10*0 

c 

98 

ANGl  10  ro 

c 

99 

ANGl 13 80 

c 

14  0 

ANGl  1090 

c 

1 0 1 

ANGL  11  00 

c 

102 

ANGl 1110 

c 

1 0 3 

ANGl  1 1?Q 

c 

10  4 

ANGL 1130 

r 

iur. 

ANGl  1140 

0 

If  G 

ANGl 11*0 

r 

1 0 7 

ANGl  UNO 

c 

1 U 8 

ANGl  1170 

c* 

119 

A NGl  1 1 fl  0 

c 

1 1 0 

ANGl  1 190 

c 

11  l 

ANGl 1200 

c 

J i 3 

ANGL 1210 

c 

113 

ANGL  1 220 

c 

1 1 4 

ANGl 1 ?3C 

c 

1 1 4 

ANGl  1?4J 

c 

lit* 

ANGl  1?*>0 

c 

11  7 

A NGL  1 ?«»0 

c 

UA 

ANGL 1 ?70 

c 

1 1 9 

ANGL 1 280 

r 

1 3 0 

ANGl  1?M0 

c 

12  1 

ANGL  1 JOu 

c 

12  7 

ANGl 1310 

c 

13  3 

ANGl 1330 

c 

134 

ANGl  1 3 3C 

c 

1 ■“> 

ANGl 1340 

c 

1 2 1* 

ANGL 1 3 *0 

c 

137 

ANGl  13  GO 

c 

12  rt 

ANGl  137  0 

l- 

129 

ANGl  t 3 80 

c 

13  0 

ANGL  l 390 

c. 

131 

ANGl 1400 

f 

132 

ANGl t m10 

c 

l 3 3 A 

ANGl 1430 

r 

1 3 *!l 

ANGl  1<*30 

c 

1 3 3G 

ANGl  l 44(1 

c. 

1 3 4 

ANGl  14^0 

c 

1 3 '■ 

ANGl 1 4 GO 

c 

1 3 h 

ANGl  1470 

c 

1 3 7 

ANGl 1480 

c 

1 3 8 

ANGl 1490 

c 

I 39 

ANGL  1 *>  0 0 

c 

1 4 U 
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j=j-i 

X 1*RE*7  (j«i  > 

IF  (J.EQ.Jl)  KI  « RF  *H1 

IF  IJ.FO.  J7.  ANO.J.NF  . Jl>  X l *RE  fH.’ 

X?sRE«7  U> 

HHIN*»)»SPHI~RE 
IF  (HMIN.LF.  0. 0)  r.O  TO  79 
IF  ( 7{ Jl .IT.HMINI  00  TO  1 8 
REF«EH(9,J> 

if  (j.ra.ja  rff*yn 

S*LF«X1«SPHT/XJ 

«LPMSIMS«LP> 

THETsXSIN(SPHI) 

0 fc  T * ALF'  “TMF T 

F0*TAN<  »LP> -TAN < TUFT) 

FHT2iFBT?«Ft 
BET  r.flt  T Jf  BET 
BM I N*  OF  T 1 ♦ BE  T 7 
Al*M,o/CA 

thi  = thf  tecs 

PRINT  A02,  J,  X?,  THE  T , «L  P . BE  T7,  Bf  T ,IM1N  ,HMIN  ,FPT  2,TB1  ,8r,  AL 
RN.RF  F/FM  (O,  J-n 
IF  (SAlP.Gr.RN)  RN*1.0 
SPHI«SALP*RN 

00  TO  J.T 

TXJ«YN1«TX(9|-EH(9,  Jtt 
YNUTXJ 

IF  (ABF  <H?-7  ( J*l)  ) .Lt  ,1.  JE-5)  YN1»IX<9> 

IF  <»PS(H1-Z(J*1  1 | , ic.l.  Gi-S)  YMI  = T X < 9> 
f’N  = l . 0 
1.0  TO  1« 

call  point  <hnin» yn,n,np, tx, IP) 

IP'102 

TX.5«T  X(4I 

IF  (J.f  O.JI.ANn.HJ.OE.HlI  CO  TO  17 
IF  (J.EQ.Jl. OR. J.FQ.J7I  T\X»YNJ.TX(9>-tH(4,J) 

IF  (HMIN.GT. H?>  TX.YiTXCBI 

IF  (J.f  q.  Jl.  AMO.  HMI  N.  r.T.  B2  ) r.O  TO  17 

RN-REF/  TX  S 

IF  (S9LF.GE.RN)  R N-'  1 . 

SPH1. SALP’RN 

X«X1»SFMI-pf 
niFiA<«s  (HMIN-X) 

H HI  Ns  X 

IF  (OIF-1. (If  >51  19,19,10 
X?*RFFhPIN 

PRINT  903,  BNIN.OIF.RN 

THE  T« AE INC  SPHII 

If  (RN. Ft).  1.0)  FBT,T«-TAN(TnET> 

IF  CRN.  £Q,  l>  GO  T(1  20 


ANGL 1010 
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14  1 

ANGL  162U 

c 

14? 

ANGl 1633 

c 

1 H 3 

An  GL  1640 

c 

144 

AnGl 1560 

c 

146 

AnGc  1661) 

c 

146 

AnGlIH^o 

c 

14  7 

angl l^ao 

c 

14  6 

anclic jo 

c 

149 

ANGl  16  00 

c 

16  0 

A nG  L .1  ti  1 1> 

c 

15  1 

A nGl 1 6 ?0 

c 

10? 

ANGl 1630 

c 

103 

AnGL 1640 

r 

104 

ANGl  ltj5u 

c 

100 

AnGI 1660 

c 

1 0 6 

ANGL 1670 

c 

15/ 

ANGL 1660 

<7 

1 06 

ANGl 1690 

r 

109 

ANGL 1700 

G 

160 

AN GL 1710 

c 

lb  1 

ANGL  17  20 

c 

1 b? 

ANGL  1730 

c 

163 

ANGL 1740 

c 

164 

ANGl  1750 

c 

l 66 

ANGL I7b0 

c 

166 

ANGL  1770 

c 

167 

ANGL  17  81) 

, c 

1 bft 

ANGL 1 7 SO 

c. 

1 69 

ANGLloOO 

c. 

1/0 

ANGL  Irt  10 

c 

1 7 1 

ANGL  1870 

c 

1 7 r 

ANGL  t *U 3C 

c 

1/3 

ANGL  1640 

c 

174 

ANGL 1350 

0 

1 70 

ANGL l o 6 0 

c 

1 7b 

ANGL  1 ft/- J 

c 

1 7 7 

ANGL  I860 

c 

17ft 

ANGL  169  0 

c 

1 79 

ANGL 1 900 

c 

1 6 0 

ANGL  1910 

c 

1 fll 

ANuL 1920 

c 

18? 

ANGl  19  30 

c 

163 

A NG  L 19  4 u 

c 

164 

ANGl 1 950 

c 

186 

ANGL I960 

c 

166 

ANGl  19  70 

c 

167 

ANGL 1980 

c 

166 

ANGl  19  90 

c 

1 ftft'l 

ANGL  7 0 00 

c 

169 

i ' <*t  A.SH « 
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nNX=(T>3-l.Q)#AL0G(<T  X3-1.  09 /(  REF-1.  0)  ) /(  X,?- XU 

A NGL  20  1 0 

C 

10  0 

FRT3s-TAN<THET>Ml.0-l.fl/(1.0+TXJ/  (X2*nNX))  ) 

ANGL2030 

c 

101 

Of  Trfl  ,*-»rl-THFT 

ANGL 20 30 

c 

is? 

'JET?*BfT2»PfT 

angi204q 

c 

103 

9M I N*  Rf  T 1 fRCT? 

ANGL20G0 

c 

)04 

If  <H?,GF.M1)  r.O  To  2 3 

ANGL2060 

c 

105 

BET=n€Tl*2.*nET? 

ANGL20 7u 

c 
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1.8b  1.61  1.1?  1.3?  1.  52  1.  70  1.90  2.01  1 . 

2.11  2.28  2.21  2 .13  2 . 00  1.91  1 ,92  1.9  7 1. 

1.  58  1,  2 7 1.  20  1.1  8 1.11  0.99  4.86  0.71  0. 
-0.35-0.99-0. 69-0.  7 9-0.  94-  1.  11-1.24-1.41-1. 
-2.  81  -2  .9  9-3. 16-3 . 37  - 3 . 56-3 . 7 5 -2 . 9 4 -4 . 1 3- 4. 
-5.  00-5.  0 0-5. 00-5.00  -5 ,00  -5 .00-5 . 0 6-5 ,00-4. 
-2.  89-2. 79-2.  74- 2. 6 3- 2.  4 7- 2.  2 9- 2. 20-2. 1 7 -2  . 
-2, 07 -1.96-1. 77-1. 70-1. 63-1. 60-1.59- 1.43-1. 
-0.  93-3. 8 7-0  . 91  -0  .88-0. 71-0. 62 -0.59  -0.58-0  . 
-u.  01-0.  0 8-0. 20-  0.  16-  0.  02  0.  18  C.32  0 .4?  0. 
0.  58  0.5  3 0 . 44  0.39  0.  36  0.35  0.23  0.26  0. 

0. 32  0,37  0.58  0.64  0.87  0,98  1.00  1.02  1. 

1.  40  1.3  2 1 . 72  1.37  1,  42  1.50  1.4  2 1.38  1 . 

1,68  1.60  1.56  1.56  1.63  1.64  1.56  1.49  1. 
1.62  1.6  3 1 . 71  1.72  1.  70  1.70  1.67  1.62  1. 

1.  26  1.2  0 1.  13  1 . 14  1,19  1.29  1.50  1.72  1 . 
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-5.03-5.0  0-5.00-5.30-5.00  -5.00  -5.00  -5  .00-5  . 
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93  1.91  2.12  2.10  2.01  2.18  1,99  8625 

88  1.91  1.91  1 . a2  1.93  1.74  L.61  89flu 

6 0 0.4  4 0.31  0.19  0.03-  0 , 07-0.  21  8975 

57-1  .73-1 .91  -?  .09-2  .27-2.45  -2 ,63  905d 
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1 3 1 . 08  1 . 08  1.16  1.16  1,  70  1.411  0 32  6 

3n  1.  38  1.49  1 , b 3 1.62  1,62  1.7ul  640l) 
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66  1.78  1.82  1.88  l.o?  1.89  1.9910d25 
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59  4*53  0.69  0.67  0*48  0,52  j, 6211 0 75 
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1 6-  2 . 3 6-  2 . 4 4- 2 . 4C  - 2. 4 9- 2 . 4 8 -2 . 4 3i  1 3 C ; S 
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49 -1.25-1. 24-1. 08-C. 90-1. u 6-0. 9 1118? 6 
38-J  .4  2-0  .48-0. 34  -Q.5  7-13.  17-0.2811900 
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51  0.6U  0.51  0.51  u.40  j.  44  0.4J122o0 

31-0. 41-0. 41-0. 79 -0.5 3-0. 69-0. 8412 2 76 
67-1. 5 1-1 , 66-1. 60 -1.69-1. 63 -1.5 11 2 350 

39- 1.33-1.40-1.35-1,37-1.39-1.4112426 

4 1-  1. 31-1. 15-1. 13-1. 20-1. 41-1. 8812  6 00 

57-1.69-1 . 70-1 . 70-1.66-1 . 84-1 .601 ?5> 6 
53-1, 48-1 . 48-1. 60-1.58-1 . 69-1 . 7912650 
09-3.  4 4-3,  79-4.  14-4,  49-4. 84-5. 191 27 ?5 
99 -2. 14-2. 41-2. 12-1. 99-  1.84-1. 7 91  3400 

22 - 0.91-0. 90-1. 31- 0.76-0. 90-0. 9013i<75 
85- 0.61-0. 61-0. 48 - 0.61-4  .92 -6  .6  31  3 55  0 

0.19-0.10  0,20  0.20  J.02  .2013625 

16  0.26  0.  0.09  0.09  0.09  0.0713?04 

26  0.07  0.09-0.29-6,69-3.69-0.7413776 
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| Appendix  B 

! Basic  Flow  Chart  for  LOWTRAN  4 

I 

i 

J A general  flow  chart  for  l,OWTl(AN  -1  is  given  in  Figure  li  1 which  shows  the 

I overall  mode  of  operation  of  the  program.  IVlo"e  detailed  flow  charts  are  also  given 

! for  the  main  blocks  in  the  program,  that  is,  where  the  equivalent  absorber  amounts 

and  refraction  calculations  are  made  (Figure  Kill,  the  transmittance/ radium  o loop 
(Figure  lid),  and  the  transmittance  calculations  (Figure  H4). 

The  notation  used  in  the  flow  charts  is  as  follows; 

1 (1)  If  a condition  stated  within  a given  block  is  fulfilled,  then  the  direction  of 

] flow  is  sideways  as  indicated  by  the  direction  m which  the  block  points  (for 

^ example  for  the  following  block  O). 

| (d)  If  the  condition  stated  within  a block  is  not  fulfilled,  the  flow  is  downwards. 

1 The  numbers  appearing  on  the  flow  charts  correspond  to  the  statement  numbers 

! given  in  the  main  program  (see  Table  Al). 


.■31 

1 


t>7 


a 

1 


1 


Appendix  C 

Symbols  and  Definitions 


AB 

AHAZE 

AHZ2 

AJ 

ALAM 

ALP 

ANGLE 

AO 

AVW 

B 

BBG 

BBIK 

BET 

BETA 

EJ 

CA 

CO 

cw 

Cl 

C2 

C3 


Absorption  at  frequency  V;  al$o  average  transmittance 
Aerosol  number  density  for  MODEL  - 7 
Aerosol  number  density  for  MODEL  = 7 
Equivalent  absorber  amount  per  km  at  level  J 
Wavelength  (pm) 

Angle  of  arrival  at  adjacent  level 

Input  zenith  angle  (degrees) 

Constant  A = (R  + Hl)n  sin  6 
o o o 

Average  wavelength  used  in  refractive  index  expression 
Storage  parameter  for  BETA 

Black  body  function  of  boundary  times  the  total  transmittance 
along  the  path 

Black  body  function  of  the  IK  layer  and  the  frequency  V 

Angle  subtended  at  the  earth's  center  a|;  path  traverses  adjacent  levels 

Total  angle  subtended  by  path  at  earth's  center 

Equivalent  absorber  amount  per  km  at  level  J + 1 

Conversion  factor  from  degrees  to  radians 

Wavelength  dependent  coefficient  used  in  refractive  index  expression 
Wavelength  dependent  coefficient  used  in  refractive  mdex  expression 
Log  absorption  coefficient  for  water  vapor 
Log  absorption  coefficient  for  uniformly  mixed  gases 
Log  absorption  coefficient  for  ozone 
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04 

05 
(,'b 
07 
OVA 

oa 

n 

DTAU 

DP 
].)S 
L)Y 
l>Z 
E(K) 
EH(1, 1) 
E11(2,I> 
EH  3.1) 
EH(4, 1) 
EH(5,  it 

EUU>,I) 

E1I(7, 1) 

Ell(fl.I) 

eu<o,  i) 

EH(10,1) 

EH(  1 1,1) 

KV 

F 

PP 

PAO 

PACTOK 

P’O 

PW 


H 

HI 

112 

HA  US 


Absorption  eueffi  eient  for  nitrogen  i~4  (in;) 

Absorption  coefficient  for  writer  vapor  continuum  (~  10  pm) 

Extinction  coefficient  for  molecular  scattering 
Extinction  coefficient  for  aerosol  models 
Aerosol  absorption  coefficient 

Absorption  coefficient  for  ozone  (UY  and  visible  regions) 

Water  vapor  amount  (pr.  em/km)  at  level  1 

Differential  transmittance  (due  to  absorption)  across  the 
IK  layer 

Dew  point  temperature  (°C) 

Path  length  from  level  1 to  level  I + 1 

Wavenumber  increment  at  which  transmittance  is  calculated 
Height  increment  from  level  1 to  level  1+1 
Equivalent  absorber  amounts  per  km  at  height  HI 
Equivalent  absorber  amount  per  km  for  11. ,0  at  level  Z(l) 

Equivalent  absorber  amount  per  km  for  l'0,+  N,;P>  etc  at  level  Z(l) 
Equivalent  absorber  amount  per  km  for  at  level  Z(l) 

Equivalent  absorber  amount  per  km  for  N9  at  level  Z(l) 

Equivalent  absorber  amount  per  km  for  H„Q  continuum  at 
leveled),  (10pm) 

Equivalent  absorber  amount  per  km  for  molecular  scattering 
at  level  Z(I) 

Equivalent  absorber  amount  per  km  for  aerosol  extinction  at 
level  Z(I) 

Equivalent  absorber  amount  per  km  for  ozone  (UV  and  visible) 
at  level  Z(l) 

Mean  refractive  index  of  layer  above  level  Z(l) 

Equivalent  absorber  amount  per  km  for  11,0  continuum  at 
level  Z(l),  (4  pm) 

Equivalent  absorber  amount  per  km  for  nitric  acid  at  level  7.(1) 
Integrated  absorber  amount  from  level  i to  level  1+1 
P unction  for  determining  saturation  vapor  density  of  water  (gm  m ‘^) 
Black  body  function  ( W/cm“-ster-pm) 

Interpolation  parameter 
Integration  weighting  parameter 

Transmission  function  logarithmic  absorber  amount  scale  for  O 

O 

Transmission  function  logarithmic  absorber  amount  scale  for  11,  0 
and  the  uniformly  mixed  gases 

Altitude  dependent  control  parameter 
Initial  altitude  (km) 

Final  altitude  (km) 

Nitric  acid  absorption  coefficient 
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117(1) 

117-1 

117:2 
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1ATM 

ICOUNT 

inv 

IEMISS 


IF1N1) 

I HAZK 

IJ 

IK 

IKI-O 
IK  MAX 
II. 


IM 

IP 


IT  Y PK 
IV 
IV 1 
IV  2 
IX  Y 
J 

J EXTRA 

J MIN 
JP 

jSTOR 

J 1 

J 2 
Iv 


Aerosol  number  density  (no.  rm  ’ ) 

Estimated  tangent  height  (km) 

Minimum  altitude  of  path  trajectory  (km) 

Nitric  acid  volume  mixing  ratio  (times  1.0  !• -i  OP)  at  the  level  /(U 
Interpolated  nitric  aeid  volume  mixint;  ratios 
Hollerith  titles  for  visibUR 

Aerosol  number  density  (no.  eni  "')  for  211  lem  visual  range 
Aerosol  number  density  (no,  em  "’)  for  f>  km  visual  ranee 

Running  integer  used  as  altitude  (level)  indicator  and  frequency 
indicator 

Number  of  model  atmospheres 
Output  page  counter 
Frequency  increment  (cm  ) 

Input  control  parameter  determining  modi'  of  program  execution 
(=0  for  transmittance,  =1  for  radiance  mode) 

Indicator  for  using  subroutine  ANG1, 

Aerosol  model  indicator 
Equals  IK 

Running  integer  used  as  layer  indicator  along  the  atmospheric  path 
I-ower  limit  of  layer  loop  (-1) 

Upper  limit  of  layer  loop 

Integer  indicator  used  to  determine  if  Hie  atmospheric  path 
intersects  the  earth 

* * 

Parameter  used  when  reading  in  a new  atmospheric  model 

Indicator  for  using  subroutine  I'OIiNT  to  calculate  refractive  index 
Only  (IP  = 0)  or  equivalent  absorber  amounts  also  (IP  i 0) 

Indicator  for  type  of  atmospheric  path 

Frequency  at  which  transmittance  is  calculated 

Starting  frequency  (equivalent  to  VI) 

Last  frequency  (equivalent  to  V2) 

Parameter  for  terminating  program  and  cycling  indicator 
Running  integer  for  altitude  identification 

Integer  indicator  used  when  111,  112,  and  HMIN  are  in  the  same 
layer  (ITYP1C=2) 

Altitude  indicator  for  minimum  height  of  path 
Print  option  parameter 

Integer  indicator  used  when  vertical  profile  changes  from 
downward  to  upward  path 

Level  indicator  for  altitude  111. 

Level  indicator  for  altitude  H2 

Absorber  indicator,  K - 1,  2,  3,  etc.  , corresponds  to  H.,0 
uniformly  mixed  gases,  0„  etc,  respectively 
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IU 

Ml, 

MODEL 

Ml 

M2 

M3 
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Nil 

Nl, 

NLL 

NP 

NP1 
NP2 
P(M,  l) 
PI! 

PHI 

PI 

l'PW 

PS 

PSI 

PT 

KADMAX 
HA  DM  IN 
HADSUM 
RANGE 
UK 


Upper  limit  of  absorber  ;mmunt  loops  (-11) 

Integer  used  in  reading  two  model  atmospheres  on  one  eard 

Integer  used  in  reading  two  model  atmospheres  on  one  eard 
and  cycling  parameter  tor  downward  looking  paths 

Kroquenev  indieator  lor  nitrogen  continuum  transmittance  calculation 
Frequency  indicator  for  ozone  transmittance  calculation 
Parameter  used  for  defining  longest  of  two  paths 

Integer  storage  for  parameter  LION,  needed  for  cases  run  in 
succession 

Frequency  identifier  for  I'V  and  visible  ozone  transmittance 
calculation 

Frequency  identifier  for  UV  and  visible  ozone  transmittance 
calculation 

Integer  used  to  identify  required  model  atmosphere 
Number  of  levels  in  radiosonde  data  input  (MODEL  - 7) 

Integer  used  to  identify  required  model  atmosphere 
Integer  for  selecting  II,, O altitude  profile  for  (M=M1) 

Integer  for  selecting  temperature  altitude  profile  for  (M_M2) 

Integer  for  selecting  O.^  altitude  profile  for  (M=M3> 

Indicator  for  level  below  given  input  altitude  used  in  POINT 
subroutine;  also  as  frequency  indieator  in  UV  and  visible 
ozone  transmittance  calculation 

Frequency  indicator  for  water  vapor  continuum  transmittance 
calculation 

Number  of  levels  in  model  atmosphere  data 
Equals  NL-1 

Indieator  for  determining  whether  HI  or  112  coincide  with  levels 
in  the  model  atmosphere 

Value  of  NP  for  altitude  HI 

Value  of  NP  for  altitude  H2 

Pressure  (mb)  at  level  1 for  model  atmosphere  M 
180°  - PHI 

Angle  of  arrival  at  112 
3.  Hl?92t>54  that  is  (ji) 

Partial  pressure  of  water  vapor  (in  atmospheres) 

Total  pressure  in  atmospheres 

Angular  deviation  of  path  from  initial  direction 

Product  of  total  pressure  (atm)  and  the  square  root  of  273 /T(M,  l) 

Maximum  value  of  radiance 

Minimum  value  of  radiance 

o 

Integrated  radiance  (W  / cm"  -ster ) 

Path  length  (km) 

Earth  radius  (km) 


illi 


PS 


REK 

Rll 

UN 

HO 

HX 

HI 

SALP 

SPH1 

SR 

SUM 

SUM  A 

SUMV 

SUMVV 

T(M,  1) 

TAUC. 

THEY  (IK) 

TBOUNU 

THET 

THETA 

TMP 

TR 

TS 

TSNKW 


TSOLD 
TS  1 
XT 

TX(K) 


TX(9) 

TX(10) 

TX1 

TX2 

TX3 

Tl 

T2 


V 

Vll(K) 

VIS 


Refractive  index  of  air  at  level  I 
Relative  humidity  (%) 

Ratio  of  refractive  indices  of  air  above  and  below  a given  level 
Karth  radius  (km)  read  in  as  input  ( = Rlf) 

Ratio  of  earth  center  distances  between  adjacent  levels 

The  product  of  the  sine  of  the  initial  venith  angle  and  the  earth 
center  distance  to  starting  altitude 

Sine  of  angle  of  arrival  at  adjacent  level  (cf  sin  a) 

Sine  of  the  local  7enit.l1  angle  at  a given  level  (cf  sin  (1) 

Slant  range  (km) 

Sum  of  the  optical  thicknesses  of  absorbers  4 through  1 1 
Accumulated  integrated  absorption 
Radiance  (W/cm  -ster-cm  ) 

Radiance  (W/cin“-ster-|i  nil 

Temperature  (°K)  for  model  atmosphere  M at  level  1 
Total  transmittance  to  the  boundary 
Average  temperature  of  the  IK  layer 
Input  temperature  of  the  boundary  ir.  °K 
Zenith  angle  at  a given  level  (in  radians) 

Zenith  angle  at  a given  level  (in  degrees) 

Ambient  temperature  ("Cl 

Transmittance  scales  for  transmission  functions 

Ratio  of  standard  temperature  (273.  0°K)  to  temperature  at  Level  1 

Transmittance  (due  to  scattering)  to  the  far  boundary  of  the  IK  layer 

Transmittance  (due  to  scattering)  to  the  near  boundary  of  the  iK  layer 
Ratio  of  296  °K  to  ambient  temperature  (°K) 

Ratio  273.  15/(TMP+273.  15) 

Equivalent  absorber  amounts  per  kin  at  a given  altitude  obtained 
from  POINT;  also  transmittance  values  at  a given  wavelength  for 
each  absorber  type  (K  = 1,  KIVlAX) 

Total  transmittance  at  frequency  V 
Absorption  due  to  aerosol  only  at  frequency  V 
Refractive  index  of  layer  above  initial  altitude  HI 
Refractive  index  of  layer  above  final  altitude  112 
Refractive  index  of  layer  above  minimum  altitude  HMIN 
Temperature  of  the  boundary  (°K) 

Temperature  of  the  1+1  boundary  used  in  index  of  refraction 
calculation 

Running  frequency  (cm  *) 

Integral  of  the  equivalent  absorber  amounts  from  111  to  level  1 
Visual  range  (km)  at  sea  level 
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V UMAX 

V R MIN 
VX 

VI 

V2 

W(K) 

\\  H(M,  I) 
WLAY  (I,  K) 

WO(M,  i) 
WPATHdK,  K) 

WS1 

WS2 

WS3 

W2 

X 

XD 

XH 

XI 
XX 

XI 

X2 

V 

YN 

YN1 

YN2 

YY 

Z(I) 


Frequency  of  the  maximum  radiance  (cm  O 

Frequency  of  ihe  minimum  radiance  (cm  ) 

Wavelength  at  which  aerosol  coefficients  are  read  in  (jam) 

Initial  frequency  for  transmittance  calculation,  cm 

Final  frequency  for  transmittance  calculation,  cm  1 

Total  equivalent  absorber  amount  for  entire  path 

Water  vapor  density  for  atmospheric  mcdel  M at  level  I (gm  m "’) 

The  absorber  amount  for  the  species,  K,  and  the  atmospheric 
layer,  I 

Ozone  density  for  atmospheric  model  M at  level  I (gm  m 3) 

The  cumulative  absorber  amount  of  the  species,  K,  for  the  IK  layer 
along  the  atmospheric  slant  path 

Transmission  function  scaling  factor  for  HgO  at  given  wavelength 

Transmission  function  scaling  factor  for  C02>  etc.  , at  given 
wavelength 

Transmission  function  scaling  factor  for  0~  at  given  wavelength 

J -3 

Water  vapor  density  for  atmospheric  model  M at  level  1+  1 (gm  m ) 

Input  height  to  POINT  subroutine 

Wavenumber  interpolation  parameter  in  UV  ozone  transmittance 
calculation 

Wavenumber  interpolation  parameters  in  H20  continuum  calculation 
Wavenumber  interpolation  parameter 

Wavenumber  identification  parameter  for  UV  ozone  transmittance 
calculation 

Earth  center  distance  of  level  I 
Earth  center  distance  of  level  I + 1 
Input  zenith  angle  in  radians 

Refractive  index  of  layer  below  input  height  from  POINT  subroutine 
Refractive  index  of  layer  below  initial  altitude  HI 
Refractive  index  of  layer  below  final  altitude  H2 
Aerosol  absorption  coefficient  at  frequency  V 
Altitude  at  level  I in  km 
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Appendix  D 


Errata  Sheet  No.  2 (September  ]977),  Atmospheric  Transmittance 
From  0.25  to  28.5  jim:  Supplement  LOWTRAN  3B(1976), 
AFGL-TR-0158,  1 November  1976, 
Environmental  Research  Papers,  No.  587 

1.  Page  52  - Line  A 3M  should  read  A 3:- 

2.  Page  53  - The  second  line  A 81  should  be  removed 

3.  Page  54  - Reverse  the  order  of  statements  A 103F  and  A 103G  and  relabel  them: 

l.  e.  , IF  <M.  EQ.  0)  Z(K)  = HI  A*  103F 

J ^ I FIX  (Z(K)  + 1.  0E-6)  + 1 A*  103G 

4.  Page  56  - Line  *A  185B  should  read: 

EH(5, 1)  = D-'PPW* EXP  (6.  08*(TS1-1.  0))  + 0.  002  ■ D(PS-PPW ) *A  18511 

5.  Page  56  - Line  A 204  is  correct,  i.  e.  IP=-1 

6.  Page  63  - Line  A 586A  should  be  removed. 

7.  Page  64  - Replace  UALENTVENT  by  EQUIVALENT  in  line  A 623 

8.  Page  78  - The  9th  and  10th  cards  from  the  bottom  of  page  7 8 should  be  inter- 
changed. The  wavenumber  identifications  for  these  cards  are  17800  and  19400 
(see  extreme  right  hand  side  of  card). 

9.  Cautionary  Note:  When  standard  radiosonde  data  arc  used  (MODEL-7  option), 
insert  a card  for  sea  level  even  though  the  required  '■.ransmittancc  path  does  not 
extend  to  sea  level.  The  reason  for  inserting  the  sea  level  altitude  card  is  tc 
correctly  interpolate  the  aerosol  number  densities  at  the  required  altitudes  for  a 
given  sea  level  visual  range.  However  the  above  docs  not  apply  if  the  user  is  in- 
serting his  o.vn  aerosol  extinction  data  for  a given  starting  altitude. 
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10.  Example:  On  page  7 9 a set  of  input  data  for  LOWTRAN  313  is  given.  The  first 
example  (represented  by  the  first  three  cards)  is  to  calculate  the  transmittance  for 
a 65°  zenith  angle  slant  path  from  altitudes  2.  5 km  to  3.  5 km  for  a 23  km  visual 
range  (rurM  aerosol)  subarctic  winter  atmosphere  covering  the  wavenumber  range 
from  2350  to  2450  cm  ' (i.  e.  4.  08  - 4.  27 Jim). 

The  four  examples  which  follow  are  to  calculate  the  average  transmittance 
from  1820  to  20,000  cm  * (i.  e.  from  0.5  to  5.5  ,um)  for  a 10  km  horizontal  path 
at  sea  level  (using  the  1962  U.  S.  Standard  Atmosphere)  for  four  different  aerosol 
models,  namely,  Maritime.  Urban,  Rural,  rr.d  Tropospheric  respectively. 

it  is  recommended  that  the  various  aerosol  model  data  sets  be  labelled  and 
stored  for  further  use. 

Tlie  extinction  coefficients  for  the  Average  Continental  aerosol  model  originally 
contained  in  LOWTRAN  3 are  included  as  an  additional  data  set  at  the  end  of 
LOWTRAN  3B  (1976)  card  deck  available  from  the  National  Climatic  Center, 

Federal  Building,  Asheville,  No.  Carolina  28801  for  a charge  of  $20.  00.  (Please 
address  requests  to  Mr.  R.  Davis.) 

11.  Page  25  - delete  the  last  three  lines  on  this  page. 

12.  Page  26  - Delete  lines  +A  494  and  *A  495,  and  change  line  ''  A 491  to  read: 

NH  = XI  + 1.  001 

13.  Page  33  - Lines  A 85G  and  A 8511  should  be  interchanged  to  be  consistent  with 
page  53. 

14.  Page  54  - Insert  - 1FUXY.GT.3)  GO  TO  8 after  A 104N. 

15.  Page  62  - A 563  should  be 

IFfVLS.GT.  0.  0.  AND.  VIS.  LT.  2.  0)  XX=0.  158 

Cautionary  Note:  The  temporary  fog  correction  will  scale  aerosol 
extinction  regardless  of  atmospheric  path.  As  presently  coded,  it  probably  should 
be  restricted  to  horizontal  paths  (1TYPE=1)  under  1 km  in  altitude. 

16.  Page  66  - Delete  card  C 21  B in  ANGL.  ^ 


*NEW  ERRATUM 


